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Abstract 



In this paper, the low-mass electron pair production in 158 AGeV/c Pb-Au col- 
lisions is investigated with the Cherenkov Ring Electron Spectrometer (CERES) 
at the Super Proton Synchrotron accelerator (SPS) at CERN. The main goal is 
to search for modifications of hadron properties in hot and dense nuclear mat- 
ter. The presented re-analysis of the 1996 data set is focused on a detailed 
study of the combinatorial-background subtraction by means of the mixed-event 
technique. The results confirm previous findings of CERES. The dielectron pro- 
duction in the mass range of 0.25 < m cc < 2 GeV/c 2 is enhanced by a factor of 
3.0 ± 1.3(stat.) ± 1.2(syst.) over the expectation from neutral meson decays. The 
data is compared to transport model calculations and seem to favor the version 
including in-medium effects. 

Furthermore, the development of a new technology to manufacture ultra- 
lightweight mirrors for Ring Imaging Cherenkov detectors (RICH) is described. 
Replacement of the RICH-2 glass mirror by a mirror almost transparent to elec- 
trons would considerably improve the performance of the upgraded CERES de- 
tector system including a radial Time Projection Chamber (TPC). 

Zusammenfassung 

In dieser Arbeit wird die Produktion von Elektronenpaaren kleiner Masse in 
Pb-Au Kollisionen bei 158 AGeV/c mit Hilfe des Cherenkov Ring Elektron 
Spektrometers (CERES) am Super Proton Synchrotron Beschleuniger (SPS) des 
CERN untersucht. Ziel des Experimentes ist der Nachweis von Veranderungen der 
Eigenschaften von Hadronen in einem dichten und heiBen Medium aus Kernma- 
terie. Der Schwerpunkt der hier prasentierten Neuanalyse des Datensets von 1996 
ist die detailierte Untersuchung der kombinatorischen Untergrundsubtraktion 
mit der Methode der gemischten Ereignisse. Die Ergebnisse der Untersuchung 
bestatigen die vorhergehenden Resultate von CERES. Die im Massenbereich 
0.25 < m ee < 2 GeV/c 2 gemessene Rate von Elektronenpaaren ubersteigt die Er- 
wartung von den Zerfallen neutraler Mesonen um einen Faktor 3.0 ± 1.3(stat.) ± 
1.2(syst.). Die Daten werden mit theoretischen Vorhersagen im Rahmen eines 
Transportmodells verglichen. Die beste Ubereinstimmung ergibt sich mit einem 
Szenario, welches Mediumeffekte einschlieBt. 

In einem weiteren Teil der Arbeit wird die Entwicklung einer neuen Tech- 
nologie zur Herstellung ultra-diinner Spiegel fur Ringabbildende Cherenkov De- 
tektoren (RICH) beschrieben. Der Austausch des gegenwartigen RICH-2 Glas- 
spiegels mit einem fur Elektronen fast transparenten Spiegels wiirde die Leistung 
des mit einer TPC nachgerusteten CERES Detektorsystems betrachtlich steigern. 
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Chapter 1 
Preface 



This thesis is about the study of dielectron production in Pb-Au collisions at 
158GeV/c per nucleon measured with the Cherenkov Ring Electron Spectrome- 
ter (CERES) at the Super Proton Synchrotron (SPS) accelerator at CERN. Part 
of the heavy ion research program at CERN, CERES is committed to the explo- 
ration of nuclear matter under extreme conditions of high temperature and high 
density 

The interdisciplinary field of ultrarelativistic heavy ion collisions combines 
the elementary interaction aspect of high-energy particle physics with the macro- 
scopic matter aspects of nuclear physics. It is focused on the investigation of 
the properties of nuclear bulk matter made up of strongly interacting particles, 
i.e. hadrons, quarks, and gluons. The prediction is that nuclear matter would 
undergo a phase transition to a quark-gluon plasma (QGP) Ij, a gas of freely 
moving quarks and gluons. This triggered not only a strong theoretical interest 
in this field but also initiated a huge experimental effort to verify the prediction. 

First evidence for this new state of matter was found at CERN in the begin- 
ning of the year 2000 |2j. CERES was one of the experiments contributing to 
this fundamental discovery. 

During the phase transition, the quarks are expected to lose their constituent 
mass which leads to the restoration of chiral symmetry. Left- and right-handed 
quarks decouple and hadronic states of opposite parity become degenerate. 

The study of deconfined or chiral matter is not only relevant for the under- 
standing of heavy ion collisions but also for astrophysics and cosmology. The 
environmental conditions of a 158 AGeV/c Pb-Au collision at the SPS acceler- 
ator resemble those encountered in the evolution of the early universe, where a 
few tens of microseconds after the big bang a transient stage of strongly inter- 
acting matter persisted at temperatures of about 10 12 K and low baryon density. 
Another extreme of high densities and low temperatures created in heavy ion 
collisions at the SIS accelerator is close to the conditions occurring in the interior 
of neutron stars [13111113, where mass densities are likely to exceed 10 18 kg/m 3 - 
about four times the density of nuclei. 



1 



Chapter 2 

Relativistic heavy ion physics 



2.1 Hot and dense nuclear matter 

The observations of particle production in heavy ion collisions are related to the 
evolution of hot and dense matter - a general question of relativistic heavy ion 
physics - and in particular to the transition of hadronic matter to a quark-gluon 
plasma and the simultaneous restoration of chiral symmetry. 

First, the fundamental physics properties to describe a strongly interacting 
system of nucleons in vacuum shall be introduced. Next, the modifications in 
the presence of a hot and dense medium are discussed with special emphasis 
on theoretical concepts applied to explain experimental data. Finally, dielectron 
production is related to the properties of hadronic processes occurring in nuclear 
collisions. In this section, the velocity of light c and Planck's constant H are set 
to c — h — 1. 

In the Standard Model of particle physics, nucleons are constituted of quarks 
and gluons. Each flavor of quark comes in three colors: red, blue, and green. The 
quark color wavefunction can be written as a vector of Dirac spinors: 



(2.1) 



The strong interactions of quarks and gluons are described by Quantum Chro- 
modynamics (QCD), a local 577(3) gauge theory. The dynamics are governed by 
the Lagrangian of QCD: 

£qcd = ^(*7 M D M - M)^ - G^G^/4 . (2.2) 

V v ' * v ' 

quarks+interaction gluons+interaction 

Considering the relevant light quark flavors, the spinor ip q is represented by u, d, 
and s quark (i.e. ip = (u, d, s)). The matrix M in flavor space is composed of the 
bare quark masses, i.e. mJJ, and m", on the diagonal. 
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The gluonic part of the Lagrangian is determined by the gluonic field strength 
tensor: 

G% = d,Al - d v Al + igf^ K\K , (2.3) 

where g and f abc denote the strong coupling constant and the structure constants 
of the group SU(3) P, respectively. 

To achieve invariance of QCD under local SU(3) gauge transformation the 
derivative takes the form: 

D^^-i^-A^. (2.4) 

It induces a coupling between the spin- 1/2 colored quarks fields and the gluonic 
spin-1 gauge fields. The coupling strength g - commonly expressed in terms 
of the strong "fine-structure" constant a s = g 2 /Air - increases with space-time 
distance or equivalently decreases with the momentum transfer Q of a given 
strong process [Zj: 

asiQ) = l + 9a s (A)/47rln(|g2|/ A 2) • ( 2 - 5 ) 

This particular behaviour is a consequence of the self-interaction of gluonic fields 
leading to an antiscreening of the strong interaction which dominates the screen- 
ing of the quark color by quark- ant iquark bubbles. The reference a s (A) is fixed by 
measurements at a certain scale A where a s (A) is small enough to justify a pertur- 
bation expansion, e.g. a s {mz) = 0.118 at the Z boson mass m z = 91 GeV/c 2 

At large distance scales, quarks and gluons are confined in colorless mesons 
and baryons. Only at short distances can perturbation theory be applied, as 
quarks and gluons are quasi-free {asymptotic freedom). The breakdown of per- 
turbation theory for a s > 1 at momentum transfer of about Q ~ lGeV/c, en- 
countered in heavy ion collisions, poses the most challenging problem in theory. 
Both the effective couplings and the relevant degrees of freedom change rapidly 
with scale. 

In the limit of vanishing quark mass, the QCD Lagrangian exhibits additional 
symmetries that can be explored. It becomes invariant under global vector X b 
and axial-vector A 6 7 5 transformations in flavor space: 

^-*exp(-(i/2)\ b a b v )iP q and % - exp(-(i/2)AH7 5 )^ q , (2-6) 

with parameters ay and a\ being arbitrary vectors in flavor space. This results 
in conserved vector and axial-vector Noether currents: 

# >b = ^ 7 M A 6 /2V and j£ b = ^ lB \ h /2^ . (2.7) 

Introducing the quark-spinor projections of right- and left-handed components: 

^R/L = ^(l±7 5 )^ q , (2-8) 
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the transformation 12.61 can be rewritten as: 

- exp(-(i/2)AH)^R , V-l-V-l, (2.9) 
^ L -> exp(-(i/2)A b c^ L , Vr-V>r, (2.10) 

which describes a global SU(3)rxSU(3)l chiral symmetry in flavor space. This 
symmetry has two important implications. First, left- and right-handed quarks 
are not mixed dynamically and their handedness, i.e. the sign of the projection 
of spin on its momentum direction, is conserved. Second, corresponding vec- 
tor and axial-vector resonances are degenerate, as the respective current-current 
correlation functions IIy/ A : 

K/a(q) = * I rf 4 xe^(0|T^ /A (x)^ /A (0)|0) , (2.11) 

which determine the spectral shape of unstable resonances, are identical [H]. 

In the physical world, chiral symmetry is apparently spontaneously broken 
because chiral partners such as p(770)- and ai(1260)-meson show a large mass 
splitting (Am = 500MeV/c 2 ). Therefore, the ground state, i.e. the QCD vacuum, 
is not invariant under chiral transformation. In particular, the vacuum state |0) 
only respects vector symmetries (jy = jR+jh) [HI], while the axial- vector symme- 
try (J a — 3r-3l) is spontaneously broken: 

Q A>b \0) = \PS a ) ^ with Q A:b = J d 3 x ^75^ , (2-12) 

where Q A ,b is the axial-vector charge corresponding to the axial-vector current 
j Ab (see Eq.inD. 

The strength of the symmetry breaking can be characterized by the vacuum 
expectation value of the Goldstone boson [TT] . which is the (nearly) massless pion 
7r: 

(0\f A ,kWM = -iSjkf^e-^ (2.13) 

where f n is the measured pion decay constant of f n = 93 MeV. It is expected that 
a transition from the asymmetric phase observed (see Eq. 12.12)1 to a phase where 
the symmetry of the vacuum is restored (Q Al b\0) — Qv,b\^) = 0) can ^ e triggered 
by external parameters such as temperature and/or pressure. The expectation 
value of the so-called chiral condensate {ip^)'- 

<W = ^(0\uu + dd\0) = (OI^r + ^lIO) , (2.14) 

is the lowest-dimensional order parameter characterizing the chiral phase transi- 
tion. The quark condensate respects all unbroken symmetries of the Lagrangian, 
as it is a scalar density, diagonal in flavor space, and carrying a baryon number 
of zero. 
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The chiral condensate vanishes in the chiral symmetric phase but becomes 
finite in the asymmetric phase corresponding to a mixing of left- and right-handed 
quarks in the ground state. In other words, there is a finite expectation value 
to create a light quark-antiquark pair from a zero-point energy fluctuation of the 
physical vacuum. 

The mixing strength of left- and right-handed quarks in vacuum (ipip) is con- 
nected to the pion decay constant according to the Gell-Mann-Oakes-Renner 
relation [T2*] : 

ml fl = -2m(^V) (m « 6 MeV) . (2.15) 

A value of f n = 93 MeV from pion decay measurements leads to a vacuum ex- 
pectation value of (ipip) — —(240 MeV) 3 = — 1.8 fm -3 , which is large compared to 
the normal nuclear density of about 0.17fm -3 [Hj, indicating a strong dynamical 
breaking of chiral symmetry. 

Until now what has been discussed is the properties of hadronic matter in 
vacuum but what is of more interest is the dynamics in the presence of a hot and 
dense medium. 

When hadronic matter is heated and/or compressed, initially confined quarks 
and gluons start to percolate between hadrons to finally be liberated. This phase 
transition to a plasma of quarks and gluons is accompanied by a melting of the 
quark condensate indicating chiral symmetry restoration. Even before the critical 
region is approached, the chiral symmetry is partially restored by the presence of 
hadrons. The valence quarks and the pionic cloud of a hadron produce a positive 
scalar density inside the hadron, thus effectively decreasing the (negative) quark 
condensate. 

The expected modification of the condensate are derived for the case of high 
temperature T and low density. The equilibrium properties of a hadron gas in 
contact with a heat bath are described by the grand canonical partition function: 

Z(V,T^ q ) = Tr(e-^^ T ) , (2.16) 

where H is the Hamiltonian of the system, N is the quark number operator, and 
fi q denotes the quark chemical potential. The expectation value of the quark 
condensate is then given by the thermal average: 

= Z~ x J2( n \^\ n )e~ {En ~^ q)/T , (2.17) 

n 

where the sum is carried out over all eigenstates of the QCD Hamiltonian. Equa- 
tion 12.171 can be solved for the simplified case of a non-interacting hadron gas. 
The resulting correction reduces the quark condensate for increasing temperature: 

«M>~i_V^E> (218) 
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A <w> p ,T 



Figure 2.1: 



T [MeV] 




P 



Expectation value of the 
quark condensate as de- 
scribed by the Nambu-Jona- 
Lasinio model The 
regions accessible by various 
accelerators are highlighted. 



Each hadron species present with scalar density p s h contributes to the reduction 
of the condensate according to its sigma commutator J2h- The latter quantity 
divided by quark mass is a measure for the integrated scalar quark density inside 
a hadron h: 



The particular temperature dependence of Eq. l2.18l can be explained as follows. In 
the physical vacuum, the color fields are squeezed into hadrons by the repulsion 
of the quark condensate which fills the vacuum. With increasing temperature 
this mechanism becomes inefficient, as thermal pions are produced as excitations 
of the quark condensate, locally changing the expectation value of the quark 
condensate. 

Increasing density also reduces the quark condensate because more and more 
space is occupied by baryons (equivalent to turning on a baryon chemical poten- 
tial). In this case, the nucleons give the dominant correction leading to a formula 
similar to Eq. 12.181 



where Xijv ^ the nucleon sigma commutator and p s N denotes the nucleon scalar 
density at a given nucleon chemical potential fijy. In nuclear matter at normal 
density, the quark condensate is already quenched by 30% according to a value 
of about 45 MeV for the nucleon sigma commutator j!3j . 

The sigma commutator and the dropping of the chiral condensate can be 
estimated in the framework of effective theories, e.g. the cx-model ^3], the 
Nambu-Jona-Lasinio model ^2], or the Walecka model Figure 12.11 illus- 

trates the melting of the condensate for the example of the Nambu-Jona-Lasinio 
model \Tj\. Several effects preceding the phase transition towards a restoration 
of chiral symmetry were predicted by effective mean-field models. The most 
important effects are dropping hadron masses (BR-scaling) [19J and mixing of 
vector and axial- vector currents [20] , both leading to modifications of the hadron 
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Figure 2.2: Lattice calculation including two quark flavors of the temperature 
dependence T ~ 6/g 2 of the quark condensate (ipip) and the associated suscep- 
tibility Xm ~ 5 (ipip) / 5m q |221(left panel) and the energy density e/T 4 and the 
pressure 3p/T A j2H] (right panel). In the limit of an ideal gas of quarks and 
gluons, the energy density should approach a value of e/T 4 = 407r 2 /30 = 17.5 
according to the number degrees of freedom in the plasma phase of 16 and 24 for 
gluons and quarks, respectively. 



spectral functions [21]- The competing models shall be discussed in detail in 
connection with the experimental results presented in Sec. 15. 9. HI 

Model-independent results are obtained by perturbative low-density expan- 
sion. However, this procedure is restricted to temperatures below 120 MeV and 
cannot address the nature of the phase transition. Of the non-perturbative ap- 
proaches, only numerical lattice QCD calculations provide a stringent framework 
even though many-body theory [22] and renormalization-group techniques |23]l24*] 
are promising developments. 

The results of a lattice calculation including two quark flavors are depicted in 
Fig. 12.21 The expectation value of the quark condensate plotted in the left panel, 
as expected, shows the sudden drop at the critical temperature indicating a phase 
transition restoring chiral symmetry. This transition is accompanied by a jump 
in energy density (see upper symbols in the right panel of Fig. I2.2j) from a low 
hadronic value to nearly the level expected for an ideal gas of quarks and gluons 
(QGP) 27 j. Up to now, lattice calculation do not allow to conclude whether 
deconfinement is a phase transition of first order, second order, or just a rapid 
crossover. 

Lattice theory including non-vanishing baryon density has been impeded by 
technical difficulties. Other methods applicable at finite densities include ran- 
dom matrix theory |2"%] EE], random phase approximation |HH] ITTT] . instanton 
models [321 EHl IM| , percolation (23 EH], and supersymmetric models [37]. But 
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Figure 2.3: QCD phase diagram for the transition of hadronic to quark matter. 
An overview of the status of the experimental mapping of the QCD phase diagram 
can be found in |3*H] l3*§] . The theoretical aspects are summarized in 



their predictive power has been rather limited. A general overview and evaluation 
of the available theoretical models is presented in [2~T] . 

The emerging picture of the transition from hadronic to quark matter is il- 
lustrated in a schematic phase diagram in Fig. 12.31 

Despite the tremendous progress in recent years, dispute remains about the 
exact circumstances, at high temperatures and/or high densities, of restoration of 
broken symmetries in the medium under extreme conditions. Whether deconfine- 
ment and chiral symmetry restoration are two distinct phase transitions, or only 
one, is a matter of current debate. Lattice calculation indicate the critical tem- 
peratures of chiral restoration and deconfinement to coincide in the low-density 
scenario. In fact, 01] argues deconfinement in the light-quark sector to be trig- 
gered by the chiral transition. 

Many probes (121 ESI HH US] have been proposed to map out the behavior of 
hot and dense hadronic matter and also to highlight its eventual transition to a 
quark-gluon plasma. Among those probes dileptons stand out for two reasons. 

First, they couple directly to vector mesons. Therefore, hadronic processes are 
expected to reveal their properties in dilepton spectra. In particular, the dilepton 
rate allows for direct measurement of the imaginary part of the current-current 
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correlation function in the medium (see Eq. 12.11)) |4*o1 147]: 

d 6 N + - 2p 2 1 1 

which determines in part the vector meson resonance jlHlHHl- Apart from the 
kinematical constants p describing the meson decay, the influence of the medium 
enters in the current- current correlation function H^ V (T, /j,b) which can be calcu- 
lated by theoretical models discussed in Sec. 15.9.31 

Second, dileptons suffer minimal final-state interaction because they interact 
only electromagnetically and are therefore likely to bring information about the 
innermost zones of high-density and high-temperature matter, formed in the early 
stages of nuclear collisions, to the detector essentially unscathed. On the other 
hand, hadrons are rescattered and carry little information about the time prior 
to the freeze-out stage of the collision. 

The schematic dilepton mass spectrum in Fig. 12.41 indicates the major dilepton 
sources in ultrarelativistic heavy ion collisions. In the high-mass region, dilep- 
tons stem from hard processes (mostly Drell-Yan annihilation) occurring in the 
early pre-equilibrium stage of the collision. Furthermore, a suppression of the 
J/ijj and ip' abundance has the potential to signal the onset of deconfinement, as 
the heavy-quark bound states are dissolved due to freely moving color charges 
(Debye screening). If not masked by an enhanced open-charm production, a 
thermal signal from plasma could be revealed by the observation of associated 
DD production in the intermediate-mass region PJEZj- The low-mass region is 
the exclusive domain of the CERES experiment. It is dominated by soft pro- 
cesses involving the light quark sector. The dilepton continuum originates from 
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Dalitz decays of neutral mesons such as 7r°, 77, 77' — > e + e~7 and a; — > e + e~7i°. 
The resonance peaks are due to direct decays p, u, — > e + e~ . This region is 
particularly sensitive to in-medium modifications of the light hadrons which can 
signal the restoration of chiral symmetry. The p-meson is of special importance 
because, once produced in a dense and hot hadronic environment, it will decay 
predominantly within the fireball due to its short lifetime. Compared to the 
other sources, the relative contribution of thermal dielectron radiation from a 
quark-gluon plasma is expected to be negligible at SPS energies |50j . 

2.2 The CERES physics program 

The production of dileptons in hadronic collisions has been of great experimental 
and theoretical interest for more than 30 years. The early dilepton measurements 
were motivated by the search for the vector mesons in pp, pA, and tt~A collisions. 
The unexpected observation of a continuous dilepton spectrum for mass below 
600MeV/c 2 motivated the development of several theoretical models which are 
still relevant. Most notably was the prediction of the quark-gluon plasma - a new 
phase of matter - in 1978 lj. The thermal radiation of the plasma comprises 
low-mass dileptons and direct photons. 

CERN was the first laboratory worldwide to systematically investigate dielec- 
tron production in ultrarelativistic hadron-nucleus and nucleus-nucleus collisions. 
The Helios-l/NA34 collaboration was the first to measure e + e~ and p + p~ pair 
production in p-Be collisions |5T) . The CERES/TAPS collaboration reproduced 
them with much greater precision. Figure 12.51 shows the measured dielectron in- 
variant mass spectrum of p-Be and p-Au collisions at 450 GeV/c in comparison 
to the expected contributions of hadron decays. The simulation of the so-called 
hadronic cocktail is based on the knowledge of the branching ratios of all rele- 
vant leptonic and semi-leptonic decays and the total production cross sections of 
neutral mesons measured in pp collisions (see App. lAj). 

As evident from Fig. 12 .5[ the hadronic cocktail accounts for the measured 
dielectron yield. Previous speculations about an anomalous source of dileptons 
with mass below 600 MeV/c 2 were found to have originated from an underesti- 
mation of 7] Dalitz yield. This was proven by the exclusive measurement of the 
77 — > e + e~7 decay jHB]- 

Most important, a reference based on pp and pA data was established to 
be used in nucleus-nucleus collisions for distinguishing between new in-medium 
effects and trivial dielectron sources. 

Recent measurements of the low- mass dilepton yield in 200GeV/c p-U colli- 
sions by the NA38/50 collaboration [HI] and in 12GeV/c p-C(Cu) collisions at 
KEK [HE] could not be explained exclusively by the decay of the known hadronic 
sources. The NA38/NA50 collaboration found a significant excess in the mass 
window 0.4-0.6 GeV/c 2 . In-medium modification of the p-meson could not ex- 
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Figure 2.5: Dielectron invariant mass spectrum of p-Be and p-Au collisions at 
450GeV/c 52 . The yield observed is saturated by the expected cocktail of 
hadronic sources. 



plain the observed enhancement. NA38/50 conjectured that it may be due to qq 
annihilations (Drell-Yan process). This production mechanism is considered to 
be negligible in the CERES acceptance, i.e. for low transverse pair momentum 
g t <lGeV/c, but may become important for lepton pairs with large transverse 
momentum predominating in the NA50 acceptance. Therefore, this measurement 
is not necessarily contradicting the CERES results. At KEK, the mass spectra 
of p-C and p-Cu collisions were found to differ significantly below the o;-meson 
peak (i.e. mass window 0.4-0.6 GeV/c 2 ). This difference was interpreted as an 
in-medium modification of the p-meson spectral shape at normal nuclear density. 
Since no such effect was found by CERES, further studies are necessary to settle 
the dispute about the role of in-medium modifications of vector mesons in pA 
collisions. 

The situation changes dramatically for nucleus-nucleus collision. The dielec- 
tron yield observed in S-Au and Pb-Au collisions at 200 A GeV/c and 158 A GeV/ c, 
respectively, significantly exceeds the expectations extrapolated from p-p colli- 
sions [23 GUI EE]- While the n° Dalitz peak is well reproduced by the hadronic 
cocktail, the local minimum expected between the 7/-Dalitz component and the 
p/cu-resonance peak at around 500MeV/c 2 is entirely filled up, as apparent from 
Fig. 12.61 The integrated yield of pairs with mass above 200MeV/c 2 exceeds the 
hadronic cocktail by a factor of 5.0 ± 0.7(stat.) and 2.7 ± 0.4(stat.) for S-Au and 
Pb-Au collisions, respectively. 

This result was confirmed by the observation of an enhanced p + p~ production 
in 200 AGeV/c S-W collisions compared to 200GeV/c p-W collisions by the 
HELIOS/3 collaboration [53] ■ 

The comparison of nucleus-nucleus collision at different bombarding energies, 
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Figure 2.6: Dielectron invariant mass spectrum of S-Au and Pb-Au collisions 56, 
loTj . The Pb-Au data plotted is a weighted average of the 1995 and the 1996 data 
sample [5H] • The data is compared to the expected cocktail of hadronic sources. 
The Pb-Au cocktail plotted was corrected and extended compared to previous 
CERES publications (see App.EJ). The integrated yield for invariant mass above 
200MeV/c 2 of (1. 13±0.16) • 10~ 5 (S-Au) and (5.4±0.9)-l(T 5 (Pb-Au) corresponds 
to an enhancement factor of 5.0 ± 0.7(stat.) and 2.7 ± 0.4(stat.), respectively. 



i.e. different initial conditions, should allow for independent interpretation of 
temperature and baryon-density driven changes of the dielectron spectrum. 

The most recent result for dileptons is the measurement of the invariant mass 
spectrum of Pb-Au collisions at 40 GeV/c per nucleon (see Fig. 12. 7j) . A recently 
resolved problem in GENESIS (see App. |XJ) resulted in a 30% increase of the 
predicted low-mass yield compared to [EH! EI] • The data plotted in Fig. 12.71 were 
taken from 151 and normalized to the expected yield of pairs with mass below 
200MeV/c 2 according to the procedure described in [50] . 

An enhancement of the dielectron yield, larger even than in the 158 A GeV/c 
data, is observed, relative to the expected yield of hadronic sources. Since the 
detector system upgraded in 1998 was not yet fully operational, the data set is 
limited in terms of statistics and momentum resolution. 

The experimental results on dilepton production in nucleus-nucleus collisions 
have experienced many responses from theoretical physicists. These were mainly 
stimulated by the prospects of chiral symmetry restoration and deconfmement 
(see fig. 12. 8j) . Aside from the focus on dielectrons, the CERES collaboration 
has also extensively studied charge particle production [S3], high momentum 



12 



-1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I- 

Pb-Au 40 A GeV 

p ± > 200 MeV/c 

°W >35mrad i 
2. 1 < r| < 2.65 

(NA, =115 




0.2 0.4 0.6 



1 1.2 
[GeV/c 2 



invariant pair mass 



Figure 2.7: 

Dielectron invariant mass 
spectrum of Pb-Au colli- 
sions at 40GeV/c per nu- 
cleon. The data in [HI] 
were normalized to the ex- 
pected low-mass yield sim- 
ulated with GENESIS (see 
App. lA"]) . The integrated 
yield of pairs with invariant 
mass above 200MeV/c 2 of 
(4.2 ± 1.1) ■ 10~ 6 corresponds 
to an enhancement factor of 
4.5 ± 1.2(stat.). 



pions [nHlIM], direct photon production fo5\ I56J. and azimuthal correlations of 
charged particles (i.e. flow) \§T\. The upgrade of the spectrometer with the TPC 
allowed to greatly extend the scope of CERES towards hadronic observables. 
Recent results (HOI IBB] include the measurement of Bose-Einstein correlations, 
mean transverse momentum fluctuations, and hadronic production of A, A, and 
K . 
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Chapter 3 

The CERES detector 



3.1 Setup 

The CERES experiment was designed for the detection of dielectrons with in- 
variant mass up to 2.0GeV/c 2 , covering a range in pseudo-rapidity of 2.1 < rj < 
2.65 [69, 70 . The initial setup consisted of two ring imaging Cherenkov detectors 
(RICH), placed before and after an azimuthally deflecting magnetic field. The 
RICH detectors provide particle identification and a measurement of the trajec- 
tory. The azimuthal deflection in the magnetic field determines the momentum. 
The measurements of p-Be, p-Au, and S-Au collisions (HSIHHIIEE] were performed 
with this configuration. 

In 1994 and 1995, the spectrometer was substantially upgraded in order to 
cope with the high multiplicity environment encountered in ultrarelativistic Pb- 
Au collisions. The original setup was extended by two silicon drift detectors 
(SDD) and a multiwire proportional counter with pad readout (PD) as illus- 
trated in Fig. 13.11 The SDDs sample each track on two additional points. This 
allows for a precision reconstruction of the event vertex, a measurement of the 
energy deposition dE/dx, and a reliable determination of the charged-particle 
multiplicity. The main purpose of the PD is to help the ring recognition in the 
RICH detectors. The enhanced track reconstruction and electron recognition 
capabilities were demonstrated with the study of high-momentum pion and di- 
electron production in Pb-Au collisions at 158 GeV/c per nucleon recorded in 
1995 [03 EH- 

The 158 AGeV/c Pb-Au collision data analyzed in this work was recorded 
in 1996 with the setup shown in Fig. 13.11 The following section gives a brief 
description of the individual components of the experimental setup. 



14 




Figure 3.1: Schematic view of the CERES setup used for the measurement of 
Pb-Au collisions in 1995 and 1996. 



3.2 Target region 

The setup of the target region is shown in Fig. 13.21 The target consists of 8 
gold foils of 600 fim diameter and 25 /im thickness. For the particles to hit just 
one of the consecutive targets disks, a space of 3 mm between targets is chosen, 
minimizing the probability of secondary interactions. The target is surrounded 
by a tungsten shield to protect the readout of the RICH UV-detectors from highly 
ionizing particles scattered backwards. 
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Figure 3.2: 

Schematic view of 
CERES target region. 
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Figure 3.3: Schematic view of a (radial) slice of the SDD (upper panel). The 
interlaced anode structure improves the single-hit resolution by charge sharing 
(lower panel). 



3.3 Silicon drift detector 



The first detector system of the CERES apparatus, located ~ 12 cm downstream 
of the target, is a doublet of silicon drift detectors (SDD) of 4" diameter [72]. The 
CERES experiment was first to successfully employ this type of radially symmet- 
ric position-sensitive detector ;73,. A detailed description of these detectors is 
given elsewhere [Hj. The basic principle of operation is illustrated in the upper 
panel of Fig. 13.31 A charged particle traversing the detector produces a cloud of 
electron-hole pairs inside the depleted region of the semiconductor. The particle 
energy required to create an electron- hole pair is about 3.6 eV [73]. The charge 
deposited by a minimum ionizing particle in a 280 /im thick Si-detector is about 
3.3 fC (20500 electrons) corresponding to an average energy loss of 74keV. 

In radial electric field generated by a set of concentric implanted voltage 



16 



dividers, the electrons drift radially towards a segmented anode at the outer 
circumference of the silicon wafer. The segmented anode, shown in the lower 
panel of Fig. I3.3| consists of 360 partially interlaced pixels to provide a precise 
azimuthal position measurement. The charge signal collected for each anode is 
digitized by a fast FADC with a sampling frequency of 50 MHz. Given the drift 
velocity known, the radial position of a charged-particle hit can be determined 
by a measurement of the drift time with respect to the first-level trigger. 



3.4 Ring Imaging Cherenkov detector 

The essential components of the 1996 CERES apparatus are two Ring Imaging 
Cherenkov detectors (RICH). The first of these is situated between the SDDs and 
a short super-conducting double solenoid, and the second is behind the solenoid. 
An electron produced in a collision emits Cherenkov photons while traversing 
the Methane filled radiator volume. A spherical mirror reflects the Cherenkov 
light to form a ring image at the mirrors focal plane. In case of RICH-2, the 
mirror is split up in 8 smaller panels of equal size for manufacturing reasons. 
The geometry of one such panel is shown in Fig. 13.41 The imaging properties are 




described in detail in Sec. 14.2.21 The photons are detected in a position sensitive 
UV-detector located at the image plane. This gas detector is filled with helium, 
methane, and Tetrakis-dimethylamino-ethylene vapor (TMAE) [ZH]. The high 
quantum efficiency of TMAE in the UV region made it the preferable detector 
gas at the time when CERES was designed. The UV-detectors of RICH-1 and 
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RICH-2 are separated from the radiator volume by a CaF 2 and a quartz window, 
respectively. 

Electrons emitted after photoabsorption in TMAE are amplified in three steps: 
two parallel-plate avalanche stages and a subsequent Multi-Wire Proportional 
Counter (MWPC). After an amplification by a factor of 2-5 • 10 5 , signals are read 
out on 50000 individual pads covering the geometric acceptance. A complete 
description is given in [7UJ [771 EE] • 



3.5 Magnetic field 

To determine particle momentum, a magnetic field of 7 T is produced between 
RICH-1 and RICH-2 detector by a pair of super- conducting coils with alternating 
currents. The currents in the additional correction coils are adjusted for the field 
lines in the RICH-2 radiator to point back to the target and the RICH-1 detector 
to become almost magnetic field free. This particular field shape is illustrated in 
Fig. ETT1 

Particles crossing the radially increasing magnetic field {B ~ 1/r) between 
RICH-1 and RICH-2 are azimuthally deflected by an angle A0: 

144 mrad „ . 
A0 « GeV/c . 3.1 

V 

The polar angle of the particles is approximately conserved because particle tra- 
jectories and magnetic field lines in RICH-2 do not cross in this direction. How- 
ever, all particles traversing the RICH-2 detector are slightly bend towards the 
beam axis according to their initial azimuthal deflection. This so-called second- 
order-field effect is illustrated in Fig. 13.51 

While RICH-2 measures the local space direction of the trajectory in the 
radiator, the Pad Chamber samples a point on the particle trajectory with respect 
to the vertex. Therefore, the azimuthal deflection observed in RICH-2 is about 
1.5 times larger than in the pad chamber. 

To simulate the particle trajectories in the presence of the B-field, the Poisson 
program package [7^] was used to calculate a field map for the CERES geometry 
because the magnetic field distribution has not been measured directly. The 
accuracy of the map was estimated to be better than AB / B ^0.5%. 



3.6 Pad Chamber 

The Pad Chamber (PD) is located downstream of the RICH-2 radiator. It con- 
sists of a MWPC with pad readout. Each of the 20000 pads has a dimension of 
7x7 mm 2 , which results in an angular resolution of better than 0.6 mrad. The 
main purpose of the PD is to limit the number of possible electron track candi- 
dates, given by the combination of rings in the RICH-1 and RICH-2 detectors, 
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by measuring an additional point on the particle trajectory. This background 
reduction is vital in the high charged-particle multiplicity environment of Pb-Au 
collisions. Furthermore, the PD provides an independent reference frame for the 
geometry calibration of the other detectors jHS|- An exhaustive description is 
given in [80 . 



3.7 Trigger system 

The CERES trigger system consists of three Cherenkov counters (BC1, BC2, 
BC3) and two plastic scintillator detectors: veto counter (VC) and multiplicity 
detector (MD). BC1 is located 60m upstream of the experiment, BC2 and BC3 
are directly before and after the target, respectively. MD is at the downstream 
end of the spectrometer. 

The minimum-bias collision trigger requires a lead signal in BC1 and BC2, 
and no signal in BC3. The central trigger additionally demands a signal in MD, 
with the threshold set at a level corresponding to 100 charged particles. This 
is equivalent to the most central 35% of the geometrical cross section. To avoid 
beam pile-up, the particle which triggers the reaction must not be followed or 
preceded by any other lead particle going through BC1 for several microseconds. 



3.8 1998 detector upgrade 

The CERES detector system was upgraded in 1998 by the addition of a new 
magnet system and a radial Time Projection Chamber (TPC) [HI]. The PD and 
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Figure 3.6: Schematic view of the upgraded CERES setup used for the measure- 
ment of Pb-Au collisions in 1998,1999, and 2000. 

the multiplicity detector were removed. This new setup is illustrated in Fig. 13.61 

There are three major benefits of the upgrade. First, the high-precision mo- 
mentum measurement of the TPC is expected to improve the mass resolution 
to dp/p = ^(0.0105)2 + (0.0103 • p ■ GeV _1 c) 2 jHE]- Second, between the RICH 
detectors no magnetic field is applied, allowing to operate both detectors in a 
combined way with doubled photon statistics for the RICH rings and increasing 
efficiency |60j . Third, considerable reduction of the combinatorial background 
can be achieved due to higher photon statistics in the RICH detector and the 
additional dE/dx measurement in the TPC. Furthermore, the spectrometer ca- 
pabilities for the study of hadronic observables are significantly extended. 
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Chapter 4 



Development of an 
ultralightweight mirror for RICH 
detectors 

4.1 Reasons for the replacement of the RICH-2 
mirror 

After the upgrade of the CERES experiment with a TPC downstream of the ex- 
isting detector as described in Sec. 13.81 a new tracking scheme has been developed 
in order to improve background rejection. According to this scheme, both RICH 
detectors are operated without magnetic field and allow a combined use for elec- 
tron identification and tracking while the momentum is measured separately in 
the TPC. All particles must traverse the RICH-2 mirror before entering the TPC. 
Multiple scattering in the mirror material changes the particle direction and, thus, 
deteriorates the invariant mass resolution. Additionally, electrons lose energy by 
Bremsstrahlung. The resulting low energy tail impedes the spectroscopy of vec- 
tor resonances. Therefore, the replacement of the thick RICH-2 glass mirror by 
an ultralightweight mirror almost transparent to dielectrons would considerably 
improve the performance of the new detector system as will be discussed in detail 
in the following sections. 

4.2 Impact of the RICH-2 mirror on the spec- 
trometer performance 

4.2.1 Interaction of electrons in matter 

High-energy electrons traversing the matter of a mirror are affected in two ways. 
First, all electrons with momentum p are deflected due to multiple Coulomb 
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E (MeV) 



Figure 4.1: Fractional energy loss per radiation length in lead as a function of 
electron/positron energy jS]. 

scattering from nuclei with a probability depending on the thickness x and the 
radiation length X of the material. The Coulomb scattering distribution is 
roughly Gaussian for small deflection angles, with a width given by: 



Second, electrons lose energy by bremsstrahlung at a rate nearly proportional to 
their energy. The cross section can be approximated in the "complete screening 
case" as jH]: 



with y=k/E being the fraction of the electron's energy transferred to the radiated 
photon. The energy loss due to scattering and ionization is negligible compared 
to bremsstrahlung for electrons with E > 50MeV/c 2 as illustrated in Fig. 14. II 8j. 
The thickness of matter in terms of electron energy loss can be conveniently 
measured in units of the radiation length X which is the mean distance over 
which a high-energy electron loses all but 1/e of its initial energy by brems- 
strahlung. The radiation length can be approximated for a material with a charge 
number Z and a mass number A as [8 u : 




(4.1) 




(4.2) 



716.4 gcmr 2 A 



(4.3) 



o — 



Z(Z + 1) ln(287/\/Z) ' 
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For compounds or mixtures each material contributes with a fraction Wj propor- 
tional to its molecular weight Af 



^ = Ef = ET^- (4-4) 

j Aj j ^comp -A-j 

The material thickness, which corresponds to 0.01X , is given in Table I4~T1 for 
several materials. 



material 


density 
in g/cm 2 


thickness 
(x/Xq = 0.01) 


comment 


soda lime float glass 


2.5 


1.4 mm 


CERES RICH-2 mirror 


carbon fiber (MAN) 


1.7 


2.6 mm 


CERES RICH-1 mirror 


carbon fiber (HEXCEL) 


1.6 


2.7 mm 


COI prototype 


Poly-Carbon ceramic 


1.54 


2.8 mm 


HADES mirror [S2| 



Table 4.1: Equivalent thickness of mirror materials 



4.2.2 Imaging properties of the RICH detector 

A brief review on the fundamentals of RICH detectors with special emphasis on 
the CERES RICH detector will precede more extensive discussion of the impact 
of the RICH-2 mirror on the performance of the CERES spectrometer. 

A Ring Imaging Cherenkov detector measures the photons radiated by a 
charged particle traversing a transparent medium (radiator) with a velocity higher 
than that of light in the medium. The emission angle 9 of the so-called Cherenkov 
photons is then determined by the index of refraction n of the radiator medium 
and the velocity (3 of the charged particle: 

cos 6> photon = . (4.5) 
np 

At atmospheric pressure gases have refraction indices close to one. The threshold 
velocity is best expressed in terms of the Lorentz factor 7 th : 

7th =(l-lp. (4.6) 

Methane (CH 4 ) with 7 t j, = 32 was chosen for the radiator gas of the CERES 
RICH [77] . It makes the detector almost blind to hadrons, except to pions with a 
momentum of more than 4.5GeV/c. The Cherenkov photons emitted along the 
trajectory of the particle in the radiator and reflected by a spherical mirror create 
a ring image at the focal plane. The particular optical scheme is illustrated in 
Fig. 14.21 A photon detector, located at the focal plane, allows to determine the 



23 



UV detector plan 



ring image . , 
♦ ♦ 



UV mirror 



Cherenkov 
photons 



electron track 



R 



Pb beam 



target 



CH4 chamber 



Figure 4.2: Schematic view of the CERES RICH detector illustrating the origin 
of ring images. 



ring position and radius, as well as the number of photons (see Fig. I4.3|) . 

The ring-center position is a measure of 9- and 0-coordinates of the original 
charged-particle track. The ring radius is related to the photon emission angle 
and, thus, to the velocity of the particle. The Cherenkov angle approaches its 
asymptotic value 9 max for electrons with a velocity close to the speed of light 

mx9 mKX = — • (4.7) 
7th 



Formula 14.71 applied to CERES geometry results in an electron ring radius of 
Roo ~ 31.25 mrad, slightly depending on radiator temperature and atmospheric 
pressure. 

The efficiency of the ring reconstruction depends strongly on the number of 
detected photon per charged particle. The total number of photons N p h emitted 
per unit path length x and unit photon energy interval E p ^ is related to the 
half-angle # p h of the Cherenkov cone: 



d 2 K 



P h 



dEph dx he 



— sin 9 ph = — (1 - n C n 4 {E ph ,T) 



I) 



(4.8) 



and, thus, to the index of refraction of the radiator gas ncH 4 which is a function 
of photon energy. 

The energy range of the detected photons is limited by the photo-sensitivity of 
TAME {E > 5.4 eV) jH3] and the ultraviolet transparency of the quartz entrance 
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Figure 4.3: 

Asymptotic electron ring re- 
constructed from 11 photon 
hits in RICH-1 |H3j. 



window (E < 7.4 eV) [HS1 [H^j which separates the UV detector from the radiator. 
In this energy range, the index of refraction is very close to one and nearly 
constant. The temperature gradient in the radiator is small due the continuous 
gas flow. Therefore, the expected total number of photons for a radiator of length 
L is given by: 

^LiV (sin 2 # ph > . (4.9) 

The factor iVo, called the figure of merit, is defined by the product of quantum 
efficiency of the UV-detector Q, total transmission of radiator gas and quartz 
entrance window T, and mirror reflectivity R: 

a r 

N = — Q(E)T(E)R(E)dE. (4.10) 
Tic J 

The emission probability for k Cherenkov photons is described by a Poisson dis- 
tribution: ^ 

P(N = k) = ^-exp(-(N)) (k = 0, 1, 2, ...) , (4.11) 

where (N), the mean number of photons, is given by Eq. 14.91 Successful detection 
of an electron requires the reconstruction of a ring image composed of k single 
photon hits. The reconstruction efficiency depends strongly on the mean number 
of photon hits per ring excluding background contributions jHZj- Thus, it is 
limited by all contributions in equation 14.101 in particular by the UV reflectivity 
of the mirror coating. Additionally, the UV detector spokes and small gaps 
between adjacent mirror segments lead to a local reduction of the number of 
detected photons. 

The right assignment of a RICH ring to the external track information is of 
great importance for efficient particle tracking. It is determined by the ring- 
center resolution crR ing and, thus, by the mean number of hits per ring (N) and 
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the single-hit resolution (T hit according to: 



^Ring = y ^ _ 2 ffhit • (4.12) 

There are four major contributions to the single photon hit resolution. First, 
multiple scattering of the charged particle within the radiator contributes on 
average cr mult « 0.26 mrad. Second, the chromatic dispersion of the radiator gas 
results effectively in a smearing of Cherenkov angle according to: 

^disp = \^zr\ B , (z 3 » i; n « i) [HZj , (4.13) 

cr n being the rms width of the index of refraction averaged over the bandwidth 
and weighted with the probability to detect a photon QTR (see Eq. 14.90 . With 
Odisp ~ 0.53 mrad, chromatic dispersion [HI] is by far the dominating contribution 
to the single-hit resolution. Third, mirror shape irregularities that occur on a 
scale of less then the radius of a Cherenkov light cone distort the ring image and, 
thus, contribute to the single-hit resolution. This is one of the main issues of this 
chapter. Large scale mirror deformations shifting the entire ring can be corrected 
to first order by local adjustment of the focal length, provided the deformations 
are continuous for adjacent mirror segments. Finally, the granularity of the UV- 
detector as determined by a pad size of (2.7mm) 2 for RICH-1 and (7.6 mm) 2 for 
RICH-2 (equivalent to about 2 mrad per pad in both cases) results in an expected 
single-hit resolution of <r pa d = 1.8(1.4) mrad for the RICH-1 (RICH-2) detector. 

To protect the UV-detectors from particles produced in the collision, the 
target is placed at 0.8 • f m i rror which leads to a small deviation from the ideally 
flat focal plane. The contribution thereof to the single-hit resolution is negligible. 

All contributions are independent and, hence, add in quadrature to the single 
photon resolution: 

°hit = \J Vmult + °disp + a mirror + a pad ■ ( 4 -14) 

4.2.3 Simulation of energy loss in the RICH-2 mirror 

The effect of bremsstrahlung was studied using a GEANT detector simulation [SB] 
including contributions of multiple photon radiation in the inhomogeneous ma- 
terial distribution of the CERES setup according to Eq. 14.21 

Figure 031 shows the relative energy loss of electrons due to bremsstrahlung in 
the CERES setup upstream of RICH-2 (« 1.3% of a radiation length), which is 
dominated by the target and the RICH-1 mirror, and the additional contribution 
of an ultrathin RICH-2 mirror (0.5% of a radiation length) in comparison to the 
present thick mirror (4.5% of a radiation length). 

Considering the logarithmic scale in Fig. 14 A\ it becomes obvious that the 
RICH-2 glass mirror is presently the dominating source of bremsstrahlung leading 
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to a significant low energy tail for every resonance in the dielectron pair mass 
spectrum. To verify this for the mass region of interest, the dielectron decay 
distributions of 0- and cu-mesons were folded with the bremsstrahlung spectrum. 
Figure 14.51 shows the result for the case of the u resonance. 

The number of dielectrons in the peak drops significantly due to the brems- 
strahlung tail. Discriminating the resonance peak from the background of other 
sources becomes increasingly difficult. An ultrathin mirror, in comparison, would 
significantly reduce the dielectron loss almost to the minimum level determined 
by the upstream material. Although not shown here, the situation of the 0-meson 
is qualitatively comparable. 

4.2.4 Quantitative estimate of the influence of the RICH- 
2 mirror on the invariant mass spectrum 

A dielectron invariant mass spectrum was produced for quantitative study us- 
ing the Monte Carlo detector simulation. The input distributions of various 
dielectron sources were obtained from the GENESIS e~e + event generator (see 
App. EJ. The appropriate momentum resolution and the CERES acceptance 
cuts were applied to the simulated mass spectrum. 

To perform meson spectroscopy, the dielectron signal S^nke has to be ob- 
tained by subtraction of the uncorrelated background pairs -B un iik e according to 
Eq. 15. 1 II in Sec. 15.61 Next, the resonance signal S meson can be extracted from the 
background of other dielectron sources -B mes0 n: 

^meson ^unlike -Bmeson • (4.15) 
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Figure 4.5: 

Impact of the RICH-2 mir- 
ror on the spectrum of the uj 
resonance. A Breit-Wigner 
function was assumed for 
the spectral shape of the uj 
resonance. 
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For an expected signal-to-background ratio of Unlike /-Bunlike ~ 1/10 in the mass 
region of u- and ^-resonance, the relative error of the signal AS me son/S m eson is 
dominated by the statistical error of the uncorrelated background of dielectron 
pairs B unhkc : 



VSunlike + 2B unlikc + \/B 

meson + 2B unlike \/8B 



unlike 



In the next step, the meson signal is determined for each mirror version inte- 
grating the pair yield within the 3<r-width of the original resonance peak (see 
Fig. 14. 5 J) . The results are summarized in Fig. 14.71 and Table l4~2l It turns out 
that the number of dielectrons in the peak would increase by as much as 30% if 
the present thick mirror was replaced by an ultrathin mirror. The impact of the 
30% difference is most profound in a low-statistics and high-background scenario 
because large width of mass bins and large statistical errors of the background 



Option 


Total number 
in u peak 


% in 
3cr width 


Total number 
in peak 


% in 
3o" width 


1.3% RICH-1 


2300 


89.8 


500 


89.2 


+0.5% RICH-2 (thin) 


2300 


83.0 


500 


81.9 


+4.5% RICH-2 (thick) 


2300 


64.6 


500 


62.9 



Table 4.2: Impact of bremsstrahlung on the number of dielectrons in the u- and 
the ^-resonance peak. 
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subtraction (see Eq. I4.16|) strongly reduce the statistical significance of the me- 
son signal. The reduction of the multiple scattering in the thin mirror results 
in an improvement of the total momentum resolution by approximately 0.5% at 
p = 6 GeV/c and 1% at p — 10 GeV/c as seen in Fig. 14.61 

It should be noted that in case of the thick mirror, shown in panel (4) of 
Fig. 14. 1\ the region between u- and 0-meson peak, which is particularly sensitive 
to possible shape changes of the p-meson, would be masked by the bremsstrahlung 
tail of the 0-meson resonance. 

To summarize, spectroscopy of u- and 0-mesons would greatly benefit from 
replacement of the present RICH-2 mirror by an ultrathin mirror. An increase 
in signal-to-background ratio of the meson resonances and better access to the 
spectral shape of the p-meson peak are the two main prospects. 



4.3 Ultralight weight RICH-2 mirror 

In the past, manufacturing of ultralightweight mirrors has proven a technologi- 
cally very challenging endeavor for two reasons [51] • First, reduction of the mirror 
thickness results in a quadratic decrease in bending stiffness making it increas- 
ingly difficult to maintain sufficient optical imaging quality. Second, achieving 
a high reflectivity for photon energies in the UV range requires not only a sur- 
face micro-roughness below 3 nm but also very sophisticated coating technol- 
ogy esi nu mm • 
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Figure 4.7: Impact of RICH-2 mirror on the dielectron mass spectrum. 

(1) cocktail of expected hadronic sources (HHI 

(2) including the limited momentum resolution 

(3) including the modifications by bremsstrahlung in RICH-1 mirror 

(4) including the modifications by bremsstrahlung in RICH-2 mirror 
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L length scale 



4.3.1 Mirror distortions 

Deviations of the mirror surface from an ideal spherical shape will result in a 
blurring of the point image produced by parallel radiated Cherenkov rays. The 
deformation of the mirror can be measured in terms of the slope error. It can 
be specified by a twofold of the angle between the actual slope of the mirror 
surface and the nominal value which corresponds statistically to a rms-width of 
distribution of the reflected light (see Fig. I4.8p . 

In a Monte Carlo simulation including real background events of the 1995 data 
set, the ring reconstruction efficiency (defined as the probability to reconstruct a 
ring for a given Monte Carlo track) and the ring-center resolution in RICH-2 were 
studied as a function of the slope error of the mirror. Local distortions smaller 
than the mirror area illuminated by a Cherenkov light cone were assumed to be 
randomly distributed. This is roughly equivalent to a random deviation of the 
reflected light from its nominal direction with a Gaussian probability distribu- 
tion. Therefore, the hit position of each Cherenkov photon in the UV-detector 
is smeared by the convolutions of the probability distributions describing slope 
error, chromatic dispersion, electron drift, and finite pad resolution. 

Figure 14.91 shows the decrease of the ring reconstruction efficiency with in- 
creasing slope error. This reflects the decrease of the number of true hits per 
ring due to the fixed size of the ring search mask in the reconstruction algorithm 
and the relative increase of misidentified hits on rings. This analysis allowed 
for setting of the limits of tolerance for slope errors at a maximum of about 
0.7mrad. It is worth stressing that improving the mirror quality beyond a slope 
error of 0.25 mrad will not lead to a gain in reconstruction efficiency. Further- 
more, it should to be noted that the ring-center resolution deteriorates quickly 
with increasing mirror deformation as seen in Fig. 14.91 It is clear that as soon 
as the errors induced by a certain slope error become larger than the contri- 
bution of the chromatic dispersion (0.53 mrad) the ring-center resolution drops 
significantly. Optimizing the tracking strategy and fine tuning of the ring fitting 
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Figure 4.9: Ring reconstruction efficiency (left panel) and ring-center resolution 
(right panel) of RICH-2 as a function of the slope error of the RICH-2 mirror. 



algorithms could slightly improve the result. Combining both studies, an up- 
per limit of 0.7mrad(rms) for the slope error was estimated. The present glass 
mirror, in comparison, has an overall slope error of 0.35 mrad. Additionally, the 
actual focal length of the mirror must not deviate by more than 1.0% from the 
nominal value of / = 4000 mm. The gap between adjacent mirror segments must 
not exceed 2 mm which corresponds to an average loss of about 1.8% of the pho- 
tons for 50% of all rings. A small gap would allow for an easy adjustment of the 
mounted mirror segments. 



4.3.2 Reflectivity for UV photons 

The only metal known to provide high reflectivity in the UV region is aluminum 
(Al). The optical properties of metals are characterized by the index of refraction 
n and the extinction coefficient k. The normal-incidence reflectance is given by: 



(n - l) 2 + k 2 

«■ = bb ■ (4 - 17) 



which holds true for Cherenkov photons because of the large focal length of the 
mirror and the small Cherenkov angle. A maximum reflectance of 95% can be 
achieved for aluminum at a wavelength of 200 nm j^H] . In case of the size of 
the surface roughness exceeding 10% of the photon wavelength A, the diffusely 
reflected component of the beam becomes large due to scattering from surface 
structures according to: 

exp — , (4.18) 



-Rn V A 
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with a being the width of the autocorrelation function (z(r)z(r')) of the sur- 
face profile. The latter describes the spatial correlation between the height of 
the surface at points z(r) and z(r') which can be approximated by a Gaussian 
distribution: 

K(r - r') = (z(r)z(r')) « a exp (-(r - r') /L) , (4.19) 

where a and L denote the micro-roughness and the average correlation length, re- 
spectively. Thus for UV photons, a micro-roughness of less than 3 nm is required, 
comparable to excellent polished float glass. 

The most widely used technique for depositing Al coatings is evaporation 
in high vacuum All coatings for this study were performed at the CERN 
coating facility. Purest grade of Al(99.9999%) was evaporated from a tungsten 
coil. It is then deposited on the rotating substrate to ensure uniform thickness 
of the coating. In order to produce highly reflecting films, extreme care must be 
taken to ensure that the evaporated coatings are not contaminated by residual 
gases present during deposition. In principle, the metal should be deposited 
at the highest possible rate and not thicker than needed to be just opaquely 
reflecting. Otherwise, the surface roughness will increase. The optimal thickness 
for Al was found to be about 7 A 97 . The small thickness of the film means a 
rather slow deposition rate and, therefore, an ultrahigh vacuum of 10~ 7 mbar was 
of utmost importance. The thickness of the coating was controlled by measuring 
the oscillation frequency change of a crystal induced by the material deposited 
onto it. 

A natural oxide film grows on the evaporated Al surface to an ultimate 40 A in 
thickness as soon as it is exposed to air. While this oxide layer prevents aluminum 
from tarnishing, it also causes a drastic decrease in reflectance in the UV region. 
Therefore, the Al film needs to be protected with a magnesium fluoride (MgF2) 
overcoating preventing oxidation. The thickness of the MgF2 film is chosen such 
that the reflectivity of the combined layer is enhanced by destructive interference 
of the light reflected on both boundaries of the film. A MgF 2 layer of 3.1 A in 
thickness was used for all coatings produced for this study. 
It is well known that evaporated MgF2 coatings absorb water when exposed to 
air. This decreases the reflectance in the UV region significantly. Therefore, once 
coated, a mirror must always be kept free from moisture, preferably in a protective 
nitrogen atmosphere. Additionally, outgassing or diffusion of a component of the 
mirror substrate can deteriorate the coating. Previous experience in the coating 
of the RICH-1 mirror showed that a replicated gold surface needs to be covered 
by a blocking layer such as SiO / Cr to prevent the gold atoms from diffusing into 
the aluminum coating [94J. 
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4.3.3 Other important aspects 

The mirror is operated in a methane (CH 4 ) environment at atmospheric pressure 
in a temperature range of 35 to 45 °C. All specifications must be fulfilled at least 
for this operating temperature interval. Furthermore, the mounted mirror should 
withstand any thermal stresses induced in a temperature interval of 10 to 45 °C 
due to shut down of the heating during the off-line period. 

Any additional contamination of the radiator gas with lppm of water or 
oxygen will reduce the UV photon efficiency by about 1% due to absorption. 
The present gas system can achieve an equilibrium concentration of O2 < 1 ppm 
and H 2 < 2 ppm. If a new mirror introduces a larger oxygen source an upgrade 
of the RICH-2 gas system would be required to maintain the present low level 
of oxygen contamination. Of further importance is the fact that cleaning of the 
CH4 radiator gas after opening the RICH-2 takes at least 3 to 6 weeks. 

4.4 Manufacturing technologies 

The manufacturing technology depends on the choice of the mirror substrate 
such as glass, metal, or composite materials. Traditionally, mirrors for optical 
applications such as the old RICH-2 mirror were made of glass for the following 
reasons: technological control over wide range of physical properties, excellent 
surface quality due to polishing, and very low production cost. For the purpose 
of reducing the thickness below 1% of a radiation length, a float glass mirror 
must be thinner than 1.4 mm according to Table 14.11 This, however, is not 
feasible because glass as a quasi-fluid rapidly loses its long term shape stability 
for large-size mirror panels with a thickness of less than 3 mm. 

In the scope of this study, three alternative approaches were considered. First, 
to construct a mirror from coated Mylar foil keeping a spherical shape by applying 
a pressure gradient between front and back side. A two Mylar foil mirror setup 
of a 100 (Urn thickness for each would correspond to 0.1% of a radiation length. 
This option was rejected because any spatial anisotropy in elasticity or substrate 
thickness will lead to large shape distortions. Second, a novel mirror substrate 
based on carbon ceramic was developed by DASA/IAGB j^H] and the Technical 
University of Munich for the HADES RICH mirror jS2j- The extremely high stiff- 
ness of the ceramic allows for minimum mirror thickness of about 2.0 mm. This 
corresponds to 0.75% of a radiation length for a substrate density of 1.54 g/cm 2 . 
Additionally, a very high local surface quality and, hence, reflectivity can be 
achieved because substrate shells are individually polished to a micro-roughness 
of about 2 nm. The major disadvantage of this type of substrate is the high risk of 
residual stresses in the material. At the time of this evaluation the mirror panels 
regularly broke during polishing or trimming indicating insufficient control of the 
manufacturing process. In combination with the prohibitive high cost of such a 
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Figure 4.10: Structure of a uni-directional carbon-fiber material. 



mirror, this led to discarding of the option. 

Third, mirror substrates made of carbon-fiber-composite (CFK) materials 
were considered. These materials had been used for many high-tech applica- 
tions in the past decade. Many times, it was successfully shown that large ul- 
tralightweight mirrors can be manufactured using CFK substrates with a replica 
technology for coating [HU El EH] • The main advantages of CFK materials are 
very low area density (pa ~ 1.7g/cm 2 ), high stiffness, and industrial mass pro- 
duction guaranteeing reasonable prices. The significant technological advances 
of recent years led to a considerable reduction of inherent disadvantages such as 
fiber printing, inhomogeneous module of thermal expansion, and moisture ab- 
sorption. In view of these promising developments and the lack of alternatives, 
all subsequent studies were focused on CFK materials. 

4.4.1 Carbon-fiber-composite mirrors 

Carbon-fiber-composite mirrors consist of bundles of carbon fibers which are 
glued together with a special resin in a baking process. Three types of sub- 
strates can be distinguished depending on the arrangement of the carbon-fiber 
bundles: uni-directional materials with irregular orientation of chopped bundles, 
materials with long parallel aligned bundles, and woven fabricates consisting of 
interlaced texture of carbon-fiber bundles. Figure l4~. lOl illustrates the composition 
of a typical uni-directional material commonly used as mirror substrate. These 
days, the properties of carbon-fiber substrates can be controlled within a broad 
range to fit specific requirements. This can be done by choosing the proper fiber 
material and resin system and by adjusting the parameters in the curing process. 
Most difficulties of CFK substrate arise from the spatial inhomogeneity of the 
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material due to anisotropic fiber structure and the thermal inhomogeneity due 
to the difference of the coefficient of thermal expansion (CTE) between resin and 
carbon fibers. The lower CTE of the resin system creates valleys between adjacent 
carbon-fiber bundles when the substrate cools down after the high-temperature 
curing process as shown in Fig. 14.101 

This effect, known as fiber printing, can be reduced by applying an addi- 
tional resin layer to the substrate surface in a replication process as illustrated in 
Fig. 14.111 This process can also be utilized to transfer a reflective metal coating 
such as aluminum and silver to a mirror surface. 

This technology was successfully applied in the manufacturing of the CERES 
RICH-1 mirror by MAN Technologie AG |H4" t 1101] . However, it cannot be uti- 
lized for UV coatings, as the protective MgF2 coating layer is very fragile and, 
therefore, cannot be transferred from master to substrate. 

Another drawback of CFK materials is their spatial inhomogeneity. This 
can be controlled by combining several uni-directional carbon-fiber sheets with 
appropriate fiber orientations. Approximate spatial homogeneity to thermal, me- 
chanical, or bending stresses can be achieved through arrangement of 2, 4, or 8 
layers of CFK, respectively. Adherence to absorption of moisture causing the ma- 
terial to expand and to act as a source of water in a dry atmosphere such as the 
CH 4 radiator gas is a further disadvantageous property of CFK materials |l()2j . 

4.4.2 Evaluation of CFK manufacturing concepts 

A survey of CKF mirrors from all major vendors was carried out to assess the 
impact of recent improvements in manufacturing technologies on the inherent crit- 
ical properties: surface micro-roughness, fiber print-through, and substrate stiff- 
ness. Mirror samples were obtained from the following companies: ARDCO 103J, 
SESO [El!, and Composite Optics Inc. (COI) jUS]. MAN Technologie AG [TTTT] . 
that had previously built several 1-m-diameter CFK mirrors for CERES RICH- 
1, had no longer the facilities to do the gold coating of the replication master 
and could not offer any alternative solution. Before assessing technologies with 
respect to the specifications, experimental methods to measure the samples will 
be described. 

UV reflectivity can be measured directly with a reflectometer or indirectly by 
determining the micro-roughness with a surface interferometer. The samples were 
evaluated in the CERN reflectometer depicted in Fig. 14.121 A grating selected a 
monochromatic light beam with a wavelength of 160 to 230 nm with a bandwidth 
of 2nm from the continuum of a deuterium lamp. The intensity of the incident 
beam I was measured with the photo multiplier PM1. After rotation of the 
central mirror the beam was sent to the center of the measured sample. The 
ratio of the intensity of reflected beam J re fi measured with photo multiplier PM2 
to the intensity of the incident beam Jo gave the reflectance of the sample. All 
measurements were performed in ultrahigh vacuum to avoid photon absorption. 
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epoxy is cured at an 
elevated temperature 




4. The substrate and the 
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seperated from the master 
by mechanical force and 
the release agent is 
removed with water 




Figure 4.11: Replication process for a sandwich-design CFK mirror developed by 
OPTICON E 
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Figure 4.12: Setup for UV reflectivity measurement at CERN |§T] . 



To evaluate the CFK material specific fiber print-through a surface profile 
needs to be measured with a resolution of less than 0.1mm as determined by 
the diameter of the carbon-fiber bundles. Figure 14.131 shows the experimental 
setup used to scan flat reflective samples with a laser beam. The UV laser excels 
because of a small beam diameter of 300 fim and a very high intensity that allows 
to measure the reflection from uncoated samples. A surface profile was obtained 
by moving the sample on a linear stage relative to a fixed laser beam and recording 
the image of the reflected laser beam on a position sensitive diode. Any surface 
deformation with a length scale larger than the beam diameter results in a shift 
of the image position on the diode. The average slope error A0 can be computed 
from the rms-width of the distribution of image positions: 



A(f) = arctan (tana + ((x — x) 2 )/s 



(4.20) 



where s denotes the distance between mirror and diode and a is explained in 
Fig. 14.131 The orientation of the carbon-fiber bundles needs to be taken into 
account for highly anisotropic uni-directional CFK fabric. 
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Figure 4.13: Setup for the measurement of the local slope error of flat mirror 
samples. 



The new substrate from ARDCO is a sandwich of CKF layers and ROHA- 
CELL foam. The ROHACELL foam stiffens the mirror considerably and increases 
the total thickness only by 0.3% of a radiation length. The mirror samples man- 
ufactured by ARDCO |l()3j represented a novel sandwich design developed for 
the RICH mirrors of the PHENIX experiment [§H]- Two uni-directional plane 
CFK shells of 1 mm in thickness are stabilized and stiffened by a 1 cm thick RO- 
HACELL foam core. A detailed analysis of the sandwich design concept will be 
presented in Sec. 14.4.31 

The composite substrates were coated by OPTICON Inc. 400] with an Al/Au 
layer by a replication process (see Fig. I4.11j) . To achieve high UV reflectance, 
both samples were coated with a high quality Al and MgF2-protective film (by A. 
Braem, CERN). This required the development of a novel technology involving 
differential pumping on the sandwich structure to prevent outgassing and collapse 
of the ROHACELL foam in the ultrahigh vacuum during coating. 

A summary of the measured reflectance before and after UV coating is shown 
in Fig. 14.141 A reflectance of 87% for a wavelength above 180 nm is close to 
the expected value of 90% jHS] ■ Furthermore, it is comparable to the float glass 
witness sample indicating that the release agent used by OPTICON during repli- 
cation did not deteriorate the micro-roughness of the replicated Al or Au surface. 
The strong decrease of reflectance of the replicated samples without UV coating 
for a wavelength below 200 nm, as apparent in Fig. I4.14[ is caused by the natural 
oxide surface layer as explained in Sec. 14.3.21 Inasmuch as Au and Al coating ex- 
hibit a comparable reflectance, the replication of Al is easier because Au coating 
requires an additional SiO/Cr blocking film to avoid diffusion into the surface UV 
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Figure 4.14: 

Reflectivity of the aluminum 
coating with MgF 2 protec- 
tive layer applied to a CFK 
mirror with aluminum or 
gold replicated coating. 



coating. It can be concluded that the UV reflectance of replicated CFK mirrors 
is not limited by the surface micro-roughness and is comparable to a glass mirror. 

Figure 14.151 shows the local slope error along a 12 mm surface profile mea- 
sured using the experimental setup of Fig. 14.1 HI The surface quality of the 
ARDCO/OPTICON sample is surprisingly poor. The print-through of individual 
fiber bundles was already evident from visual inspection. The deformations have 
a maximum amplitude of the order of 3.5mrad resulting in a slope error of ap- 
proximately 1.6mrad. OPTICON claimed that the large fiber printing was solely 
due to the poor quality of the ARDCO CFK substrate. Rough surface struc- 
tures result in an uneven gap width between substrate and replication master 
and hence in variations of the thickness of the epoxy resin layer. These variations 
induce mechanical stresses when the resin shrinks during the curing process. 

The 1 m diameter SESO mirror |lU4j . originally manufactured for the HADES 
experiment [82J, was made of woven CFK fabric which is very cheap but gives 
inferior stiffness. Its thickness of 2.0 mm corresponds to 0.8% of a radiation 
length. The surface profile shown in Fig. 14.161 is much smoother and exhibits 
a slope error of only 0.11 mrad (rms), one order of magnitude better than the 
ARDCO sample. Detailed study of the mirror shape by HADES |106j revealed 
a large astigmatism due to the lack of stiffness. According to SESO, only an 
increase of substrate thickness up to 4 mm would allow to achieve a sufficiently 
high stiffness. 

The sample provided by COI consisted of a uni-directional pan fiber sub- 
strate (M55J/954-3 made by HEXCEL [T7T7] ) of 0.4 mm in thickness. A novel 
polycyanate ester resin was used to glue the laminate layers for minimal differ- 
ence in CTE between laminate and resin system. An Al surface coating was 
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Figure 4.15: Spatial variation the local slope error (left panel) and average slope 
error assuming spatial invariance (right panel) of a surface profile of the ARDCO 
mirror sample. 

replicated onto the substrate by OPTICON using the replication technique de- 
scribed in Fig. 14.111 The measurement of the slope error along a surface profile 
is summarized in Fig. 14.171 The slope error of 0.1 mrad is remarkably small and 
comparable to that of the SESO sample. 

In conclusion, an average slope error of less than 0.15 mrad can be achieved 
by combining best available CFK substrates with novel efficient manufacturing 
approaches. This result invalidates the widespread opinion that CFK mirrors 
have poor quality due to fiber print-through. 

COI Inc. was chosen for further mirror development based on the measurement 
and its excellent reputation as the largest manufacturer for advanced carbon-fiber 
substrates for all kinds of optical and spacecraft high-tech applications pi] . 



4.4.3 Mechanical stability of the mirror 

Generally, very thin large-size CFK mirrors possess insufficient strength and stiff- 
ness to hold proper shape under their mass unless mounted or otherwise sup- 
ported. Fixing the mirror at its outer perimeter to a mount, as done for the 
RICH-1 mirror, is clearly not possible due to the segmentation of the RICH-2 
mirror. 

There are two concepts for stabilizing a thin mirror shell. First, the sandwich 
design: the CFK face skin bonded to a lightweight core material such as RO- 
HACELL foam and CFK back skin to counterbalance thermal stresses. Second, 
the isogrid design |108j : the thin mirror shell is supported with an isotropic core 
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Figure 4.16: Spatial variation the local slope error (left panel) and average slope 
error assuming spatial invariance (right panel) of a surface profile of the SESO 
mirror sample. 



■a 

1 



o 

K 
u 

u 
o 




10 
[mm] 




0.2 0.4 
[mrad] 



Figure 4.17: Spatial variation the local slope error (left panel) and average slope 
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Figure 4.18: Isogrid design support structure of the second COI prototype mirror. 

structure achieved by assembling strips of CFK flatstock in an egg-crate fashion 
forming equilateral triangle core cells (see Fig. I4.18|) . 

The sandwich approach was first adapted to RICH mirrors by the PHENIX 
collaboration j^J. In Sec. 14.4.21 measurements of a test sample were described. 
The main advantages for RICH mirror application are: very low density foam 
allowing a high core thickness and resulting in an excellent mirror stiffness, low 
material and assembly cost, and technological feasibility for mounting at the back 
skin. An unsealed ROHACELL foam core, used for the PHENIX detector, was 
ruled out because ROHACELL has a strong affinity to take up moisture also 
causing the foam to expand. As a result, an additional sealing of the panel edges 
is necessary to prevent outgassing into the radiator. Replacing ROHACELL by 
a reticulated vitreous carbon core (RVC) with similar properties but less affinity 
to moisture would reduce the outgassing rate but sealing of the core would still 
be necessary. 

The development of the isogrid design was driven by the need for more ther- 
mally stable satellite reflectors for use in telecommunication and in remote sensing 
of the atmosphere. From the fabrication point of view, isogrid mirrors provide 
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higher bending stiffness and are less prone to thermal distortion because the en- 
tire structure is constructed from a single material. An exhaustive comparison of 
sandwich and isogrid design can be found in |108| . For the application to RICH 
mirrors the main advantage of this solution is that the CFK face skin can be as 
thin as 0.76 mm assuming 12 layers of pan fiber tape or about 0.5% in terms of 
a radiation length. Additionally, 1% to 5% of the electrons would traverse the 
material of the support structure made of ribs with 2.5 cm height and 0.76 mm 
thickness. This would correspond to a thickness of 0.6% to 3.8% of a radiation 
length depending on the incident angle of the electrons. From general consider- 
ations it becomes clear that this solution is favorable compared to an equivalent 
(thicker) homogeneous mirror because most electrons traverse the thinner face 
skin and remain within the narrow line width of the u- or 0-meson whereas the 
electrons traversing the much thicker back structure are anyway "lost" in the 
bremsstrahlung tail. The higher number of almost unaffected electrons improves 
the background discrimination considerably. In this design without a foam core, 
there is no need to worry about outgassing or special coating arrangements which 
leads to substantial time and cost savings. The critical issue of this design is the 
possible print-through of the support structure caused by the gluing of the back 
structure to the face sheet. 

The apparent advantages of the isogrid design led to the decision to further 
pursue this novel solution to ultralightweight UV mirrors. 

4.5 Prototype measurements 

The first full size prototype was fabricated by COI and delivered in October 1999. 
The substrate consists of 12 layers of uni-directional pan fiber tape (M55 J/954-3 
by HEXCEL |107j ) with a combined thickness of 0.76 mm. The specific substrate 
was selected to provide a near zero CTE to insure dimensional stability over a wide 
temperature range and for its exceptionally high stiffness. The unitape layers are 
rotated consecutively by 30 degrees each to obtain best thermal isotropy. To 
bond the layers, a polycyanate ester based resin system was chosen to match 
the CTE of the substrate. The face skin restrained to a glass master was cured 
in an autoclave at a temperature of approximately 140 °C. In the same step, 
the mirror surface was replicated by use of a newly developed release agent of 
COI making any additional replication with a metal coating unnecessary. This 
improvement resulted in considerable time and cost savings. After curing, the rib 
structure was bonded to the back side of the face skin which was still restrained 
to the glass master. The isogrid structure of the second prototype including five 
additional vertical ribs is illustrated in Fig. 14.181 After the second curing at room 
temperature to bond the rib structure, the composite mirror is released from the 
glass master. From now on the back structure is responsible for preserving the 
shape of the mirror. 
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Figure 4.19: Setup for the optical mirror quality measurement with a HeNe laser 
(632 nm). 



The reflectance of the replicated mirror surface was sufficiently high even 
without an additional metal coating to allow for laser measurement of the mirror. 
Figure 14.191 shows a schematic view of experimental setup used for mirror shape 
measurements. The mirror was illuminated by a point source of monochromatic 
light (632 nm HeNe laser beam focused by a microscope lens) placed at the center 
of the mirror sphere. The radius of curvature of the mirror panel was measured 
as 7880 ± 50 mm. 

A summary of the optical measurements is presented in Fig. 14.201 It shows 
the focal image of the point source when the whole panel is illuminated. The 
apparent structures indicate rather large surface distortions. The slope error 
is about 1.4mrad (rms). This is equivalent to an image size of about 2.5cm 
(diameter). Heating the mirror up to the operating temperature of 35° C did not 
improve the slope error. A closer examination by illuminating individual surface 
regions of 10 cm diameter, the size of a Cherenkov light cone, reveals a large 
spatial variation of the focal point as evident by the position of the individual 
focal images. This astigmatism is produced by a larger radius of curvature along 
the fiber direction of the front CFK layers and can be directly attributed to the 
higher bending stiffness along the fiber direction. 

The average slope error of the spots, shown as the diameter of the circles, 
varies between 0.40 mrad and 0.91 mrad (rms). It turned out that the surface 
area within each back structure triangle is especially distorted. Most likely it 
is due to stresses induced by the shrinking of the adhesive used to bond the 
isogrid. The direct print-through of the back structure was visible but fairly 
small compared to those other distortions. Although the optical quality for small 
areas was encouraging, deficits in manufacturing and design were apparent. 

Based on this experience the following steps were taken to improve the local 
and the overall mirror shape while increasing the thickness by as little as possible. 

• 50% increase of number of layers for front panel to gain more resistance 
against plucking stresses 
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Figure 4.20: Variation of position and width of the focal image of the first COI 
prototype overlayed with the focal image of the full mirror. 
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• 100% increase of rib structure height to improve overall bending stiffness 

• addition of five perpendicular ribs to the isogrid as shown in Fig. 14.181 to 
enhance spatial isotropy 

• reduction of adhesive used for bonding to decrease plucking 

The second prototype was delivered by COI in September 2000. Figure 14.211 
shows the focal image of the full illuminated mirror using the setup shown in 
Fig. 14.191 The most striking feature of Fig. 14.211 is the absence of substructures 
like the ones in Fig. 14.201 The overall slope error was determined to 0.16 mrad 
by the width of the image projections shown in Fig. 14.221 It is worth stressing 
that design improvements resulted in tenfold increase of the mirror quality. 

The spot size is already close to the limit of the setup defined by the divergence 
of the laser beam and the resolution of the CCD camera. The remaining surface 
distortions can only be seen in the Foucault image of the mirror surface shown 
in Fig. I4.2HI The image was captured on a white screen placed at a distance 
of about 5 m from the mirror. The brighter diagonal lines indicate a remaining 
print-through of the long continuous isogrid rib possessing the same magnitude 
as the local plucking distortions. 

To evaluate the local surface errors quantitatively, the mirror was measured 
with phase shifting interferometry by COI using a CO2 laser with a wavelength of 
10.6 /nil. The phase shift was observed with respect to the position of a reference 
mirror. Figure 14.241 shows the interference pattern for 15 horizontal and for 15 
vertical fringes. The mirror surface structure shown in Fig. 14.251 was reconstructed 
by deconvolution of both interference signals. 

A surface irregularity of 4.2/xm(rms) was measured which is equivalent to a 
peak-to-valley deviation of 25.8 fim. Besides, a small edge cusping at the corners 
was evident resulting in first-order astigmatism of 4 /im peak deviation for an 
average angular profile. It is worth stressing that the central 160-mm-radius zone 
exhibits a very flat average radial profile with less than 1 /xm peak deviation. 

4.6 Summary and outlook 

By use of the Monte Carlo detector simulation based on a realistic hadronic cock- 
tail it was demonstrated that the meson spectroscopy capabilities of the CERES 
detector system would greatly benefit from a new ultrathin RICH-2 mirror re- 
placing the original thick glass mirror. The properties and technical requirement 
of a such a mirror were specified. 

A survey of available mirror substrates led to the conclusion that only mirrors 
made of carbon fibers could fulfill these expectations. The evaluation of CFK mir- 
ror samples from different vendors proved that high UV reflectivity and excellent 
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Figure 4.21: Focal image of the second COI prototype mirror. 
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igure 4.22: Projections of the focal image of the second COI prototype. 
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Figure 4.23: Surface image of the second COI prototype mirror captured with a 
HeNe laser illumination |109| . White spots indicate convex and dark spot concave 
distortions. 

imaging quality, as such the most critical requirements, can be achieved by the 
appropriate choice of carbon-fiber substrate and mirror fabrication technology. 

A novel approach to the stabilization of ultrathin CFK mirror shells by an 
isogrid support structure was developed. It is especially suitable for meson spec- 
troscopy because the relative fraction of dielectron pairs remaining in the reso- 
nance peak compared to those shifted in mass due to bremsstrahlung is enhanced 
and, thereby, the meson signal-to-background ratio is improved. 

The first prototype incorporating the new isogrid design was manufactured 
by Composite Optics Inc. in October 1999. Optical measurements revealed large 
local deformations and a significant astigmatism indicating insufficient substrate 
thickness and support structure stiffness. Based on these findings and bound to 
a tight schedule, the CERES collaboration decided that the considerable risks 
involved in the implementation of the mirror for run time in 2000 could not be 
justified by the expected benefits. Furthermore, the ongoing study of mirror 
technologies was to be finalized. 

Several measures were taken in the fabrication of the second prototype to 
optimize the mirror design. As a result, the overall slope error of the mirror 
was reduced to 0.15 mrad, i.e. ten times. The residual astigmatism of 4 /xm is 
negligibly small. Such an excellent optical quality has never been achieved for an 
ultrathin carbon-fiber mirror. The novel replication technique and the simplicity 
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Figure 4.24: Static interferogram for 15 horizontal (left panel) and 15 vertical 
(right panel) fringes across the aperture using a CO2 interferometer 109 . 




Figure 4.25: Contour plot of the surface error profile in microns 109]. Posi- 
tive/negative values represent concave/convex aberration feature, respectively. 
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of the isogrid design result in low manufacturing cost. Thus, significant cost sav- 
ings of more than 50% could be gained in comparison to the HADES RICH mirror 
which is fabricated on the basis of carbon ceramic substrate technology [106 . 

The second prototype arrived only shortly before the main production run in 
Fall 2000 ruling out the implementation of a new full-size mirror in the CERES 
RICH-2 detector. Nevertheless, the development of an ultrathin carbon-fiber 
mirror in only three years is an accomplishment considering that comparable 
projects took more than 5 years until completion j^S]. Future applications are 
to be expected driven by the recent renaissance of RICH detectors in heavy 
ion physics and particle physics [HEJ USE 11101 HHj - A summary of the mirror 
development is in preparation for publication |112j . 
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Chapter 5 



Analysis of Pb-Au collision data 
at 158 GeV/c per nucleon 

5.1 Introduction 

Any experiment attempting to measure low-mass dielectrons in ultrarelativistic 
heavy ion collisions has to deal with major experimental challenges. First, the 
large number of produced particles leads to a high detector occupancy resulting 
in a serious load on the detectors. Distinction of electrons from the 2 to 3 orders 
of magnitude more abundant hadrons requires a detector system with superior 
particle identification capabilities. Only a small fraction of all electrons come from 
nontrivial sources other than 7-conversion and Dalitz decay n° — > e + e _ 7. Finally, 
uncorrelated electrons and positrons originating from a large fraction of partially 
reconstructed pairs form a huge combinatorial background when combined to 
pairs. 

Although the CERES experiment has been designed to detect electrons under 
such conditions, a sophisticated data analysis is a necessity in order to extract a 
statistically significant dielectron signal. 
The data analysis consists of the following steps: 

• calibration of detector raw data 

• reconstruction of hits in each detector 

• combination of hits of all detectors to particle tracks including particle 
identification and momentum determination 

• rejection of accidentally matching track segments and reduction of combi- 
natorial background 

• single-track efficiency correction by means of Monte Carlo detector simula- 
tion 
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• subtraction of the combinatorial background from the pair distribution 

This Chapter is focused on the results of new studies and developments. The 
key points of the chain of analysis will be explained and details will be given where 
necessary for general understanding. A thorough description of the previously 
used detector calibration and hit reconstruction algorithms can be found in S63, 
IfiTH I7T| 1113] . New detector calibrations developed in this paper are included at 
pertinent places. 

5.2 Reasons for re-analysis of the 1996 data set 

The 1996 data set has been analyzed twice before [SU} 1113j . The third analysis 
was motivated by the following reasons. 

First, in previous analyses the signal of correlated electron-positron pairs was 
extracted by subtracting the mass distribution of like-sign pairs from the unlike- 
sign spectrum. This procedure assumes the combinatorial unlike-sign background 
to exactly resemble the like-sign spectrum. Improper subtraction of combinato- 
rial background as a cause of the dielectron enhancement observed was a seri- 
ous concern, since the signal-to-background ratio is very small (i.e. 1 : 13 |113j ). 
Furthermore, signal and background distribution are similar in shape. The en- 
hancement observed could be explained by a 3.6% increase of the background as 
noted in |113j . Such a situation could occur in case of an unrecognized asymme- 
try of the detector with regard to detection and reconstruction of unlike-sign and 
like-sign pairs. 

Second, a smoothing procedure was applied to the like-sign combinatorial 
background to reduce statistical error. In fact, the statistical errors of the indi- 
vidual mass bins were shifted to a systematic error of the integral distribution. 
A certain level of subject iveness resulting from the particular choice of the fit 
function (see Fig. 5.12 in |71j ) increased the difficulties in understanding the 
systematic error of the background subtraction. 

In this paper, these particular questions will be addressed in the discussion 
of the application of the mixed-event technique, an alternative method for the 
construction and subsequent subtraction of the combinatorial background. The 
comparison of the two background subtraction techniques will allow a test and 
validation of the assumptions inherent in both methods. 

The first attempt to construct a mixed-event combinatorial background was 
based on the results of the previous first-stage analysis |113j . Figure I5~T1 shows 
the mixed-event distribution to deviate strongly from the same-event like-sign 
background. In absence of other plausible explanations, a new first-stage analysis 
had to be performed to exclude any artifacts of the hit and track reconstruction 
or of the event selection applied for data reduction. 

Several minor problems of previous analyses were addressed in the process 
of the data re-analysis with the mixed-event technique. Most notably, a very 
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Figure 5.1: Comparison of mixed-event and same-event combinatorial back- 
ground based on the previous raw data analysis [HHJ EE1 EH Hl^j . Mixed-event 
and same-event background differ by up to 25% for mass above 0.8 GeV/c 2 . This 
deviation is an indication for an artifact of the analysis procedure. 



complex background rejection strategy had been used in (7TJ 1113] which was not 
easy to reproduce. The main focus of this paper was placed on the essentials and 
clarity with special emphasis on the rejection strategy. 



5.3 Hit and track reconstruction 

5.3.1 Technical aspects of the raw data analysis 

The first of two steps of analysis, hit and track reconstruction, requires processing 
8000 GB of raw data, collected in the recording of 6.0 ■ 10 7 events. It was carried 
out at the CERN EFF PC-farm using 40 CPUs for two weeks. The general results 
are summarized in Table 15.11 

Recent progress in data storage technology and the increase of available com- 
puting power allowed to loosen the event selection and to store all events with at 
least one electron or high-momentum pion track. Each event was characterized 
by the centrality and by the orientation of the reaction plane used only in the 
hadronic flow analysis by |114j . 
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Analysis stage 


Number of events 


Recorded on tape 
Available events 
Successfully analyzed 


60.000.000 
41.694.200 
40.418.548 


Events with 7r-track 
Events with e-track 
Events with dielectron 


34.523.208 
25.170.135 
2.380.071 



Table 5.1: Results of first-stage raw data analysis. 



5.3.2 SDD-hit reconstruction 

The major change at the first stage of the analysis was the implementation of 
a novel algorithm to reconstruct hits in the SDD developed by |114j . Before 
evaluating the impact of new software on the dielectron analysis, the general 
concept of SDD-hit reconstruction is explained with special emphasis on the 
differences between the old and the new software version. 

SDDs play a crucial role in the reconstruction of event vertex and tracks and 
in the rejection of pairs of close tracks by either resolving those or using the 
deposited-energy information for discrimination. The SDD-hit reconstruction 
software is used to determine the hit position taking into account effects such 
as the noise of the electronics, pulse shape variations, and the saturation of the 
pulse height. Overlapping hits are resolved as far as possible. 

The strategy previously applied was based on the clustering of pulses of ad- 
jacent anodes to hits and determining the hit position by a center of gravity 
method Overlapping hits were split only in anode direction in case of a 

local minimum between adjacent pulses in a hit cluster. As noted in [H], this 
method produces a large fraction of artificially split single hits deteriorating the 
single-track resolution and the close-track rejection power. The new method was 
developed to improve this situation |114| . It employs a Gaussian pulse fit to find 
the hit center as well as a special logic to recognize overlapping hits based on a 
double-Gaussian fit of the pulse shape in time-bin direction and a simple local 
minimum splitting in anode direction. 

A Monte Carlo simulation of the double-hit reconstruction efficiency by [114 
is plotted in Fig. 15.21 With the old software a large excess of reconstructed 
pairs appears in time-bin direction, an indication of artificial splitting of single 
hits. This peak has disappeared after introduction of the new hit reconstruction 
software. 

The double-hit resolution, defined as the distance of reconstruction of both 
hits with 50% probability, is improved by 60% for anode and time-bin direction. 
The latter is verified in this paper by an evaluation of the measured distribution 
of hits in the vicinity of a reconstructed track which gives a double-hit resolution 
of 5.6mrad (see Fig. 15.30 . 
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Figure 5.2: Monte Carlo simulation of the SDD double-hit reconstruction ef- 
ficiency achieved by means of the new (red) and the old (blue) SDD-hit soft- 
ware |114j . 
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Figure 5.3: Double-hit reconstruction efficiency of SDD-1. Local constancy of 
the hit density provided, the distance between a hit belonging to a track and the 
next closest hit can be described by a Gaussian probability distribution. The 
finite double-hit resolution results in a drop in the observed distribution at the 
point of overlap of close hits. Reconstruction efficiency was determined from 
the relative difference between the distance distribution observed and a Gaussian 
distribution fitted to a distance range of non-overlapping hits. The effects of 
artificially split hits were removed by a linear approximation. 
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Figure 5.4: Density distribution of the next closest hits in SDD-1. The old version 
of the SDD software (left panel) splits overlapping hits only in anode direction. 
Hence overlapping rings are not completely resolved and the double-hit resolution 
in time-bin direction is very poor. The new software version (right panel) splits 
overlapping hits in both anode and time-bin direction invoking a double-Gaussian 
fit. The peak shown in yellow and red color levels indicates artificial hit splitting. 

The Monte Carlo simulation of the double-hit reconstruction efficiency does 
not address the problem of artificially split single hits because a split single hit 
could be mistaken for a properly reconstructed double hit. A new approach based 
on experimental data was taken in this paper to resolve this issue. 

Without artificially split hits, the SDD-hit density must be uniform close to 
a hit belonging to an electron track. Only information about the next closest 
hit was stored in the first stage of analysis. The density distribution of the next 
closest hit rapidly decreases with increasing distance to a given track. It can be 
approximated by a Gaussian probability distribution. One part of a hit artificially 
split belongs to a track while the other results in an enhanced density of next 
closest hits. The finite double-hit resolution causes a depletion in the vicinity of 
a hit. The observed density distribution of next closest hits is plotted in Fig. 15.41 
The old software (left panel) shows a huge excess due to artificial splitting at a hit 
distance range of 1 to 2 anodes and 1 to 3 time bins. Overlapping hits in time-bin 
direction were not split which leads to a depleted hit density up to a distance of 5 
time bins. In contrast, the new software shows this area filled resulting from the 
additional splitting in time-bin direction. However, a significant peak indicating 
artificially split single hits is still present. 

The fraction of artificially split hits was estimated for a quantitative study 
by the integrated yield of the observed excess relative to the number of regular 
background hits. Contrasting both software versions in Fig. 15.51 (left panel) shows 
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Figure 5.5: Impact of artificial hit splitting in the SDD. Fraction of artificially 
split hits as a function of the matching quality cut of SDD-1 and SDD-2 (left 
panel). Contrast of the pair efficiency loss caused by artificial hit splitting in the 
SDD for new and old software version (right panel). The efficiency loss estimated 
in data was verified by an overlay Monte Carlo simulation. Considering the differ- 
ence of the SDD-matching cuts applied in the previous analysis (0.9mrad) [113 j 
and this paper (1.3mrad), an additional 8%/16% track/pair reconstruction effi- 
ciency loss was estimated for the new version. 



the total fraction of artificially split hits to be reduced from 25% (old version) to 
19% (SDD-1) and 16% (SDD-2) (new version). 

The artificial hit splitting in both SDDs occurs independently. In particular, 
it is independent of the quality cut applied to the matching of SDD-1 and SDD-2 
hits unless the position of artificially split hits is shifted. Figure 15.51 (left panel) 
reveals the relative fraction of split hits in the new version to decreases rapidly 
with increasing matching quality exactly as one would expect if an increasing 
fraction of artificially split hits were rejected due to poor matching. This feature 
results in an additional pair efficiency loss as illustrated in Fig. 15.51 fright panel). 
Although the new software reduces the fraction of split hits considerably, those 
that are still artificially split exhibit a larger position shift (i.e. comprise smaller 
hit fractions and/or more single-anode hits) which in turn results in a 8%/16% 
track/pair reconstruction efficiency loss compared to the old software. 

The study of this effect is not yet completed to date but is considering the 
following possible explanations: the Gaussian double-hit fit is more sensitive 
to distorted pulse shapes compared to the center of gravity method and the 
parameters controlling the double-hit splitting in anode direction need to be re- 
adjusted to accommodate small variation in the drift-time calibration or the pulse 
shape between adjacent anodes. 



58 



5.3.3 Outline of the tracking strategy 



After completion of the reconstruction of hits from the detector signals, the spa- 
tial hit positions are utilized to construct particle tracks. The tracking strategy 
applied in this analysis closely follows the approach taken in |113j . It encompasses 
the following successive steps: 

• matching of SDD-1 and SDD-2 hits 

• vertex determination as interception of all SDD track segments 

• reconstruction of RICH-1 Cherenkov rings utilizing SDD track segments as 
pointer to select proper ring candidates 

• combination of fitted RICH-1 rings and RICH-2 rings/PD hits to RICH-PD 
track segments taking into account the inhomogeneous B-field 

• final matching of SDD and RICH-PD track segments to complete particle 
tracks 

The particular purpose of this strategy is to aid the ring recognition algorithm 
in both RICH detectors by external tracking information of SDDs and PD. All 
tracking steps were put under strict scrutiny in search for artifacts that could 
lead to a biased track reconstruction. The resulting improvements and other 
important modifications are documented below. The matching of the detectors 
measuring angles (RICH-1/2) and those measuring spatial coordinates (SDD, 
PD) has been changed everywhere. It is now based on solid angles rather than 
polar and azimuthal angles not reflecting spherical symmetry. 

5.3.4 Reconstruction of SDD track segments 

The most likely position of the event vertex is determined by an optimization pro- 
cedure (Robust Vertex Fit The vertex is defined as the spatial point where 
the weighted sum of its squared distance to all trajectories of matching SDD-1 
and SDD-2 hits is minimal. Next, SDD track segments are created. A simple 
predictor pointing to SDD-1 is calculated for each hit in SDD-2 (see Fig. 15. 6|) . 

Previously, a 4° tilt of both SDDs [HE] with respect to the beam axis was 
neglected in the predictor determination. To include this effect, an improved 
predictor method was developed using analytic geometry. Next, a binary-search 
algorithm is applied to find the SDD-1 hit closest to the SDD-2-vertex predictor. 
This algorithm was modified because it was failing in certain rare cases. The size 
of the search window of the best-matching hit was adjusted to always cover the 
maximum matching window of 7mrad. 

To avoid ambiguous rejection of track segments within the limit of the detector 
resolution, tracks sharing a common hit in SDD-2 (so-called reversed VsDD-tracks) 
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SDD-1 SDD-2 



Figure 5.6: Illustration of the SDD-vertex tracking. The vertex is defined as the 
spatial point where the weighted sum of its squared distance to all trajectories 
of matching SDD-1 and SDD-2 hits is minimal. 

were allowed in addition to track segments sharing a common hit in SDD-1 (Vsdd- 
tracks). 

These changes also apply to the second stage of SDD-vertex tracking where 
the z-position of the reconstructed vertex is refined to the exact position of the 
target disk. An optimized x-y vertex position is obtained by minimizing the 
sum residual distance of the initially found SDD-vertex tracks pointing to the 
target. After refinement of the vertex position, an optimized set of SDD-vertex 
track segments is obtained following the steps described above. This procedure 
improves the SDD-vertex pointing resolution by 20%. Further optimization of 
the code implementation results in a 20% reduction of the overall execution time. 

5.3.5 Reconstruction of RICH-PD track segments 

Particle tracking is closely interlocked with electron identification and momentum 
determination as described in Sec. 13.51 because of the particular setup of the 
CERES detector. Electron track reconstruction in the RICH-PD detector system 
therefore has to meet all of the following principal objectives: highest possible 
tracking efficiency, sufficiently good pointing accuracy for minimal probability 
of accidental matches with fake rings, and precise momentum determination. 
Independent fulfillment of these requirements is restricted by: the momentum 
dependent multiple scattering (see Eq. l4.ip . the 6*-deflection caused by the second- 
order-field effect (see Fig. 13.51 in Sec. 13. 5 J) , and small deviations of the magnetic 
field from the nominal B(r) ~ 1/r dependence. 

Best reconstruction efficiency for low-momentum tracks is vital to reduce the 
combinatorial background. Higher-order effects are important and are to be 
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Figure 5.7: 

GEANT simulation of the 
difference between the (f>- 
measurement in PD and 
RICH-2 depending on the radial 
track positions R. The average 
deflection measured in the 
PD was subtracted: A$ = 

A0SDD-PD — O.64A0 SDD „ RICH _2- 

The functional dependence 
was fitted by a fourth-order 
polynomial to obtain the new 
PD-0-predictor function. 

The previous tracking strategy of 
[BUI IB!)! I7T| 1113] was replaced by a new tracking strategy based on a detailed 
simulation of the CERES detector including all higher-order effects. It grounds 
in several basic ideas as follow. First, all previously fitted RICH-1 rings and all 
RICH-2-ring candidates are combined provided they fall within a butterfly-shaped 
matching window described by: 

A0(RICH-1- RICH-2) = ^resolution + (^scattering Jp) 2 , (5.1) 

accounting for multiple scattering and detector resolution. Moreover, 0-deflection 
in the magnetic field has to be less than 0.9 rad (corresponding to p > 150 MeV/c). 

The expected ^-deflection in RICH-2, caused by the second-order-field effect, 
is calculated for each combination of rings as a function of the 0-defiection and 
the radial position of the track. The latter dependencies reflect the residual non- 
linear contribution of the magnetic field. The radial shift of the vertex with 
respect to the radial symmetry of the magnetic field is also included. 

After this correction, the remaining difference in ^-direction between RICH-1 
and RICH-2 is attributed to multiple scattering and detector resolution. A follow- 
your-nose approach was applied to achieve maximum tracking efficiency With 
multiple scattering occurring mainly in the RICH-1 mirror, a correlated shift in all 
downstream detectors, namely PD and RICH-2, must follow. Multiple scattering 
in azimuthal direction however cannot be distinguished from the 0-deflection by 
the magnetic field and is directly contributing to the momentum resolution. Next, 
a predictor is computed into the PD. It includes multiple scattering in ^-direction 
and the second-order- field effect. The predictor function relating RICH-2 and 
PD coordinates is depicted in Fig. 15.71 It was obtained by tracking 10 6 particles 
through the CERES detector setup using the GEANT simulation package [88J. 




0.1 0.2 0.3 0.4 

[rad] 

<j) deflection (RICH-2) 



treated carefully especially for these tracks. 
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New predictor (o" = 1.5 mrad) 
Old predictor (o" = 3.6 mrad) 



Figure 5.8: 
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Deviation of new and old 
predicted 0-position from the 
true hit position simulated in 
GEANT including multiple 
scattering. The A0 distribution 
of the old version is not only 
much broader but also exhibits a 
general offset of about —4 mrad 
with respect to the true hit 
position. 



The expected ring center position in the RICH-2 detector was determined from 
the average orientation of the particles momentum vector while traversing the 
RICH-2 radiator volume rather than from the center of the fitted RICH-2 ring. 
The latter method becomes inaccurate for low momentum tracks because the ring 
shape observed is distorted by the second-order-field effect. 

Figure 15.81 illustrates the accuracy of the new PD phi-predictor A$pd com- 
pared to the previously employed version. The A0 distribution of the old version 
is not only much broader but also exhibits a general offset of about —4 mrad 
with respect to the true hit position. The matching of PD hits and RICH- 1/2 
segments was little affected because the 0-offset was partly counterbalanced by 
a 3 cm misalignment of the PD z-position. 

In the next tracking step, all RICH-2 ring candidates are fitted, provided a 
PD hit was found close to its predictor. All complete RICH-PD track segment are 
stored for subsequent matching to the previously constructed SDD-vertex track 
segments. 

The maximum search window for the closest PD hit was modified to a fixed 
solid angle (as opposed to the formerly fixed rectangular matching window in x-y 
pad coordinates which does not reflect the appropriate symmetry of the detector) . 

Sometimes more than one track segment shares either the same RICH-1 or 
RICH-2 ring. The first case represents unresolved dielectrons (called V-tracks) 
with a small opening angle for that the segments have opposite charge. They 
are kept for further studies. In previous analyses, all other multiple matches 
were rejected by choosing the track with the best #-match between RICH-1 and 
RICH-2 detector. This procedure favors reconstruction of high-momentum tracks 
which are less affected by multiple scattering. The charge determination is refined 
at a later stage of the analysis, based on the more precise measurement of the 
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Figure 5.9: 

GEANT simulation of 
the relative error of the 
momentum determination 
excluding multiple scat- 
tering. The momentum 
reconstructed with the old 
software exhibits an overall 
2.5% offset which increases 
for small momenta. 



momentum p 



deflection between SDD and PD instead of the initial combination of RICH-1 
and RICH-2 detector. The pair-charge dependence of this rejection introduces a 
subtle bias for high-momentum tracks. 

With the new tracking, all multiple matches were kept for later evaluation 
based on an improved deflection determination and a further rejection of fake 
rings by additional quality cuts. The surviving V-tracks of the same charge were 
used to evaluate the fraction of accidental matches of unlike-sign V-tracks. The 
few multiple matches that remained unresolved after the quality cut were finally 
rejected to avoid a tedious selection procedure. 

By means of GEANT detector simulation, the relation between the momen- 
tum of the electron and its azimuthal deflection in the magnetic field was de- 
termined including higher-order corrections. Figure 15.91 shows the momentum 
resolution contrasting new and old version of the momentum fit. The relative 
error of the new method of 1% is small compared to the observed resolution (see 
Sec. I5.4.3j) . The 2.5% momentum offset for the old version is caused by a previ- 
ously unnoticed misalignment of the PD z-position. This shift was independently 
confirmed by a measurement of reconstructed D-meson mass |116j . 

Which of the software modifications finally led to the convergence of same- 
event and mixed-event background (see Sec. I5.6.4J) could not be verified because 
the necessary repetitions of first-stage analysis were prohibited by their excessive 
need of computing power, storage space, and time. 
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5.4 Global calibration of the spectrometer 



5.4.1 Intercalibration of detectors 

High-precision alignment of all detectors is of crucial importance for efficient 
tracking and subsequent rejection of accidental matches. The external distance 
measurements of individual detector with respect to a fixed laboratory frame, 
made during assembly, do not have sufficient precision. But they provide ini- 
tial values and constraints for an intercalibration with a data sample of recon- 
structed high-momentum pions. Compared to electrons, high-momentum pions 
(p > 4.5 GeV/c) are no much affected by multiple scattering and, therefore, more 
suitable for precise detector alignment. Starting with the autocalibration of the 
SDDs and the reconstruction of the event vertex, described in Sec. 15.3.41 all iden- 
tified pion tracks are matched to the independent reference of PD hits and RICH 
rings. Each detector was aligned in x-y coordinates with respect to the center 
of SDD-1 which was used as reference point. Any residual offsets in ^-direction 
can be removed by applying a radial correction to the local drift velocity of the 
SDDs. 

Using this calibration strategy, a high-precision alignment of the '96 detector 
setup was carried out prior to this paper. Details of the calibration are doc- 
umented in [BUI inn] ■ The attempt to reuse the calibration parameters for this 
analysis failed, as they were unable to reconstruct the vertex at the nominal tar- 
get positions. An investigation revealed the positions of the SDD hits, obtained 
with the new fitting algorithm, to have been subject to a systematic shift due 
to the non-Gaussian pulse shape. The intercalibration of SDD-detector system 
had to be redone. A simple study of the calibration parameters was carried out 
to avoid an elaborate recalibration of the Silicon-vertex telescope. The correct 
spacing of the target disks was regained by decreasing the outer-radius parameter 
of SDD-1 by 24 /zm resulting in a systematic expansion of the radial scale. The 
residual time- dependent offset between the nominal target and the reconstructed 
vertex position was used to compute a small correction factor for the local drift 
velocity in both SDD detectors for each run. After this correction, the recon- 
structed vertex distributions were centered at the nominal position of the target 
disks as shown in Fig. 15.101 (left panel). The individual vertex distributions are 
clearly separated. This allows for unequivocal identification of the target disk in 
which the interaction took place. The SDD-track reconstruction can be refined 
with the knowledge of the exact vertex z-position. 

The higher double-track resolution of the new SDD-hit reconstruction software 
(see Fig. 15 .2j) improved the resolution of the vertex z-position by 20 /im to 260 /im 
(compare [HH])- For geometry reasons explained in |117j the vertex z- resolution 
is about 15% better at the far side of the SDD as seen in Fig. 15.101 (right panel). 

The vertex z-resolution could be further improved by 15% in case of best cal- 
ibration of the SDD- vertex telescope |117j . Figure 15.111 shows the time variation 
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Figure 5.10: Reconstructed vertex z-positions for run 230 (left panel). The vertex 
distributions of the individual target disks have a width of o = 0.24-0.28 mm 
(right panel). The width is small compared to the separation of adjacent targets of 
about 2.4 mm. The vertex z-resolution decreases with increasing distance between 
target and SDDs for geometry reasons. 



of the mean and the width of the reconstructed-vertex distribution for each tar- 
get disk. The residual offset of about 100 /zm leads to a very small systematic 
shift of the radial(#-) matching of 0.15 mrad which is negligible compared to the 
resolution of the SDDs. 



5.4.2 Matching distributions 

Internal consistency and quality of the readjusted calibration can be evaluated 
by looking at the residual offset of the centroid, the shape, and the width of the 
matching distributions of all detector combinations. 

Figure 15.121 shows the matching distributions of high-momentum electrons 
(p>2GeV/c). The large background fraction of the SDD matching distribution 
(first column in Fig. I5.12j) can be attributed mainly to single anode hits with 
very poor resolution in anode direction. In case of a SDD-vertex track segment 
with a single anode hit in either drift chamber, the (^-coordinate of the combined 
track was determined by the 0- value of the other detector rather than the average 
of both drift chambers in order to improve the matching with the downstream 
detectors. The limited reconstruction capability for partially overlapping rings in 
RICH-1 originating from SDD-1 conversions and close ir° Dalitz pairs (here not 
rejected by other means) results in large fractions of unrecognized double rings 
which exhibit a poor matching quality of SDD and RICH-1 (second column in 
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Figure 5.11: Run-to-run stability of the vertex reconstruction. The a-width of the 
vertex fit is indicated by the red dots. The few runs exhibiting a larger cr-width 
contain less than the average number of events per run. The residual difference 
between reconstructed and nominal vertex z-position is small compared to the 
width of the distribution and is therefore neglected. 
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Figure 5.12: Matching of SDD-1 and SDD-2 (first column), combined SDDs and 
PD (second column), RICH-1 and PD/RICH-2 (third column) of high-momentum 
electrons (p>2GeV/c). The ^-positions of hits/rings in PD/RICH-2 were cor- 
rected for the second-order- field effect. The matching distributions were fitted 
with the sum of two Gaussians representing the distributions of true- and fake- 
tracks. The fake tracks were rejected in a later stage of analysis. 

Fig. EH). 

The centroid offset is a measure of the quality of the geometrical inter- 
calibration of individual detectors. The observed offset of less than 0.5mrad 
for all detector combinations is small compared to the width of the matching 
distributions and, thus, confirms the excellent quality of the calibration. 

The remaining small miscalibration of the SDD-vertex telescope previously 
mentioned resulted in a run-to-run variation of the centroid of SDD-RICH-1 and 
SDD-PD matching distribution shown in Fig. l5.lHl The similarity of both distri- 
butions is prove that the variation is indeed caused by the SDD. The offset of the 
matching was not corrected because the SDD-RICH-1 offset is small compared 
to the width of the distribution (see Fig. I5.12|) and the SDD-PD ^-match was 
not used in the analysis (see App. lBj) . 

In the limit of high momentum, multiple scattering is negligible and, hence, 



67 



1 1 1 



t I... 


1 1 

.' 1 


1 1 1 

'. 1 ■■■ 


i 1 i - 

* : ^SDD~^RICH-1 : : 
• ■ - 

1 i 1: 








a : n n 


; w '■ U SDD "PD 


»• • » ■** 






^..M-fcC^....*...: 


y'f 




*. .; 











-0.5 



150 



200 



250 300 
run number 



350 



400 



Figure 5.13: Run-to-run variation of the centroid of SDD-PD and SDD-RICH 
^-matching distribution. The similarity of both distributions is prove that the 
variation is caused by a small miscalibration of the SDD-vertex telescope. The 
residual SDD-RICH-1 offset is negligible compared to the width of the matching 
distribution of 2.1mrad. The SDD-PD #-match was not used for rejection (see 

App.ng. 



the width of the matching distribution a is determined by the contributions of 
the single detector resolutions: 



0"Match Detector 1-2 — \j ^resolution Det 1 + ^resolution Det 2 ■ (5-2) 

Given the matching distributions of any combination of three independent detec- 
tors, Eq. 15.21 can be resolved to extract the single detector resolution: 



. /^Match Detector 1-3 ^Match Detector 2-3 ^Match Detector 1-2 /r n\ 
O resolution Det 1 — \ « • \ a -^J 



Table 15.21 summarizes the extracted detector resolutions. The observed resolu- 
tions are close to the expected values calculated from detector properties [113, 

EH]. 
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Detector 


resolution [mrad] 


SDD-1 


0.28 


SDD-2 


0.28 


RICH-1 


2.03 


RICH-2 


1.66 


PD 


0.54 



Table 5.2: Single-track detector resolutions extracted from the width of the 6- 
matching distributions for various detector combinations. The values given for 
the RICH-1 and the RICH-2 detector represent the ring center resolution. 



5.4.3 Momentum resolution 

The experimental momentum resolution is determined by the accuracy of the 
measurement of the azimuthal deflection between detectors before (SDD and 
RICH-1) and after the B-field (RICH-2, PD): 

dp d((j)sBD, RICH-1 — 0RICH-2,Pd) /r . N 

_ = ~A 7 • ( 5 - 4 ) 

P <PSDD, RICH-1 — </>RICH-2,PD 

It is composed of the single detector resolution and the momentum dependent 
multiple scattering. The latter can be inferred from the experimentally accessible 
^-matching distribution: 



A0(p) « A0(p) = y/o-Lolution + Scattering / ■ (5-5) 

It should be noted that Eq. 15.51 slightly underestimates the multiple scattering 
contribution because it is partly absorbed in the follow-your-nose tracking ap- 
proach applied to the RICH-PD track segments in ^-direction (see Sec. I5.3.5|) . 

The left panel of Fig. 15.141 shows the width of the 6— matching distribution 
of various detector combinations as a function of momentum. To obtain the 
momentum resolution, Eq. 15.41 can be expressed in terms of Eq. 15.51 



dp . . A/ ^resolution P ^"scattering , , 

— (P) = T~7 • (5-6) 

P A0 Dct i_Dct2 

The accuracy of the momentum determination according to Eq. 15.61 is depicted 
in the right panel of Fig. 15.141 for various detector combinations. Examination 
reveals the momentum measurement based on the SDD-PD combination to be 
best for high momenta while the RICH-l-RICH-2 deflection measurement is best 
for low momenta. This ground in the fact that RICH-2 measures the local 0-angle 
of the particle trajectory in the RICH radiator after the azimuthal deflection in 
the magnetic field. The 0-defiection observed in the RICH-2 detector is about 
56% larger than in the PD detector. For small deflection, i.e. high momentum, 
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Figure 5.14: Momentum dependence of the ^-matching of various detector combi- 
nations (left panel). A contamination of high momentum pions and 7 conversions 
leads to a seeming deterioration of the resolution for high momenta. Momentum 
resolution of various combinations of detectors before and behind the magnetic 
field (right panel). 

the resolution is dominated by the intrinsic detector resolution and favors the 
SDD-PD combination. To optimize the mass resolution, the latter effect was 
taken into account by simply switching from the RICH-l-RICH-2 momentum 
measurement to the SDD-PD combination for momenta smaller than 0.8 GeV/c. 
The combined momentum resolution given by: 

^(p) = y/(0.022 ± 0.001) 2 ■ p 2 ■ (GeV/c)" 2 + (0.041 ± 0.002) 2 (5.7) 

agrees with the result of previous studies |113j . 

5.5 dN^/drj measurement and centrality deter- 
mination 

Global observables of a relativistic heavy ion collision such as the multiplicity, 
i.e. total number of emitted particles N^, or the transverse energy E t carry 
important information about the reaction dynamics [1 19|. HL20 . The centrality of a 
collision in particular can be inferred from the particle yield. The characterization 
of collisions in terms of centrality and iV c h forms the basis for comparison among 
various collision systems and different experiments. 

In the CERES experiment, the number of charged-particle tracks iVch is mea- 
sured with SDDs in the pseudorapidity range from 2 to 3. The distribution of N c ^ 
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Multiplicity distribution for 
SDD tracks. 
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per event obtained is shown in Fig. 15.151 The cutoff towards lower multiplicity 
corresponds to the trigger threshold of 100 hits in the multiplicity detector MD. 
The peak at iV c h = 60 is caused by non-target interaction. These event are later 
removed by an offline centrality cut. 

Figure IS". 161 shows the position of the trigger slope (left edge) as well as of the 
central slope (right edge) of the multiplicity distribution to change with time. 
The apparent decrease of the central-slope position was caused by the gradual 
deterioration of one high- voltage sector in the SDD (see Fig. 15.371 in Sec. 15.7.2)) . 
The remaining variation was attributed to a slight temperature dependence of the 
SDD reconstruction efficiency. The data was corrected for both effects. The large 
fluctuation of the trigger slope is most likely caused by an unstable voltage supply 
for the multiplicity detector. It results in a variation of the initial centrality 
selection and, therefore, does not directly impact on the measured multiplicity 
except for the weighted multiplicity average (less than 1%). The limited statistics 
of the data sample did not allow to enable full exclusion of the low-multiplicity 
range affected by the trigger fluctuation. 

The measured multiplicity distribution needs to be corrected for the recon- 
struction efficiency of the SDD to obtain the true number of emitted charged 
particles. The efficiency correction was derived from a Monte Carlo simulation of 
10000 realistic UrQMD events |121j . It describes passage of all particles through 
target and detectors taking into account energy loss, 7-conversions, ^-radiation, 
particle decays, and all detector properties including electronic noise [113) 1114] , 
All simulated events were reconstructed as genuine raw data events. 

Delta electron tracks do not point to the vertex and can be recognized and 
removed to the level of less than one track per event [122 . The beam pile-up was 
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Figure 5.16: Run-to-run variation of the position of the trigger slope (lower panel) 
and the central slope (upper panel). The position of the slope was defined as iV meas 
at half of the maximum (see Fig. 15.15]! . 

considered to be small. 

The comparison of the number of reconstructed tracks with the number of 
initial tracks, illustrated in Fig. 15.171 fleft panel), shows the obtained correction 
factor of 5% to be smaller then the 20% upward scaling found in the previous 
analysis |113j . The difference is attributed to a use of a relatively large matching 
window between SDD-1 and SDD-2 (7mrad) and the inclusion of reversed Vsdd- 
tracks (i.e. two tracks sharing the same hit in SDD-2) in the SDD tracking which 
results in a 15% increase of random background matches as shown in Fig. 15.171 
(right panel). 

The average charge-particle multiplicity obtained with the Monte Carlo cor- 
rection for the two cases with and without reversed VsDD-tracks differs by 10% 
which suggests that Monte Carlo simulation does not fully describe these subtle 
differences in the tracking. Therefore, it seemed best to account for the re- 
versed VsDD-track contribution by applying an additional —15% correction to 
the measured multiplicity. The fully corrected multiplicity distribution is shown 
in Fig. E3E1 

The systematic error of this method is very difficult to evaluate, as the cor- 
rection relies solely upon the Monte Carlo simulation. A rough estimate can be 
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Figure 5.17: Charged-particle reconstruction efficiency simulation (left panel). 
The true number of charged particles is a linear function of the mea- 
sured track multiplicity N meas : N ch = (1.05±0.01) • iV meas + (8. ±2.). The efficiency 
correction of iV meas is larger for the old version without reversed VsDD-tracks 
(^V"ch = 1.13 ■ iV meas + 15.). Impact of reverse VsDD-tracks on the measured aver- 
age track multiplicity as a function of the SDD-l-SDD-2 matching cut (right 
panel). Reversed VsoD-tracks are a source of fake tracks for a matching window 
larger than 2 mrad. 
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Figure 5.18: Efficiency corrected distribution. The multiplicity distribution 
was divided into four bins of equal statistics for the study of the centrality de- 
pendence of the dilepton production. 
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obtained by attributing all variations in the measurement of various track 
reconstruction methods to the systematic error. A comparison of the new and old 
SDD software gives a difference of 5%. Additionally, the Monte Carlo simulation 
lacks a description of the observed efficiency variations (4%), the VsDD-track con- 
tributions (8%), and realistic SDD pulse shapes including artificial hit splitting. 
The sum of all contributions gives an estimate of the upper limit of about 11% 
for the relative systematic error. 

The centrality of a collision refers to the fraction X of the total geometric 
cross section a P 



' geom • 



Cgeom — [R projectile + -R target) j R ~ 1-2 fm^S , (5-8) 

where R and A denote the radius and the mass number of the colliding ions, 
respectively. The total cross section observed can be related to the detected 
number of charged particles for a minimum bias event selection: 

-^V minbias / r n \ 

a = . 5.9 

AT AT v 1 

-''target 1 'beam counter 

Here iV target , iV beamcounter , and iV minbias are the number of target nuclei per unit 
area, the number of beam particles, and the number of observed minimum bias 
events, respectively. Integration of Eq. 15.91 gives the centrality X for a certain 
iV c h cutoff: 

X= —F, ijf) dN ^- ^ 

Cgeom J (cutoff) \dl\ ch J 

The resulting relation is plotted in Fig. 15.191 The N ch cutoff of 104 applied in the 
off-line analysis corresponds to a relative cross section of (32±6)%. The selection 
of most central events according to iV c h > 310 equals a relative cross section of 
(3.5±0.6)%. The relative cross section exhibits a large uncertainty because the 
radius of the colliding nuclei in Eq. 15.81 is not precisely known. 



5.6 Simulation of the combinatorial background 

5.6.1 Sources of combinatorial background 

There are two dominant sources of dilepton in the mass region of m ce < 2 GeV / c 2 : 
first, leptonic and semi-leptonic decays of scalar and pseudo-scalar mesons ■§] and, 
second, photon conversions in the target and the first silicon drift chamber with 
7T° — * 77 being the predominant source of photons. 

These decays lead to the production of electron-positron pairs, henceforth re- 
ferred to as unlike-sign pairs. Production of electron-electron or positron-position 
pairs, so-called like-sign pairs, requires higher-order processes. The strongest of 
them is the 7r — > e + e + e~e~ decay, which is not only suppressed by a factor of 
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~ 380 [Hj relative to the 7r Dalitz decay but also charge symmetric. It is negligible 
for this analysis. 

The experiment measures a certain number of positron tracks n + and of elec- 
tron tracks n_ in each event. Exclusive measurement is not possible because most 
dileptons are produced within the collision zone. Considering all combinations 
of observed electron and positron tracks of an event, it is impossible to decide 
which unlike-sign pair originates from a single decay and which is an accidental 
combination of individual tracks of separate decay processes. 

Therefore, two classes of unlike-sign pairs can be distinguished: the actual 
signal of correlated dielectrons S + _ and the so-called combinatorial background 

pairs The total observed unlike-sign pair distribution can be expressed 

as: 

jytotal = ^corr. + £Comb. _ (5 n ) 

Both signal and background exhibit a continuous spectrum. The combinatorial 
unlike-sign background can be estimated by the same-event like-sign method or 
by the mixed-event technique. 

The first method is based on the fact that the same-event combinatorial like- 
sign background is identical to the combinatorial unlike-sign background in the 
absence of correlated like-sign pairs from physics origin and under the assumption 
of acceptance and efficiency for electrons and positrons being the same. In the 
mixed-event technique, tracks from different events are combined to yield the 
combinatorial unlike-sign background. 

All published dielectron invariant mass spectra of the CERES collabora- 
tion [SZ1 1123[ I124[ I125[ I126j were obtained with the same-event like-sign back- 
ground method. Previous attempts to employ a mixed-event background were 
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Figure 5.20: 

Track multiplicity of run 
144. The distributions of 
electrons and positrons 
are almost identical. Both 
are well described by a 
Poissonian distribution with 
a mean value of fi = 0.515. 
The variation of the track 
multiplicity distribution 
within a charged-particle 
multiplicity bin is ne- 
glected [T27j . 



unsuccessful. The first detailed comparative study of both methods will be pre- 
sented in Sec. l5.6.2l and l5~(Hfl The comparison of both simulated backgrounds can 
yield valuable insight, test inherent model assumption, and ultimately increase 
confidence in the background subtraction procedure. 

5.6.2 Same-event combinatorial background 

The multiplicity of electrons N_ or positrons N + produced in a collision can be 
described by a Poisson probability distribution P (see Fig. I5.20J) : 

]y fc 

P(iV +/ _ = k) = -±p exp (-N +/ _) . (5.12) 

The probability B to observe n + /_ tracks out of N + /_ initial particles is dis- 
tributed binomially: 

B{ "+<- = k) = *^~-*)i (1 - '+/->"♦'-' ■ < 5 - 13 > 

The probability to observe a lepton track £ is a product of the probability of a 
particle falling into the acceptance and the reconstruction efficiency depending 
on various single-track selection cuts. e + and e_ are treated separately to account 
for possible charge asymmetries of the acceptance or the reconstruction. Making 
use of Eq. 15.1 HI one obtains the average number of reconstructed tracks: 

= e + N + ,{^ = 6 + (l-s + )N + + 6 2 + N 2 + , (5.14) 
nl = £_iV_ , JnZf = e_(l -£_)iV_ +e 2 _N 2 . 
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In first-order approximation, the mean number of pairs per event with N + positrons 
and N_ electrons is given by: 



N + b(b — \\ i 
n++ = «++ E 2 B ( n + = k ) = 2 K ++ £ l N +( N + ~ X ) > ( 5 - 15 ) 

k=0 " 

— = K __ ^ k ( k ~ X ) B (n_ = k) = -K__e 2 _iV_(iV_ - 1) , 

fc=0 ^ 

N + N- 

— = k.\ E E k B(n + — k) I B(n_ = I) = «_| e + E-N + N_ . 

k=0 1=0 

The factor k denotes two-track efficiency, introduced by physics correlation, de- 
tector, or analysis for each sort of charged pairs. Making use of Eq. 15.151 the 
number of pairs averaged over all events becomes: 

oo 

(n++> = E ^P(N+) (5.16) 

N + =0 

i oo 

= E N+{N + - 1) P{N+) 

Z N+=0 



oo 1 2 

(n__) = E nZZP(N^) = -K^e 2 _(Nl) , (5.17) 





V_=0 



oo oo 



("+-> = E E - K +- e + £-N + AL . (5.18) 

N + =0N-=O 

Equation 15. 181 represents the unknown unlike-sign combinatorial background. By 

comparison of (n++) and (n ) with it becomes obvious that the geometric 

mean of the like-sign background is an excellent approximation of the unlike-sign 
combinatorial background: 



= 2yW) (n-\ 

unlike-sign bg like-sign bg 



K + _s + 6_N + AL = y/K ++ e + N + N_ . (5.19) 

This proposition is fulfilled if the two-track efficiency k is pair charge independent. 
A charge asymmetry of the single-track reconstruction probability e, on the other 
hand, will not alter this result. A more detailed derivation including an extension 
to the case of the simultaneous occurrence of correlated and independent lepton 
sources can be found in |128j . 

For the detector the assumption about the pair charge independence of the 
efficiency is certainly true for large pair opening angles. In this case, the tracks are 
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Figure 5.21: The magnetic field breaks the symmetry of like-sign and unlike-sign 
pairs, for unlike-sign pairs comprise a "cowboy" and a "sailor" configuration. 



well separated in all detectors. It can, however, be contested if the pair-charge 
symmetry of the detector is broken by the magnetic field. While the like-sign 
pairs with a small opening angle remain always close in space, the unlike-sign 
pairs form either a so-called "cowboy" or "sailor" configuration as illustrated in 
Fig. 15.211 Thus, a finite two-track resolution would affect like-sign and unlike-sign 
pairs differently. 

The second disadvantage of the like-sign background estimation method is 



the large statistical error a(^jA(n ++ )(n )) of the simulated combinatorial back- 
ground: 



a(V4(n ++ >(n__>) = 

V 1 \J a (n ++ ) 



(n+ + )0-(n__) 



(5.20) 



a 



(n— ) 
((nuke) 



2<n ++ )^2(n__». 



A Poissonian distribution for the statistical error and statistical independence of 

the (++) and ( ) pair samples were assumed. By use of Eqs. 15.111 and 15.211 

the relative statistical error of the correlated dilepton signal AS + -/ S + - can be 
expressed as: 



AS 4 



a 



N 4 



a 



B+- 



v2 



Or, 



B A 



(5.21) 



\/2(ni ike ) 



(n + _)-2 v /(n ++ )(n__) 
To circumvent the statistical limitations of the same-event like-sign background, 
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Figure 5.22: Ratio of the like-sign background to the smoothed background (left 
panel). The 4 th -order-polynomial correction function for the smoothed like-sign 
background distribution was included. Comparison of the observed like-sign dis- 
tribution and the corrected smoothed background (right panel). All rejection 
cuts were applied. 



a smoothing procedure was applied in previous analyses [ZH I113j . An approxi- 
mation of the functional shape of the same-event combinatorial background was 
obtained by sampling the final 8-, <ft-, and p^-distribution of single tracks and cal- 
culating the invariant mass of each pair of simulated tracks. Using this procedure, 
a background distribution can be obtained that is basically free of statistical er- 
rors. To account for impact parameter dependence, the tracks were divided into 
four multiplicity and four theta bins. Finally, a mass dependent correction factor 
was applied to the smoothed invariant-mass background distribution. It is ex- 
tracted by fitting the ratio of the observed combinatorial like-sign background to 
the smoothed background while keeping the integral of the measured distribution 
as shown in Fig. 15.221 (left panel). The resulting invariant- mass distribution of 
the smoothed background is plotted in the right panel of Fig. 15.221 

Effectively, this procedure hides the statistical error by shifting it to the sys- 
tematic error of the background subtraction. An additional drawback of this 
method results from the fact that any scattering observed between adjacent mass 
bins is solely attributed to statistical fluctuations. This is a rather weak assump- 
tion because artifacts of the RICH-ring reconstruction algorithm are known to 
alter the touching-ring configuration (see Fig. I5.23| . This leads to structures in 
the pair-opening-angle distribution which are directly translated into localized 
variation in the invariant-mass spectrum. Any correlation of this type would be 
masked by the smoothing. 
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Figure 5.23: 

The pair-opening-angle 
distribution shows nearly 
touching or slightly over- 
lapping RICH rings to be 
"pulled" or "pushed" to a 
touching-ring configuration 
with an opening angle of 
60mrad. Additionally, the 
reconstruction efficiency 
increases as the photon hits 
of one ring promote the 
reconstruction of an other 
touching ring. 

To exclude these potentially dangerous effects and to assess the systematic 
errors involved, it would be of utmost importance to independently verify this 
method. 

5.6.3 Mixed- event technique 

Alternatively to the use of the independent like-sign pairs of each event, unlike- 
sign pairs, obtained by combination of opposite charged tracks of different events, 
are inherently independent. This procedure is commonly known as the mixed- 
event technique. The straightforward modification of Eq. 15.151 gives the mean 
number of observed mixed unlike-sign pairs of two randomly selected events A 
and B with initial multiplicity N+i_ and N^_: 

n^ ed = ^2Y, klB ( n + = k ) B ( nB = +J212 klB ( n - = k ) B ( n + = , 

k=0 1=0 k=0 1=0 

(5.22) 

where each track of event A is combined with all tracks of event B with opposite 
charge. Again, the averaging of many pairs of events with different initial particle 
multiplicity yields: 

oo oo 

^^ n mi_xed pA(i)pB( k } 
i=l k=l 

£*Nf NE + e% e*N* N± . (5.23) 

It is important to note the single-track events to be contributing to this average. 
In contrast, the same-event technique requires at least two initial tracks per event. 




I mixed\AB 
\ n -\ — / 
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To obtain the mixed-event background, the mixed unlike-sign distribution (see 
Eq. I5.23J) has to be normalized with the number of mixed-event pairs A^" mixed : 

(n™* cd ) = J^( n t °? ld ) AB = — eZ 7V+ Nl , (5.24) 

^ 'mixed 

where N exp denotes the total number of fully analyzed events including those 
where no electron track was found. Equation 15.241 defines the mixed unlike-sign 
background. 

All attempts to construct the mixed-event unlike-sign background on the basis 
of the previous raw data analysis failed because the single-track events have been 
left out. Furthermore, only events containing at least one dielectron recognized 
either as a so-called Dalitz or open pair were stored for further analysis (see [7T] 
for detailed description) . Inspection of Eq. I5.23l and l5.24l reveals that the complete 
rejection of single-track events together with a signal pair dependent background 
discrimination must lead to a biased event sample and, hence, to a potentially 
distorted mixed background (see Fig. 15 .If) . 

A comparison of Eq. 15.181 and 15 . 241 shows (n™ 1 * ed ) to be equivalent to the much 
sought-after independent unlike-sign background: 

(n + _) = KL xcd ), 



/c+_ e + s_N + AL = e~^e~N + N_ , (5.25) 

provided: 

e + E- = , (5.26) 

k + _ = 1 . (5.27) 

Experimentally, the condition 15.261 can be approximated by restricting the event 
mixing to classes of events with similar properties and subsequent averaging of 
all classes. Technically, all events were divided into sub-samples of 4 multiplicity 
bins and 158 time bins. The temporal sub-samples were found to be important 
because pressure, temperature, and detector parameters changed considerably 
during the 6 weeks of data recording. 

The condition 15 . 271 means an infinite two-track resolution. While this assump- 
tion holds for large pair opening angle (i.e. a pair >80mrad), it is obviously not 
fulfilled for close pairs. For those the correlation factor k + _ drops below one and, 
hence, the mean number of background pairs is overestimated by the mixed-event 

background. Note, that the correlation factor k h should not depend much on 

charged-particle multiplicity because the two-track resolution is dominated by 
the double- ring resolution of RICH-1 and RICH-2 detector, and the observed 
background in both RICH detectors is related to electronics noise and scattered 
beam particles but not to the event multiplicity. 

In general, there are two main advantages of the mixed-event technique. First, 
it is also applicable to the deconvolution of like-sign correlated signals (e.g. open 
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charm detection |129j) which, however, is of no importance for this analysis. 
Second, the statistical error of the mixed-event background distribution can be 
reduced simply by increasing the number of mixed events iV m i xe d in Eq. 15.241 It 
is important that iV m i xe d is chosen in such a way that the probability to select the 
same-event more than once remains small. Otherwise, the result will be hampered 
by the auto-correlation of the mixed events. Consequently, the ultimate increase 
of statistics is limited by the size of the sub-samples used for mixing. 

For technical reasons the ratio of the number of mixed pairs to the number 
of same-event like-sign pairs was fixed rather than the event ratio. A simulation 
based on a toy model confirmed that a pair-mixing ratio of up to n mixe d/niik c = 20 
is safe with respect to any auto-correlation )| 

In case of a small pair-mixing ratio, the statistical error of the mixed back- 
ground a((nf* ed )) can be expressed in terms of the statistical error of the same- 
event like-sign background: 



Compared to the same-event like-sign background, the statistical error of the 
mixed background is reduced by about a factor of 5 for a standard mixing ratio of 
20. For this example the statistical error of the signal is reduced by 30% compared 
to the same-event method. If the background had no statistical uncertainty at 
all, a maximum reduction of 42% could be achieved according to Eq. 15.221 

The mass dependence of the relative statistical error was estimated by the 
statistical fluctuation of the mixed-event background for randomly selected event 
sub-samples. Figure 15.241 shows the average deviation of the mixed-background 
distribution from the mean value. Taking the example of a mixing ratio of 20 
the relative statistical error of the mixed background increases continuously from 
0.5% in the mass range of 0.2 < m cc <1 GeV/c 2 to about 5% at 2 GeV/c 2 . 

According to condition 15. 2fi[ the shape of the mass distribution of the mixed 
unlike-sign background is very sensitive to variations of single-track reconstruc- 
tion efficiency or acceptance with time. It was carefully checked that the proper- 
ties of the single tracks contained in the mixed-event and the same-event selection 
match well. The single-track distributions of same-event and mixed-event back- 
ground agree well with each other, as apparent in Fig. 15.251 The small differences 
in the (^-distribution are restricted to low-efficiency regions in the RICH detec- 
tors and acceptance holes of the SDDs. These can be attributed to the inherent 
averaging of discontinuities by the mixed-event technique. The impact of this 
effect can be neglected because the occupancy is low in these regions. The close 
resemblance of both p^-distributions of mixed and same-event selections, as seen 
in Fig. I5.2fi[ is proof for local efficiency changes in RICH-2 or PD detector not to 
alter the momentum determination. 
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Figure 5.24: 

Estimate of the statistical 
error of the mixed-event 
background for mixing ra- 
tios of 10, 20, and 40. 



5.6.4 Comparison of mixed-event and same-event back- 
ground 

After the construction of the same-event and the mixed-event background in 
Sec. I5.6.2l and l5.6.3l the inherent assumptions of both methods can be verified by 
comparison of the invariant-mass distributions plotted in Fig. 15.271 The agree- 
ment of the shape of both distributions is striking considering that the number 
of background pairs varies by more than 2 order of magnitude in the observed 
mass region. Combining Eqs. 15.191 and 15.241 gives: 

same-event bg mixed bg 

s + s_ ^ k ++ k__N + N_ = —~N + JV_ . (5.29) 
Thus, the observed equivalence requires that: 

e + E-=~~ and / dK++ — _^ (5.30) 

V dm ee dm ee 

According to Eq. 15.301 the mixed-event average of the single-particle detection 
probabilities resembles the same-event average, and the two-track-correlation fac- 
tor k must be very close to one. This is a remarkable result. It proves the 
particular choice of mixing-sub-sample size to be appropriate to ensure sufficient 
temporal stability of the track efficiency and of the pair acceptance with respect 
to pair mass. Note, that the relative error of the mixed- to same-event ratio is 
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Figure 5.25: Comparison of same-event (blue line) and mixed-event (red line) 
6- and ^-distribution of single tracks shown in the upper and the lower panel, 
respectively. 
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Figure 5.26: Comparison of same-event (blue line) and mixed-event (red line) 
Px-distribution of single tracks. 

dominated by the statistical error of the same-event background as expected for 
a relative mixing ratio of 20. 

The 10% drop of the ratio in Fig 15.271 for pairs with mass below 350MeV/c 2 
is caused by the finite two-track resolution of the same-event background. The 
finite spatial resolution of both RICH detectors [7T| I113| and the artifacts of the 
ring reconstruction for touching rings introduce a correlation for pairs with small 
opening angles, as seen Fig. 15.28} that depends on the RICH-ring distance. 

Although the different behaviour of like-sign and unlike-sign pairs in the mag- 
netic field (see Fig. 15.21}) does not allow to directly conclude that: 

K + _(m ee ) = \J K ++ (m cc )K — (m ee ) = 1 , (5.31) 

it is still a very good approximation because the pair efficiency k is dominated 
by the RICH-1 detector due to its comparably low spatial resolution (Tab. l5~2"j) . 
Most important, this effect is limited to the low-mass region. In principle, the 
double-track reconstruction efficiency can be implemented into the mixed-event 
background but this procedure is obstructed by the complex experimental condi- 
tions (e.g. multiple-anode hits in SDDs and overlapping RICH-rings), eluding a 
precise description by the Monte Carlo simulation. Not correcting for the double- 
track efficiency results in a systematic error of the dilepton signal induced by the 
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Figure 5.29: 
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mixed invariant mass distri- 
bution. 



mixed-background subtraction. However, it is limited to 10% because of the large 
signal-to-background ratio (S : B ~ 1) in the low-mass region. 

It is worth mentioning that the mixed like-sign invariant-mass distributions, 

namely (n ++ ) and (n ), agree very well with each other, as apparent from the 

ratio plotted in Fig. 15.291 Rewriting Eq. 15.241 in terms of (++)- and ( )-pairs 

gives the condition for equivalence: 

(n™ 1 ) = (n mixcd ) 



eTe^N + N + = —~N_N_. (5.32) 

Provided that charge symmetry enforces identical initial multiplicity of iV_|_ and 
iV_ (see Sec. 15.6)1 . the average single-track detection probability e must be iden- 
tical for both charges, exactly as one would expect for the 0-symmetry of the 
CERES detector. 

While the arguments above were solely based on the shape of the mixed- 
event invariant-mass distribution, its total normalization is indispensable for a 
correct background subtraction. As noted in |113j . an underestimation of the 
independent background of about 4% would exhaust the strength of the observed 
dilepton excess. As already explained in Sec. I5.d3| the integrated mixed-event 
background is expected to slightly overestimate the same-event background due 
to the two-track efficiency losses. In contrast, the ratio plotted in Fig. 15.301 
shows that the same-event background overshoots the mixed-event background. 
Correcting both distributions for the single track efficiency determined in Sec. 15.81 
is not sufficient to recover the correct total normalization. There is an indication 
that this effect is related to large localized event-to-event reconstruction efficiency 
changes in both RICH detectors. If a continuous pad-wise gain calibration could 
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Figure 5.30: 

Ratio of the mixed-event 
unlike-sign background, nor- 
malized by the number of 
mixed-event pairs, to the 
same-event like-sign back- 
ground for 4 multiplicity 
bins. Both spectra were cor- 
rected for single-track effi- 
ciency 



152. 212. 275. 359. 
^ch'acc. 

resolve this problem, remains to be seen in future analysis. Unfortunately, a 
slight multiplicity dependence of the effect (see Fig. I5.30|) does not allow for 
a sufficiently accurate re-normalization. It should be mentioned that, if this 
effect were an artifact of the same-event background, it would result in a twofold 
increase of pair signal which can be clearly excluded. 

To avoid these difficulties, the normalization of mixed-event background was 
fixed with respect to the total number of observed same-event background pairs 
with mass above 0.35GeV/c 2 . This selection explicitly excludes the low-mass 
range, where the mixed-event background is not identical to the same-event back- 
ground (see Fig. 15.28)1 . As a consequence, 20% of the gain in statistics of the 
mixed-event background had to be sacrificed. The error of the total normaliza- 
tion is now limited by the statistical error of the same-event like-sign background 
integrated above m ee >0.35GeV/c 2 . 

To summarize, it is demonstrated that the mixed-event background resembles 
the same-event like-sign background. Theoretical and experimental consideration 
lead to the conclusion that both accurately simulate the much sought-after inde- 
pendent unlike-sign background. It is worth stressing that the observed likeness 
excludes the existence of any significant pair correlations and, therefore, strongly 
supports the validity of the background subtraction procedure. 
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5.7 Reduction of combinatorial background 



5.7.1 Rejection strategy 

As set out in Sec. 15 .6| the correlated-dilepton signal has to be extracted from 
the observed number of unlike-sign pairs by the subtraction of the independent 
combinatorial background pairs (see Eq. 15.11(1 . In fact, the contributions of the 
7To Dalitz decay, the 7 conversion decay, and the single tracks of partially recon- 
structed pairs would overwhelm the number of signal pairs by three orders of 
magnitude for no further rejection. Moreover, the huge relative statistical error 
of the signal, given by Eq. I5.22| would invalidate any measurement owing to the 
very small number of expected signal pairs. For this reason the recognition and 
subsequent rejection of tracks stemming from the above mentioned sources is of 
utmost importance for the reduction of the combinatorial background. 

Any rejection of tracks based on a certain cut criteria must balance the ob- 
tained background rejection power g and the unavoidable loss of efficiency e. 
Both are related to the relative statistical error of the pair signal by: 



dS VS + 2B cmnh y/2B comh 

T = s w ~ s - (5 - 33) 



V 1 - Q V^iStial 



£ ^initial 

where g denotes the probability to reject a combinatorial background pair and 
£ the pair reconstruction efficiency for a particular set of cuts. According to 
Eq. 15. 34[ it is possible to optimize the statistical significance of the observed signal 
by minimization of the ratio (1 — g)/e 2 with respect to all applied cuts. Carried 
out on data, this procedure has the inherent danger of selecting a statistical 
upward fluctuation [80J. Therefore, the reconstruction efficiency was determined 
independently in a Monte Carlo simulation of the detector system |113j . 

As already noted in Sec. 15 .2\ the rejection strategy chosen in the previous 
analyses |%fH I113j involved several correlated multi-parameter cuts. The resulting 
optimization procedure is very complex and difficult to reproduce. One of the 
goals of this work was to simplify the rejection strategy by focusing on a few 
powerful cuts and to gain better understanding of those. 

It is useful to contemplate the characteristics of 7To-Dalitz and 7-conversion 
decays for the following discussion. The yield of electron tracks originating from 
the 7r decays is given by: 

dN n „ rfALo dN ee o 

BR,o^ e+e - 7 • — f- ■ — - = 0.43 ■ 0.01198 ■ 2. = 1.1 • 10~ 2 . (5.34) 



dN ch " ^ e 1 dN ch dN^o 

Here BR 7r o_ i . e + e - 7 and -^f^- are the branching ratio and the ratio of the expected 
number of 7r mesons to the number of observed charged particles within the 
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acceptance of the CERES spectrometer, respectively. Equation 15. 341 results in an 
average of 1.4 electrons tracks per event for a mean value of N c h = 138 (77 = 2.1- 
2.65). The additional contributions of the 7/0 and rf Dalitz decay are negligible 
(see Sec. EJ. 

All photons created either by meson Dalitz decay or in the initial collision 
can convert into dilepton pairs while traversing, first, the target and, then, the 
downstream detectors. The induced yield is determined by the dominating meson 
decays: 

dN e dN-o dN~ 7 X 7 X X , 

— = 2 ■ — f- ■ — ^ • - ■ — = 2 -0.44 -2.25- - - — = 1.5- — 5.35 

dN ch dN ch dN W0 9 X 9 X X 1 ; 

0.0025 cm „ . 

= "(135cm - ~ (target) 

= 1.5.^^ = 4.5.10-' (SDD-1), 
9.36 cm 

where X/X is the thickness of the material in terms of its radiation length. The 
conversions in the target and in the SDD-1 lead to a mean track- multiplicity of 1.5 
and 0.6, respectively, for a mean value of N c h = 138. Any conversions occurring 
downstream of SDD-2 are rejected by requiring a particle hit in both SDDs. 

Both conversion and Dalitz pairs are distinguished from all other sources by 
their very low mass of m ee < 0.2 GeV/c 2 (see Fig. IA.4I in App. lA"j) . As the low- 
momentum tracks with p± < 60 MeV/c cannot be reconstructed in the magnetic 
field (see Sec. 15.3.5)) . many of those pair are only partially reconstructed and, 
hence, the mass remains unknown. Owing to this, the rejection is primarily 
based on the pair-opening-angle characteristic. Figure IH.31I illustrates that most 
Dalitz and conversion pairs have an opening angle of less than 35 mrad which 
makes it an equally distinguishing feature as the pair mass. Thus, even partially 
reconstructed Dalitz and conversion pairs can still be recognized as such, if an 
additional hit in the SDDs or a ring in the RICH-1 detector was found in the 
close vicinity of the reconstructed track. Figure IH.32I illustrates the notations for 
the common pair configurations. 

In anticipation of the following detailed study of each rejection cut, a sum- 
mary of the complete background rejection strategy is already given in Fig. 15.331 
In contrast to the previous analysis (see Fig. 6.5 in |113j ). the new strategy en- 
compasses only five rejection cuts discarding two, but at the same time improves 
the rejection power by a factor of 1.4. 

The so-called Dalitz cut was discarded for two reasons. First, it rejected little 
background and, second, single tracks were rejected relying on the properties 
of low-mass unlike-sign "signal" pairs which cannot be simulated by the mixed- 
event technique because all mixed pairs are uncorrelated by definition. The cut 
also induces a subtle background correlation explained in the following. All pairs 
of a given event that share a track with an unlike-sign Dalitz pair, identified 
by a pair mass below 150MeV/c 2 and an opening angle of less than 50 mrad, 
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Figure 5.31: Simulated opening-angle distribution of target and SDD-1 conver- 
sions (left panel) and of Dalitz and vector meson decays (right panel). The 
spectrum of conversion pairs was obtained by Monte Carlo simulation of Pb-Au 
collisions created by the URQMD model |121j . The opening-angle distribution 
of pairs originating from Dalitz and vector meson decays was modeled with the 
GENESIS event generator 
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Figure 5.32: Schematic view of the distinguished dielectron configurations. 
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Figure 5.33: Combinatorial background reduction by rejection cuts. The values 
of the previous analysis were included for comparison (see Fig. 6.5 in |f 13j ). 

are rejected. Only 50% of the identified pairs are truly Dalitz pairs, as the 
signal-to-background ratio is about 1:1. A misidentified Dalitz pair rejects one 
unlike-sign and one like-sign pairs in the example of a three-track event. In most 
cases, the unlike-sign pair rejected must have a larger mass than the misidentified 
Dalitz pair. However, no such restriction applies to the rejected like-sign pair. 
As a result, the unlike-sign background is slightly overestimated/underestimated 
by the like-sign background for large/small mass, respectively. For the second 
discarded cut, namely the close-partial-track cut, it was found that the gain in 
rejection power did not offset the decrease in efficiency, thereby, lowering the 
significance of the signal. 

A comprehensive summary of all rejection cuts is presented in Table IB. II in 
App. El Note, that the order of the rejection-cuts is important as most of those 
are partially correlated. 
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5.7.2 Double-d£/d;c rejection in SDD-1 and SDD-2 

The energy loss of electrons E\ oss in the SDDs can be approximated by a Landau 
distribution [Hj: 

(/ ~ A expt-^ + e"*) with x = Eloss - Emax , (5.36) 



dEdx 2 o 

where E max and a denote the most probable energy loss and the width of the 
distribution, respectively For a thin detector the distribution is skewed towards 
high energies (the so-called Landau tail) due to the large fluctuations of the 
number of collisions involving large energy transfer. As apparent in Fig. 15.341 the 
Landau distribution describes the measured data very well, even though the ob- 
served shape is a convolution with various electronics and detector characteristics 
as described in [74"j . 

Dilepton pairs with an opening angle of less than 2.5 mrad cannot be resolved 
in the SDDs and therefore deposite twice the energy of a single track. The corre- 
sponding double-amplitude peak is clearly seen Fig. 15.341 The target conversions 
and close Dalitz pairs distinguished by their small opening angle are most effi- 
ciently rejected by a correlated cut requiring high hit amplitude in both SDDs 
as outlined in Fig. 15.351 In SDD-1, the pair-opening-angle range to be rejected 
was artificially increased by addition of the amplitude of the next closest hit in a 
5 mrad range to further enhance the rejection power. 

In view of the fact that in the previous analysis the double-d-E 1 / dx peak (in 
Fig. I5.35J) was already rejected at the first stage of the data analysis, it might 
be conceivable that important calibrations were casuistical or simply overlooked. 
Most important, the so-called ballistic-deficit correction accounts for the fact 
that the measured number of FADC counts decreases as a function of the radial 
position of the hit. The radial width of the electron cloud increases with the drift 
time due to diffusion according to [132] : 

v 2 r(r)=v 2 rjro + 2Dt dTift (r) , (5.37) 

where a rjro and D are the initial radial width of the hit and the diffusion constant, 
respectively. The response function of the SDD pre-amplifier depends on the 
width of the input signal and, therefore, the measured output signal decreases 
for hits at the lower radius. This effect is further amplified by fact that the 
amplitude threshold, applied in the hit reconstruction for the noise suppression, 
introduces a relative amplitude loss that is largest for hits with a large width. 
This implies that the SDD-hit amplitude has to be re-calibrated inasmuch as the 
amplitude threshold is changed. 

For the previous calibration shown in Fig. 15.361 the mean hit amplitude de- 
creases by as much as 30% towards the outer radius of the SDDs, indicating an 
over-correction of the ballistic deficit and the amplitude threshold effect. As the 
relative change in ^-direction varies also with the absolute amplitude value, the 
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Figure 5.34: Double Landau fit of the dE/dx distribution in SDD-2. The dE/dx 
distribution of single and double hits is indicated by a blue and red line, respec- 
tively. The apparent difference for tracks with dE/dx below 700 FADC counts is 
caused by artificially split hits which are shifted to lower dE/dx (see Sec. 15.3. 2*j) . 




Figure 5.35: SDD-1 vs. SDD-2 dE/dx distribution. The hatched area is rejected 
by the double-d-E/dx cut. 
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Figure 5.36: Most probable number of FADC counts of SDD-hits measured as 
function of the radial(0-) hit position for single and double hits. To measure 
dE/ dx correctly, the number of FADC counts needs to be corrected for the ballistic 
deficit and the amplitude threshold effect. However, a significant decrease of 
the dE/ dx measured is observed towards large radii after applying the correction 
function of the previous analysis |113j (blue line). The constants of dE/dx desired 
was regained after the recalibration (red line). 

new ^-dependent calibration function was interpolated for hits with amplitudes 
between those of single- and double-amplitude hits. An additional correction 
factor was applied for each anode to account for small amplitude variations il- 
lustrated in Fig. 15.371 Small temperature variations of the SDD system during 
the 6 weeks of data acquisition altered the observed mean amplitude by up to 
2% |133j and were corrected for each run. 

As a result of the recalibration, the width of the Landau distribution decreased 
by about 30%, thereby improving the separation of single- and double-d-E/da; 
distribution. 

To optimize the cut values of the double-di?/ da; rejection, a simple method 
solely based on the track properties was developed. In contrast to the Monte 
Carlo simulation used in |113j . this approach sidesteps the difficulties involved 
in the simulation of the very complex SDD characteristics. Assuming that the 
amplitude distributions of single- and double-amplitude hits were exactly known, 
Eq. 15.341 could be employed to find the optimum values of the dE'/dx-contour cut 
as shown in Fig. 15.351 First, the probability to observe a track with single or 
double amplitude is defined as: 

jytrack jytrack 
Psin S lc = J^H aIld Pdoublc = jytrack" • ( 5 - 38 ) 
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Figure 5.37: Anode-wise recalibration of the dE/dx distribution. A hardware 
problem in SDD-2 appearing in run 308 resulted in an increasing deterioration of 
the dE/dx measurement of 20 anodes (red line). 



Applying the cut, the probability of a single track surviving that cut defines 
the efficiency e and the rejection power r is given by the probability to reject a 
double-di?/ dx track: 



^sYngt( survivin g) 

7ytrack 
single 



and r 



double (rejected) 



7ytrack 
double 



(5.39) 



The probability that a signal pair survives the cut can be expressed in terms of 
the single-track efficiency as: 



pee _ 2 p2 
signal fc single " 



(5.40) 



All combinatorial pairs containing at least one double-amplitude track are re- 
garded as background. The probability to find such a pair is given by: 



pee 
background 



1 r) 2 -PhouMp + 2 £ (1 r) P S i ne i e P, 



single -* double 



(5.41) 
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Figure 5.38: 2-dimensional-double-Landau fit of dE/dx distribution. 

Expressing Eq. 15.341 in terms of Eqs. 15.381 IS~4T)l and 15.411 yields the optimization 
function: 

m ^), m J^ + U0-)\^0\ . „ 42) 

\ O / \ \ iV doublc/ £ ly double/ 

Next, the observed dE/dx distribution of Fig. 15.351 was fitted with a two-dimen- 
sional double-Landau distribution. The fit result is shown in Fig. 15.381 The 
wide peak at a double amplitude in SDD-2 and about single amplitude in SDD- 
1, seen in the distribution of the residual difference between the fit and the 
data, is obviously related to 7 conversion decays in SDD-1. These were excluded 
from the fit to be treated in a separate rejection cut. Although the fit function 
underestimates the tails of the data distribution, it was made sure, that the most 
important region between the peaks of both distributions is described accurately 
enough and the residual does not affect the optimization. 

The contours of cut values with equal efficiency and equal rejection are plot- 
ted in Fig. 15.391 The value of the optimization function (see Eq. I5.42j) also 
shown in Fig. 15.391 peaks at the cut values of dE/dxsDD-i > 1130 ± 20 and 
dE/dx SDD _ 2 > 1240 ±20 corresponding to a rejection power of about 95% and an 
efficiency of 99%, respectively. The maximum is relatively broad and, therefore, 
the significance of the signal must be insensitive to slight variations of the cut 
values. Applying the dE/dx rejection cut to the data improves the signal-to- 
background ratio by a factor of two. 
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Figure 5.39: Efficiency and rejection contours of the double-d-E 1 / dx cut overlaid on 
the optimization function. The value of the optimization function (see Eq. 15.42)1 
illustrated by color levels peaks at the cut values of d-E/dxsDD-i > 1130 ±20 and 
dE/dxsDD-2 > 1240 ±20 corresponding to a rejection power of about 95% and an 
efficiency of 99%, respectively. 

5.7.3 SDD-1 conversion rejection 

The main characteristics of SDD-1 7-conversions are the following: a single- 
amplitude hit in SDD-1, a double- amplitude hit in SDD-2 or a second hit in the 
close vicinity, and an overlapping- or double-ring in RICH-1 detector. Therefore, 
the SDD-1 conversion cut rejects tracks with a double amplitude in SDD-2 (in- 
cluding a summation of the amplitude of the next closest hit within 7.5 mrad) and 
a sum amplitude in RICH-1 larger than 1550 counts. An investigation revealed 
that the rejection power of this cut is limited mainly by the poor separation of 
isolated- and overlapping-rings in RICH-1 detector (only 50% see Fig. 15.40)1 . In 
parts, this can be attributed to the fact that the selected double rings are indeed 
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Figure 5.40: RICH-1 sum-amplitude of isolated- and overlapping-rings. Isolated 
rings were selected by requiring a single-track dE/ dx value in both SDDs and no 
V-track signature in the RICH detectors (see Fig. I5.32j) . V-tracks exhibiting a 
double-track dE/ dx value in both SDDs were regarded as overlapping rings. 

only partly overlapping ring with an opening angle of less than about 10 mrad. 
Depending on the actual center of the reconstructed ring, a certain fraction of the 
overlapping rings will not be covered by the summation mask, which is a ring with 
about 6 mrad width. A simple calculation shows the maximum sum-amplitude 
separation between isolated- and overlapping-rings to be about a factor of 1.7. 

The mean amplitude per ring varies by 25% with time as apparent in Fig. 15.411 
indicating that the gain of RICH-1 detector was not always properly readjusted 
to account for variations of atmospheric pressure. The mean number of hits per 
ring depicted in Fig. 15.411 is an alternative observable to distinguish isolated from 
overlapping rings. It is much less sensitive to gain variations but the separation 
of isolated- and overlapping-rings is equally poor. If a correlated cut including 
the number of hits per ring and the amplitude sum shown in Fig. 15.421 would 
improve the situation remains to be seen. The above described version of the 
SDD-1 conversion cut rejects about 60% of the like-sign background. 

5.7.4 Track quality 

Several track quality criteria help to reject so-called fake tracks reconstructed 
from accidentally matching hits and rings. A detailed description can be found 
in [Oni EH! H13| - The cut values chosen in this analysis are largely identical to the 
previous analyses except for those affected by the improved SDD resolution. 

Most important, the SDD-l-SDD-2 matching cut was refined to 1.3 mrad 
including the dependence on z-position of the event vertex |117j . The data plotted 
in the left panel of Fig. 15.431 shows the 15% resolution increase expected between 
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Figure 5.41: Run-to-run variation of the most probable number of hits per ring 
and the most probable ring sum-amplitude of isolated- and overlapping-rings in 
the RICH-1 detector indicated by blue and red points, respectively. Isolated 
rings were selected by requiring a single-track dE/ dx value in both SDDs and no 
V-track signature in the RICH detectors (see Fig. !5.H2jl . V-tracks exhibiting a 
double-track dE/dx value in both SDDs were regarded as overlapping rings. 
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Figure 5.42: 

RICH-1 sum-amplitude vs. 
mean number of hits per ring 
of isolated- and overlapping- 
rings. Isolated rings were se- 
lected by requiring a single- 
track dE/dx value in both 
SDDs and no V-track signa- 
ture in the RICH detectors 
(see Fig. I5.r>2j) . V-tracks ex- 
hibiting a double-track dE/ dx 
value in both SDDs were re- 
garded as overlapping rings. 
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Figure 5.43: Width of the matching distribution of SDD-1 and SDD-2 as a func- 
tion of the target disk (left panel). Matching distribution of SDD-1 and SDD-2 
of single- and double-anode hits (right panel). All tracks with a SDD- matching 
quality of less than 1.3mrad were rejected. 



the first and the last target disk. It is worth stressing that the SDD-matching 
cut plays an important role in the rejection of target conversions and close Dalitz 
decays because it acts like an opening-angle cut for partially reconstructed pairs 
(i.e. one SDD hit missed). The great disadvantage of this cut is to reject a 
substantial fraction of tracks comprising single-anode hits which exhibit a very 
poor matching resolution (right panel of Fig. I5.43j) . 

The rejection of displaced artificially split hits results in a pair-efficiency loss 
of about 16% which cannot be recovered without sacrificing rejection power. As 
a result, the SDD-matching cut had to be loosened by 50% compared to the 
previous analysis (0.9mrad) in order to maximize the statistical significance of 
the open-pair signal. 

Additionally, misidentified charged pions contaminating the sample should be 
rejected. Considering that only pions with a momentum of more than 4.5 GeV/c 
produce rings in the RICH detectors, high-momentum pions are rejected by the 
characteristics of a small deflection in the magnetic field in combination with a 
smaller ring radius compared to electrons. Figure 15.441 shows that the misidenti- 
fied pions can be clearly distinguish from high-momentum electron tracks. The 
outline of the applied cut is overlayed in Fig. 15.441 

It was noticed for the first time that the nominal ring radius in both RICH 
detectors changes by about 0.8% due the variation of the atmospheric pressure as 
illustrated in Fig. 15.451 The ratio of the measured ring radius to the nominal ring 
radius of both RICH detectors is added for the rejection cut, thereby amplifying 
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Figure 5.44: Rejection of misidentified high-momentum pions (see App. iBj) . 
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Figure 5.45: Run-to-run variation of RICH-ring radius overlayed with the relative 
change of the atmospheric pressure. 



102 



Detector 


0-mm [mrad] 


6-max [mrad] 


77-min 


77-min 


SDD-1 


90. 


300. 


1.86 


3.1 


SDD-2 


90. 


300. 


1.86 


3.1 


RICH-1 


141. 


290. 


1.90 


2.65 


RICH-2 


141. 


240. 


2.1 


2.65 


PD 


135. 


252. 


2.05 


2.69 


Combined 


141. 


240. 


2.1 


2.65 



Table 5.3: Summary of the acceptance of all CERES detectors. The value cited 
for RICH-1/2 detector corresponds to the 2/3 ring acceptance. The second-order- 
field effect restricts the inner edge of the RICH-2 acceptance for pj_ < 150 MeV/c. 

the pressure dependence. If not corrected, this effect reduces the rejections power 
of the cut by about 12% and the reconstruction efficiency by 1.3%. 

5.7.5 Detector acceptance 

Table 15.31 summarizes the geometrical acceptance of the individual detectors. 
Even though the final pair acceptance is limited by the RICH-2 detector to 
the pseudo-rapidity range of 2.1 to 2.65 (corresponding to a #-range of 141 to 
240 mrad), the larger acceptance of the SDDs and of the RICH-1 detector is use- 
ful for the rejection of conversion and Dalitz pairs that fall only partially into 
the final acceptance. It is worth mentioning, that the second-order-field effect 
deflects tracks traversing the magnetic field towards the beam line, i.e. lower 
theta (see Fig. 13.51 in Sec. 13. 5|) . The resulting momentum-dependent restriction 
of the pair acceptance becomes significant for tracks with a momentum of less 
than 150MeV/c. 

The track reconstruction limits the acceptance to a momentum range of 
0.17 <p < 9 GeV/c. The lower limit is imposed by the second-order-field effect, 
which distorts the RICH-2 rings, and the rapidly rising probability to pick up 
accidental matches. The high-momentum limit results from two effects: first, 
the charge determination becomes ambiguous for very small deflection in the 
magnetic field because of the finite detector resolution and, second, the particle 
identification is lost due the contamination of high-momentum pions. 

A further restriction of the low-transverse-momentum acceptance is the most 
powerful tool to reduce the combinatorial background. The transverse-momentum 
distribution of the most interesting vector-meson decays peaks at about 350 MeV/c 
while the distribution of the trivial Dalitz decays rises exponentially for small mo- 
menta as depicted in Fig. 15.461 

The generally applied transverse-momentum cut of > 200 MeV/c rejects 
more than 60% of the background at a cost of less than 10% efficiency of the 
vector mesons. To study the low-momentum aspects of the dielectron spectrum, 
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Figure 5.46: 

Transverse-momentum dis- 
tribution of Dalitz decays 
and vector-meson decays in- 
cluded in the hadronic cock- 
tail simulated by GENE- 
SIS [H|. 



the transverse-momentum cut was lowered to p± > lOOMeV/c. 

Finally, pairs with an opening angle of a e e > 35 mrad were selected to further 
suppress dielectrons stemming from 7-conversion and Dalitz decays. According 
to Fig. 15.311 the opening angle of these background sources is predominately 
below 35 mrad. 



5.8 Reconstruction-efficiency determination 

The probability to observe a collision-produced dilepton is limited by the effi- 
ciency of the detection in each detectors, the track reconstruction, and the back- 
ground reduction. To determine the absolute yield of dileptons, the observed 
number of correlated unlike-sign pairs given by Eq. 15.111 must be corrected with 
the actual reconstruction efficiency e ee of each pair: 

N N 1 . 

" even t pai r -1 

^ee = E E — • (5-43) 

i=l k=l £ ee 

It is worth stressing that the efficiency correction of each pair cannot be generi- 
cally replaced by the average reconstruction efficiency because: 




(5.44) 



In the following it is assumed, that the pair efficiency can be factorized by the 
single-track detection probability and a small correction factor accounting for 
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pair correlations as is indirectly proven in the discussion of the mixed-event back- 
ground in Sec. 15.6.31 In view of the careful calibration of the detector properties 
with respect to the pressure and temperature variations during the run time, the 
residual time-dependence of the reconstruction efficiency is neglected. Therefore, 
the pair reconstruction efficiency can be expressed in terms of the single-track 
phase space parameters (here 9, <fi, and p) and the number of observed charged 
particles which is related to the centrality of the collision: 



- ee 



e e = £ e (N ch ,9,(f),p) 

K ee = K ee (a ee ,N ch ,6i,<j)i,Pi,62,<j)2,P2) ■ 

There are two ways to determine the pair reconstruction efficiency. First, it can 
be estimated by a comparison of expected- and observed-yield of pairs with mass 
below 200MeV/c 2 . The dilepton production in this mass range is dominated by 
the Dalitz decay of ttq-, tj-, and ^'-mesons. The yield of these contributions was 
measured with an accuracy of better than 10% jS3 E3| and is well described by 
the simulation of the cocktail of the hadronic sources including the acceptance 
and the momentum resolution of the CERES detector (see App. IA"]) . Although 
this method does not allow extraction of differential pair efficiency distributions, 
it is reckoned as a valuable reference for the average pair efficiency. 

The second option to determine the pair reconstruction efficiency is a Monte 
Carlo simulation of the complete detector system, including the track reconstruc- 
tion and the background rejection. It involves the following steps: 

• simulation of a huge number of dielectrons with realistic kinetic properties 
according to the decays of the known hadronic sources (App. EJ) 

• use of the GEANT software to simulated the passage of the generated 
pairs through the CERES-detector system. GEANT simulates all par- 
ticle interactions with detector materials such as multiple scattering and 
bremsstrahlung. The hit positions of all particles as well as their energy 
deposition was determined by the response function of the individual detec- 
tors taking internal physics and electronics effects into account. For the first 
time, the emission and subsequent tracking of Cerenkov photons included 
the optical properties of the RICH detectors. 

• embedding of detector responses obtained into the raw data of genuine 
events 

• performance of the first-stage data analysis on so-called overlay events 

• application of all background rejection cuts 
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After the last step, the differential reconstruction efficiency is determined by the 
probability to identify the simulated track in the vicinity of its original position 
in each detector. 

Extensive study of this method had proven that all relevant detector char- 
acteristics are reasonably well described by the simulation |113| . Since then, 
the simulation was refined in many details. Small inconsistencies between the 
GEANT detector emulation and the Monte Carlo simulation could be resolved. 

As a consequence of the modified tracking strategy and the new SDD-hit 
reconstruction, the Monte Carlo simulation had to be readjusted to reproduction 
of the distributions observed in data (see |114j for details). 

A number of 2 • 10 6 simulated tracks were embedded into raw data events on 
a one-to-one basis for the study of the single-track reconstruction efficiency. The 
raw events were chosen from different parts of the run time to average detector 
aberrations. A simulated particle was counted as successfully detected, if all hits 
of a reconstructed track were within a range of three times the detector resolution 
to the particle's original direction. Deciding if a particle was truly lost or merely 
scattered is ambiguous to some extend because genuine and embedded hits cannot 
be distinguished on the detector level. A comparison of events, analyzed with 
and without an additional embedded track, showed this systematic error to be 
in the order of 2%. The statistical error of the expected yield of dielectons being 
large, a relative error of less than 10% is acceptable for the differential shape of 
the efficiency distribution. 

The projections of the multi-dimensional efficiency function e e of Eq. 15 .451 were 
plotted for following discussion. The 0-dependence depicted in the left panel of 
Fig. 15.471 is almost flat, except for a large hole at 1.3 rad. This is caused by a 
region of dead anodes in the SDD-2. The 0-dependence of the efficiency was 
not corrected apart from the influence of the dead anodes, as the complicated 
interplay of several effects is not well understood. 

Figure E33 (right panel) shows the ^-dependence to be roughly constant above 
170mrad but to decreases by more than a factor of two towards the inner edge 
of the acceptance. This efficiency drop reflects the increasing hit/ring density 
close to mid-rapidity at llOmrad and the decreasing number of hits per ring in 
RICH-2 detector. At the inner edge about one third of the area of a RICH-2 ring 
is already outside of the detector acceptance. 

Figure 15.481 (left panel) demonstrates the reconstruction efficiency to be ap- 
proximately independent of momentum for a deflection smaller than 230 mrad 
(i.e. 1/p < 1.6GeV _1 c). It drops rapidly towards larger deflection correspond- 
ing to pj_ < 200 MeV/c, which however is only important for the p± > 100 MeV/c 
selection. The slight decrease towards smaller deflection, i.e. larger momentum, 
can be attributed to the momentum dependence of the butterfly-shaped matching 
window of the tracking between RICH-1, RICH-2, and PD detector. 

The efficiency of the track reconstruction additionally depends strongly on 
the charged-particle multiplicity as shown in the right panel of Fig. 15.481 The 
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Figure 5.47: Single-track reconstruction efficiency as a function of the azimuthal 
angle <fi (left panel) and the polar angle 9 (right panel). The variations in the 
distribution result from dead or inefficient anodes in the SDDs and efficiency 
variations in the RICH detectors caused by discharges and hot spots. The hole 
at 1.3 rad is related to a region of dead anodes in the SDD-2. The ^-dependence 
of the efficiency was not corrected apart from the influence of the dead anodes, as 
the complicated interplay of several effects is not well understood. The efficiency 
drop towards small ^-values as apparent in the right panel reflects the increasing 
hit /ring density close to mid-rapidity at 110 mrad and the decreasing number of 
hits per ring in RICH-2 detector at the inner edge of the acceptance. 



efficiency obtained by the Monte Carlo simulation agrees well with the values 
extracted with the first method, i.e. the ratio of the observed number of low- 
mass pairs to the expected hadronic yield of each multiplicity bin. 

In general, the multi-dimensional efficiency function Eq. 15.451 cannot be de- 
scribed by a simple product of its single-parameter projections (i.e. e(0), s(9), e(p), 
ande(iV c h)) discussed above because some of those may not be independent. In 
the following, the relation between the single-parameter efficiency projections is 
studied to find a simple representation for the reconstruction efficiency. 

Given the finite detector resolution, the reconstruction efficiency of a track is 
subject to the density of close hits or rings in the respective detectors. The hit 
density is determined by multiplicity- and ^-distribution of the charged particles 
produced in the collision, comprising mostly pions. The rapidity density dN ch /dy 
was observed to be constant in the CERES acceptance |134j . The energy of the 
pions produced in Pb-Au collisions at 158 GeV/c per nucleon is much larger than 
their rest mass. Therefore, the rapidity y of a particle can be approximated by 
the pseudo-rapidity rj which in turn depends on only on the polar angle 9 135J. 
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Figure 5.48: Single-track reconstruction efficiency as a function of the inverse 
momentum 1/p (left panel) and the charged-particle multiplicity (right panel). 
The inverse momentum is proportional to the azimuthal deflection A0 in the 
magnetic field between RICH-1 and RICH-2 detector (see Eq. 13. The blue line 
in the left panel shows a polynomial fit to the simulated data to be used later 
for correcting the momentum dependence of the reconstruction efficiency. The 
MC-simulated multiplicity efficiency distribution (blue line) depicted in the right 
panel agrees well with the average single-track efficiency (red points) determined 
by the ratio of the measured low-mass yield to the expected yield of the hadronic 
sources. 

As a result, the local hit density per unit area is given by a simple product of two 
functions of the charged-particle multiplicity and the rapidity, respectively. In 
general, it is not possible to disentangle multiplicity- and rapidity-dependence of 
the reconstruction efficiency accordingly because the efficiency is assumed to be a 
non-linear function of the local hit density. The 2-dimensional contour depicted 
in Fig. 15.491 (left panel) shows the efficiency to drop more rapidly at the inner 
edge of the acceptance. This is more clearly seen in Fig. 15.491 (right panel) - the 
slope of the ^-dependence increases with rising multiplicity. This effect has the 
important consequence of the pair efficiency to be increasing with large opening 
angles because these pairs fall only into the acceptance if both tracks are close 
to the upper ^-limit of acceptance. 

Figure 15.501 (left panel) shows the momentum dependence of the efficiency 
to be correlated weakly with the charged-particle multiplicity. As a result, the 
applied single-track efficiency correction can be decomposed as: 

£ e = e(p)e(6,N ch ) . (5.46) 

To investigate the residual pair efficiency (see Eq. I5.45|) . a number of 1 • 10 6 
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Figure 5.49: Single-track reconstruction efficiency as a function of 6 and SDD 
multiplicity (left panel). Relative change of the ^-dependence of the efficiency 
with SDD multiplicity (right panel). The 9- and the multiplicity-dependence of 
the efficiency are correlated as the slope of ^-dependence increases with rising 
multiplicity. 




Figure 5.50: Single-track reconstruction efficiency as a function of inverse momen- 
tum and SDD multiplicity (left panel). The inverse momentum is proportional to 
the azimuthal deflection A<j) in the magnetic field between RICH-1 and RICH-2 
detector (see Eq. 13.11 Overlay Monte Carlo simulation of pair reconstruction ef- 
ficiency as a function of the opening-angle cut (right panel). The simulated pairs 
consisted of 7r° and rj' mesons. Only pairs with an opening angle larger than a ee 
were accepted. 
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simulated dielectron pairs, produced in ir and 77 Dalitz decays, were embedded 
into genuine events. If two tracks were found according to the criteria mentioned 
above, the pair was counted as successfully reconstructed. The residual pair 
efficiency correction was determined as: 

£ MC 

Kee = .-track 1 .-track 2 ' ( 5 " 47 ) 

where e e was computed by Eq. 15.461 

A residual efficiency correction depending on pair properties is expected to 
depend at most on the pair opening angle for the case of touching or overlapping 
RICH rings. However, this effect shown in the right panel of Fig. 15.501 turns out 
to be small. This was already expected from the discussion of pair correlations 
in the context of the background subtraction. 

The observed yield corrected for Monte Carlo pair efficiency is about 30% 
lower than the expected yield of all hadronic sources of pairs with mass below 
200MeV/c 2 . This discrepancy is not surprising. Although the Monte Carlo 
method allows to study all aspects of the analysis, its benefit to the absolute 
efficiency determination is limited by the multitude and high complexity of the 
involved detector and analysis dependencies. Particularly, the RICH-ring re- 
construction algorithm and the hit reconstruction in the SDD are sensitiv to 
marginal changes of the parameters of the read-out electronics and of the envi- 
ronment conditions and are obstructing an adequate description by the Monte 
Carlo simulation. 

A solution was to correct all data first with the Monte Carlo pair efficiency to 
account for the differential efficiency dependencies and then to normalize the mass 
spectrum to the expected dielectron yield of pairs with mass below 200MeV/c 2 
simulated by GENESIS jHO]. 
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5.9 Physics results and discussion 



5.9.1 Results of the new analysis of the 1996 data set 

Background rejection and subtraction was applied to 4.1 ■ 10 7 recorded events, 
as described in Sec. 15.61 and 15.71 The results were a final data sample com- 
prising 3537±103 dielectrons for m ee <0.2GeV/c 2 and 1305±210 dielectrons for 
m ee > 0.2 GeV/c 2 , reconstructed with a signal-to-background ratio of 1:0.96 and 
1:13.9, respectively. The signal refers to pairs with an opening angle larger than 
35 mrad and a minimum transverse momentum of p± > 200 MeV/ c of both tracks. 
The spectrometer acceptance covers the pseudo-rapidity range of 2.1 < 77 < 2.65. 

The resulting invariant-mass spectrum of Pb-Au collisions at 158GeV/c per 
nucleon is presented in Fig. 15.511 The pair yield plotted has been divided by the 
charged-particle multiplicity in the CERES acceptance. It was measured to be 
N ch = 250±30 (see Sec. l5.5J) for the larger pseudo-rapidity range of 2 < 77 < 3 which 
translates into an average multiplicity of (-/V c h)acc. — 138 ± 16 for the nominal ac- 
ceptance. The spectrum was corrected for single-track efficiency as described in 
Sec. l5.8l and. subsequently, normalized to the expected yield from hadronic decays 
for m cc < 0.2 GeV/c 2 using the GENESIS event generator (see App. lXJ) . Compar- 
ison of observed and expected pair yield results in an average 6% reconstruction 
efficiency of low-mass pairs (m ee < 0.2 GeV/c 2 ), assuming that the total yield is 
saturated by the contribution of known hadronic decays. This pair efficiency is 
30% lower than in the previous analysis |113j . The newly implemented splitting 
of overlapping hits in the SDD results in an unexpected efficiency loss of pairs 
estimated to be 16% (see Sec. l5.3.4|) . The insufficient size of the matching window 
between RICH-1, RICH-2, and PD detector applied in first-stage analysis caused 
an additional pair efficiency loss of about 10% (see Fig. 15.48)) . 

The statistical error of the signal includes the combinatorial background con- 
tribution which corresponds to the like-sign same-event background for pairs with 
mass below 0.35 GeV/c 2 and to the mixed-event background for pairs with mass 
above 0.35GeV/c 2 (see Sec. l5~Q|) . 

Three sources contribute to the systematic uncertainties of the data: the 
extrapolated low-mass yield of the hadronic cocktail, the average charged-particle 
density (N ch ), and the pair reconstruction efficiency. 

Uncertainty in the low-mass yield of the hadronic cocktail arises from error 
of the branching ratios and relative production cross sections of 7r , 77, and 7/, as 
well as the error related to the parameterizations of the p± input distributions 
estimated to be about 25% in App. El 

The error of the average charged-particle density measurement results from 
the following factors: the accuracy of the reconstruction efficiency determined by 
Monte Carlo simulation, the linearity of the rapidity distribution in the range 
of 2 < 77 < 3, the beam pile-up, the production of 5 electrons, and the run-to-run 
variation of the trigger threshold and of the reconstruction efficiency. All other 
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Figure 5.51: Inclusive dielectron-invariant-mass spectrum of Pb-Au collisions at 
158GeV/c per nucleon divided by the number of charged particles. The solid 
line represents the expected yield of all hadronic sources [SU] . 

contributions are absorbed into the error of the relative normalization to the 
cocktail, except the time variations. The latter is estimated to be about 5%. 

The systematic error of the pair reconstruction efficiency is given by the uncer- 
tainty of the Monte Carlo description, namely: the detector response functions, 
the alignment of the overlay tracks to the event vertex, run-to-run variations of 
the gain in the UV detectors, and run-to-run temperature dependent changes of 
the SDD and RICH detector properties. As the data is normalized to the low- 
mass yield of the hadronic cocktail, the uncertainty in the pair reconstruction 
efficiency affects the results via its variation with invariant mass. This contribu- 
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tion is estimated to be about 10%. 

The combined systematic uncertainty of the absolute yield is about 40%. It 
was verified that this value is in accordance with the systematic error estimated 
by the change of the dielectron yield with respect to small variations of the 
rejection cuts. The yields of pairs with mass above m ee > 250 MeV/c 2 obtained 
from independent analyses of the two data sets of opposite field direction agree 
with each other within the limit of the statistical error: N ee = (5.3 ± 3.2) • 10~ 6 
and N ee = (4.8 ± 2.9) ■ 10~ 6 for positive and negative B-field, respectively. 

Following previous analyses, the data is compared to the expected yield of 
hadronic decays as simulated with an improved version of the GENESIS event 
generator (see App.[XJ). The hadronic cocktail was folded with the experimentally 
measured spatial and momentum resolution. Note, that for all previous CERES 
publications, the predicted yield is 30% too low for m ee < 0.2 GeV/c 2 due to a 
heretofore unnoticed computing problem within GENESIS (see App. lAjl . 

The most striking feature of Fig. l5~5Tl is the large excess of observed dielectrons 
with respect to the contributions of the hadronic decays. Starting at an invariant 
mass of about twice the pion mass, the data begins to deviate from the cocktail 
encompassing all the range up to the u resonance. For even higher mass the 
statistical error of the data becomes large, but the observed spectrum seems 
to concur with the expected decay contribution of the 0-meson. Integration of 
the measured yield of pairs with mass above m ee > 250 MeV/c 2 renders a relative 
enhancement by a factor of 3.0±1.3(stat.)±1.2(syst.) with respect to the expected 
hadronic cocktail. 

In |71[ 1113] . a strong increase of the dielecton enhancement was reported 
for small transverse pair momentum (q± < 0.5 GeV/c). However, the statistical 
significance was limited due to the transverse-momentum cut of p± > 200 MeV/c 
imposed on single tracks. To verify and substantiate this remarkable observation, 
the data analysis was extended to very low momentum tracks with p±_ > 100 MeV/c 
for the first time. This extension also made possible the study of a previously 
inaccessible region of phase space. 

The net signal increases dramatically for a lower transverse-momentum cut of 
p ± >100MeV/c: 19212 ± 291(stat.) and 2018 ± 382(stat.) for m cc < 0.2 GeV/c 2 
and m ee > 0.2 GeV/c 2 , respectively. The signal-to-background ratio, however, 
deteriorates by almost a factor of two to 1 : 1.74 and 1 : 36 for low-mass and high- 
mass pairs, respectively. Figure 15.521 shows the efficiency corrected dielectron 
invariant mass spectrum for p± > 100 MeV/c. Compared to the p± > 200 MeV/c 
spectrum, the abundance of low-mass dielectrons, mostly stemming from tt° and 
rj Dalitz decays, increases by a tenfold, while the high-mass region of the spec- 
trum remains little affected. The observed excess relative to the hadronic cocktail 
is similar to the p± > 200 MeV/c selection, as expected for an inclusive spectrum. 
Integration of the measured yield above m ee < 0.250 GeV/c 2 gives a relative en- 
hancement factor of 3.7 ± 1.0(stat.)±1.5(syst.). 

Study of the multiplicity dependence provides insight into the dielectron pro- 
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Figure 5.52: Inclusive dielectron invariant mass spectrum for p± > 0.1 GeV/c di- 
vided by the number of charged particles. The solid line represents the expected 
yield of all hadronic sources PU] . 



duction mechanism. For pairs originating from the decay of hadrons in the final 
state the yield must scale linearly with the number of particles produced upon 
freeze out of the fireball. Dielectron production, resulting from ittt annihilations 
in the hot and dense hadron gas formed in the early stages of nuclear collisions, 
is expected to increases quadratically with particle density [136 : 



dN e( 
drj 



'dN, 



cli 



dri 



with 



a 



2 . 



(5.48 



Other proposed collision scenarios involve a scaling behaviour characterized by 
values of a = 1.1 |137j or a = 1.3 jlHj for the dependence on charged-particle 
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multiplicity. 

The 1996 Pb-Au data sample, encompassing a centrality range of the top 
32% of the geometric cross section, was divided into 4 multiplicity bins of equal 
statistics as shown in Fig. 15. 181 of Sec. 15.51 The upper panel of Fig. 15 .531 shows the 
invariant mass spectrum measured for each multiplicity bin. The dielectron yield 
exhibits a strong increase with multiplicity for the mass range of 200-700 MeV/c 2 
indicating a non-trivial origin of the enhancement. 

To illustrate this dependence more clearly, the relative enhancement factor 
defined as the ratio of the dielectron yield to the hadronic cocktail is plotted in 
the lower panel of Fig. l5.53l for three invariant mass bins. The enhancement factor 
stays constant for pairs with mass below 250MeV/c 2 as expected for dielectrons 
originating from 7r°, 77, and rf Dalitz decays. The linear increase (a ~ 2), apparent 
for the mass range of 250-680 MeV/c 2 , provides strong evidence for two-body 
annihilation processes as the major source of the observed pairs. Again, an 
almost constant enhancement factor is observed for pairs with cu-meson mass and 
above, typical of final-state hadron decays. The mass spectra including very low 
momentum tracks with p± > 100 MeV/c display the same feature (see Fig. I5.54J) . 
even though the conclusion of a non-trivial origin of the enhancement is less 
compelling due to the large statistical errors. 

The dielectron transverse momentum q t , i.e. the total momentum of the pair 
perpendicular to the beam axis of the colliding nuclei, is an additional observable 
for discrimination of different production mechanisms. Figure 15.551 presents a 
comparison of the measured lorentz- invariant q t spectra with the hadronic cocktail 
for three invariant mass bins. Inspection of the mass range of 200-725 MeV/c 2 
reveals that the excess is most pronounced for small transverse pair momenta 
below lGeV/c. Most interesting, the limitation of the pair acceptance by the 

>200MeV/c cut, as apparent from the dip of the cocktail spectrum at small 
q t , is also visible in the data. Given that the excess persists for the extension 
of the acceptance towards smaller transverse momentum (p± > 100 MeV/c), it 
cannot possibly be related to a deficient understanding of the pair acceptance 
of the spectrometer. No significant deviations from the hadronic cocktail are 
observed for pairs with lower/higher invariant mass. 

For an alternative representation the invariant mass spectrum is presented 
separately for transverse pair momentum below and above 500MeV/c. Both 
spectra plotted in Fig. 15.561 are distinctly different. The difference between the 
data and the hadronic cocktail becomes more pronounced for the q t < 500 MeV/c 
selection - now extending up to the 10 resonance, while the excess is largely reduced 
for q t >500MeV/c - but still significant in the mass region of 500-680 Me V/c 2 . 
It is worth stressing that the spectrum for g t <500MeV/c is greatly improved 
by the reduced transverse-momentum cut of p± > lOOMeV/c. Integration of the 

> 200 MeV/c data sample above 250MeV/c 2 invariant mass produces an en- 
hancement factor of 3.6 ± 2.4(stat.) and 1.6 ± 0.8(stat.) for q t < 500MeV/c and 
q t > 500 MeV/ c, respectively. 
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Figure 5.53: iV^-dependence of the invariant mass spectrum for p± > 0.2 GeV/c 
(upper panel). The solid line represents the expected yield of all hadronic 
sources [90 . iV c h-dependence of the enhancement factor for three different in- 
variant mass bins (lower panel). 
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Figure 5.54: iV^-dependence of the invariant mass spectrum for p± > 0.1 GeV/c 
(upper panel). The solid line represents the expected yield of all hadronic 
sources [90 . N^- dependence of the enhancement factor for three different mass 
bins (lower panel). 
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Figure 5.55: Transverse-pair-momentum spectra for p±_ > 200 MeV/c (upper 
panel) and p± > lOOMeV/c (lower panel). The solid line represents the expected 
yield of all hadronic sources [HQ] , 
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Figure 5.56: Transverse-pair-momentum dependence of the invariant-mass spec- 
trum for p± >200MeV/c (upper panel) and p± > lOOMeV/c (lower panel). The 
solid line represents the expected yield of all hadronic sources |5D] . 
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5.9.2 Comparison to other CERES results 



In this section, the results presented in Sec. 15.9. II will be discussed and compared 
to the results of previous CERES studies of Pb-Au collisions at 158GeV/c per 
nucleon [ZE1 EDI EHl EEHE] • 

The most prominent observation of a strongly enhanced dielectron produc- 
tion for invariant mass above 200MeV/c 2 was confirmed by this re-analysis. A 
collation of all measurements in Fig. 15.571 shows the observed dielectron yields to 
be consistent within the statistical errors. The important discovery of a stronger 
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Figure 5.57: 

The observed yield [HJ I126j 
exceeds the expected contri- 
butions of hadronic sources 
by ca. a factor of 2.5 for 
m ee >200MeV/c 2 . As ex- 
pected, the dielectron yield 
measured at different stages 
of the analysis does not 
depend on the signal-to- 
background ratio. 
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than linear rise of the dielectron production rate with multiplicity was substan- 
tiated for the mass range of the largest enhancement (0.25-0.725 GeV/c 2 ). The 
level of agreement is remarkable given the fact that this new analysis was based on 
a completely different background subtraction technique, an improved tracking 
strategy, and a refined efficiency determination. 

The following detailed comparison is focused on the latest and most advanced 
previous study |113| which also served as a starting point of this paper. Subject- 
ing the invariant mass spectra plotted in Fig. 15.581 to direct comparison reveals 
differences that need to be addressed. 

First, the spectra differ by about 35% for low-mass pairs (m ee < 200 MeV/c 2 ). 
This is because the efficiency used to correct the data in the previous analysis 
was obtained from a Monte Carlo simulation while in this study the spectrum 
was normalized to the expected yield at the ir° peak (m ee < 200 MeV/c 2 ). The 
resulting discrepancy is still within the range of the systematic error of both 
measurements. The results for pairs with larger mass agree well with each other 
except for the data points at m ee = 400(650) MeV/c 2 . The apparent differences 
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Figure 5.58: Comparison of the invariant-mass spectrum (p± >200MeV/c) with 
the results of previous analyses by the Heidelberg group (HD) [7TJ I113j . The 
solid line represents the expected yield of all hadronic sources |90j . 

are statistically significant. The previous analysis presented in and this 

paper are based on the same data set; therefore, the observed discrepancy is a 
measure for the systematic error of the analysis. 

Second, the statistical errors plotted are larger for the new analysis because of 
the lower reconstruction efficiency, but similar to those of the old analysis before 
applying the background smoothing which involved a certain level of subjective- 
ness in the choice of the best background fit function. 

The factorization of the and 9 efficiency dependence and the 2.5% offset 
in the momentum determination result in a systematic error of the old analysis 
which is small compared to statistical uncertainties. 
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Figure 5.59: Comparison of the (^-spectra for p± >200MeV/c with the previous 
analysis (upper panel). Comparison of the invariant mass spectrum for 

>200MeV/c with previous analyses by the Heidelberg group (HD) [7TJ 1113] 
for q t < 500 MeV/c and q t > 500 MeV/c (lower panel). 



Figure 15.591 (upper panel) compares the transverse pair momentum spectra 
for p± > 200 MeV/c. In contrast to the previous analysis, an excess is already 
observed for a pair momentum of about lGeV/c in the mass region of 200- 
725 MeV/c 2 which is even more pronounced for the p± > 100 MeV/ c sample. The 
difference is most likely an artifact of the factorization of the N c ^- and ^-efficiency 
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dependence. In central collisions, where the largest excess is observed, factorizing 
underestimates the slope of the 9 efficiency dependence (see Fig. 15.49)) . Therefore, 
tracks at small theta are suppressed and tracks at large theta are enhanced. 
Consequently, the pairs contributing to the excess are artificially enhanced for 
very low transverse pair momentum and suppressed for larger transverse pair 
momentum. 

The same feature can be seen in the comparison of the invariant mass spectra 
for transverse pair momentum smaller and larger than 500MeV/c in Fig. 15.591 
(lower panel). The enhancement observed in |113j is larger (smaller) than what 
was found in this study for q t < 0.5 GeV/c (q t > 0.5 GeV/c), respectively. 

5.9.3 Theoretical interpretations 

In Sec. 15.9. 1| it was demonstrated that the hadronic cocktail does not suffice 
to explain the observed dielectron yield. The stronger than linear rise of the 
enhancement with the number of charged particles in the final state points to 
a two-body annihilation process. The high abundance of pions at SPS energies 
makes pion annihilation 7r + 7r~ — > p — > e + e~ the most likely explanation. In a hot 
and dense hadronic medium, the annihilation process is subject to modifications 
induced by interactions with surrounding hadrons and/or partial restoration of 
chiral symmetry as set out in Sec. 12.11 The system spends most of the time 
in its hadronic phase, even though evidence was found for a phase transition in 
Pb-Au collisions to a quark-gluon plasma [2 . The expected contribution from 
quark-quark annihilation is very small compared to conventional sources of di- 
electrons |139j . 

A complete description of dilepton production in heavy ion collisions requires 
modeling of the time evolution of the collision system. In general, theoretical col- 
lision models can be divided into hydrodynamical approaches |140[ I14l( I142| [143 J . 
transport models |144[ I145| PUT)} I147| 1148] . and thermal fireball models jl49|H50| 
115 H 1152"] . The main advantage of hydrodynamical simulations is the capability 
of incorporating phase transitions in a well defined way via the equation of state. 
In contrast, transport models are better at the implementation of rescattering 
and absorption processes. The phenomenological fireball models allow for simple 
comparison of underlying microscopic models. 

To explain the CERES dielectron data of nucleus-nucleus collisions, various 
options were proposed including Brown- Rho scaling [THl I153| I154| I155j . colli- 
sion broadening of the p-meson spectral function |15fi| 1149] , open charm produc- 
tion |157j . chiral meson mixing |158| 1159"! llfiOj . quark-quark annihilation [148] . 
and thermal plasma radiation |161j . A recent review can be found in {21]. 

The recent version of the thermal fireball simulation by Rapp [162] was chosen 
as a representative model to compare to the experimental results. It is generally 
recognized as a comprehensive and reliable simulation of the dilepton produc- 
tion in heavy ion collisions. Several scenarios of in-medium modifications for 
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invariant pair mass 

Figure 5.60: Comparison of the invariant mass spectrum for p± > 0.2 GeV/ c with 
theoretical models fLW2\. The dielectron yield predicted for nir annihilation was 
added to the standard cocktail of the hadronic sources (without p-meson contri- 
bution). Experimental acceptance as well as momentum resolution were applied 
to the model calculations. 

the 7T7T annihilation process are among the specifics of the model used. The 
time evolution is treated in a thermal fireball approach. The experimentally 
determined initial conditions (T ini = 190 MeV, ft n i = 2. 55 go), the hadro-chemical 
freezeout (Tf Q = 115 MeV, Q{ = 0.33^o), as well as a finite pion chemical potential 
fi w are included. 

Figure l5~.fi0l compares the inclusive invariant mass spectrum for p± > 0.2 GeV/ c 
with three different theoretical scenarios. First, the irir annihilation with vacuum 
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spectral function gives too much yield at the p/uj-peak and hardly fills the hole 
between m ee = 0.2-0.7 GeV/c 2 . It can be concluded that in-medium modification 
must play an important role. 

Second, the dropping p mass scenario according to Brown-Rho scaling (or 
Hatsuda-Lee sum rules) fits much better to experimental data but underesti- 
mates the observed yield at the peak of the free vacuum p-meson. This scenario 
entails a reduction of the p-meson width as well as a sharp threshold at twice 
the pion mass for the onset of the enhancement. Brown-Rho scaling is based 
on phenomenological implementation of the restoration of chiral symmetry in 
the framework of an effective field theory. In this case, the dependence of the 
in-medium p mass on temperature T and density q-q is given by: 



with C = 0.15, T} = 200 MeV, and a = 0.3 (QCD sum rule estimate). 

Third, both 7r and p properties are modified in the medium due to rescatter- 
ing (collisional broadening of the spectral function). The resulting spectrum is 
very similar to the dropping p mass scenario for mass below 600MeV/c 2 . How- 
ever, more strength is expected at the vacuum p-meson peak resulting in a better 
agreement with the data. The ingredients of this model are chiral reduction [49J, 
many-body calculation of the p- meson spectral function |149j . rhosobar exci- 
tations on thermally excited baryon resonances, and a complete assessment of 
mesonic contributions. 

Both the dropping p mass scenario and the in-medium broadening give reason- 
able account of the dielectron enhancement in the 0.3-0.6 GeV/c 2 region. This is 
true also in the case of extension of the acceptance to single track p± > 0.1 GeV/ c 
as apparent from Fig. 15.611 Again, the scenario of in-medium broadening of the 
p spectral function seems more plausible. 

Although the first CERES data for the centrality dependence of the mass 
spectrum was presented several years ago, no systematic theoretical calculations 
are available yet. 

Figure 13.621 compares the invariant mass spectrum for two distinct transverse 
pair momentum selections with model calculations. The free irir scenario without 
in-medium modifications clearly fails to account for the increase of the dielectron 
yield for low transverse momenta. 

The observed transverse-momentum dependency can arise from the fact that 
Lorentz invariance is broken in the thermal frame. Therefore, the in-medium 
propagator, describing the dynamics of the meson, can depend on energy and 
momentum separately. Transverse and longitudinal modes emerge as polarization 
states that are no longer isotropic. Moreover, the thermal occupancy is sensitive 
to a reduction of the p-meson mass: 




(5.49) 




(5.50) 
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Figure 5.61: Comparison of the invariant mass spectrum with theoretical models 
for p± >0.1 GeV/c 162]. The dielectron yield predicted for tttt annihilation was 
added to the standard cocktail of the hadronic sources (without p-meson contri- 
bution). Experimental acceptance as well as momentum resolution were imposed 
on the simulated cocktail data. 

The three-momentum dependence becomes more pronounced for smaller in-medium 
p mass. It leads to a relative enhancement of the p-meson for small three-momenta 
or equivalent small q t . 

Although the model predictions for the dropping p mass and the collision 
broadening scenario differ significantly in this particular representation, the large 
statistical errors of the data do not permit to distinguish both. The apparent 
difficulties of both scenarios to account for the large yield observed at the free 
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Figure 5.62: Comparison of the transverse pair momentum dependence of the 
invariant mass spectrum with theoretical models for p±_ > 0.2 GeV/c (upper panel) 
and p^>0.1GeV/c (lower panel) |162j . The dielectron yield predicted for tttt 
annihilation was added to the standard cocktail of the hadronic sources (without 
p-meson contribution). Experimental acceptance as well as momentum resolution 
were imposed on the simulated cocktail data. 
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Figure 5.63: Comparison of the transverse pair momentum spectra with theo- 
retical models for p_i_ > 0.2 GeV/c |162| . The dielectron yield predicted for tttt 
annihilation was added to the standard cocktail of the hadronic sources (without 
p- meson contribution). Experimental acceptance as well as momentum resolution 
were imposed on the simulated cocktail data. 

vacuum p-peak could also point to an underestimation of the tu-contribution 
which is not well determined so far. 

Finally, a comparison of the transverse momentum spectra for p± > 0.2 GeV/c 
with the model calculations is presented in Fig. 15.631 

To conclude, the present data excludes the most simple scenario with vacuum 
p-mass. The statistical errors, however, do not permit to distinguish the two 
approaches which include in-medium modifications. 

The difference between both models is most evident in the region of the u- 
resonance, namely a factor of two. Excluding one or the other with a confidence 
level of 96% would require a relative statistical error of less than 10(15)%. To 
achieve this level of accuracy for the 3.3 • 10 7 central Pb-Au collisions recorded in 
2000 would require a signal-to-background ratio of larger than 1:5(1:10) assum- 
ing a realistic pair-reconstruction efficiency of 22%. This is a very challenging 
task but the reduction of the systematic errors, i.e. multiplicity and efficiency 
determination, to same level of accuracy might prove even more difficult. 

5.10 Summary and outlook 

The most important result of this new analysis of the dielectron production in 
Pb-Au collisions at 158GeV/c per nucleon is the independent confirmation of 



128 



the previous result: a significant excess of dielectrons observed in the mass range 
of 200 < m ee < 700 MeV/c 2 compared to the expected yield of hadronic sources. 
In the same mass region, there is evidence for a stronger than linear rise of the 
yield with charged-particle multiplicity. The transverse-momentum spectra show 
an enhancement for the transverse pair momentum below 1 GeV/c that increases 
towards small q t . The contribution of very low momentum pairs, i.e. single-track 
momentum of pj_ > lOOMeV/c, was studied for the first time. It was found that 
the excess increases towards small transverse pair momentum. The use of the 
mixed-event technique for background subtractions has ruled out the possibility 
of artifacts of the same-event background subtraction as a source of the dielectron 
enhancement. The refined calibration of SDD and RICH allowed a better rejec- 
tion of the combinatorial background. However, this improvement was partially 
cancelled by a 30% efficiency loss due to new software for SDD-hit reconstruction 
and tracking. 

The comparison with theoretical model calculations shows that the observed 
dielectron yield cannot be explained by the known hadronic sources including 
contributions from free pion annihilation. Indeed, only scenarios invoking in- 
medium modification of the p- and/or 7r-meson can account for the observed 
yield as well as the spectral shape. Both the dropping p mass and the collision 
broadening scenario are viable for the present data, even though the second option 
seems more plausible. 

A precision measurement of the yield and spectral shape of the heavy vec- 
tor mesons is indispensable to distinguish the different scenarios. It requires 
a high statistics data sample combined with a much better momentum resolu- 
tion. The upgrade of the CERES experiment with a radial TPC is expected 
to fulfill both requirements. First studies [HB] achieved a momentum resolution 

of dp/p = \J (0.027) 2 + (0.024 ■ p ■ GeV _1 c) 2 for the reconstruction of A- mesons. 
Further improvements are expected after a refined calibration of the TPC. 

In the fall of 2000, the upgraded CERES experiment was operated with a very 
good performance. A sample of 3.3 • 10 7 central Pb-Au collisions at 158GeV/c 
per nucleon was recorded. The progress made so far promises to fulfill the high 
expectations for a precision measurement of the low-mass dilepton spectrum. 

Other regions of the nuclear matter phase diagram will be explored by ex- 
periments measuring dilepton production in heavy ion collisions, most notable 
HADES [H2] and PHENIX [EE]. Hopefully, the combination of all results will 
soon allow to conclude about possible restoration of chiral symmetry in hot and 
dense nuclear matter and its relation to the QGP phase transition. 
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Appendix A 

New GENESIS event generator 



The GENESIS event generator [HJ EH QUI EH EEHS1 is a tool to simulate the 
relative abundance of dielectrons produced by hadron decays in proton-proton 
(pp), proton- nucleus (pA), and nucleus-nucleus (AA) collisions. The invariant 
mass range covered by the CERES acceptance (m ee < 2 GeV/c 2 ) is dominated 
by the decay of light scalar and vector mesons comprising 7r°, rj, rf, p°, u, and 
(j) jH]. Open charm production is negligible |157j . To create this so-called hadronic 
cocktail, pA and AA collisions are treated as a superposition of individual nucleon- 
nucleon collisions. The hadronic cocktail provides a reference for the comparison 
with the yield observed in pA and AA collisions. Any deviations would indicate 
a violation of the scaling behaviour and/or in- medium effects. 

The simulation requires prior knowledge of the differential production cross 
section, the widths of all decays including dielectrons in the final state, and 
a description of decay kinematics for all relevant particles. Differential cross 
sections are unknown for most light mesons (except 7r°, tj, and T]'). The absolute 
meson yield of pA collisions can be inferred from measurements in pp collisions 
at comparable impact energies (for a compilation see |KH]). 

Proton-nucleus Collisions are modeled by a superposition of nucleon-nucleon 
collisions and the yield thereof is assumed to scale with the mean charged-particle 
multiplicity of a collision system. The relative cross sections (o~/o~ n o) for Pb-Au 
collisions are taken from a thermal model [166 . The model describes parti- 
cle production in heavy ion collisions accurately, as demonstrated in Fig. IA.1I 
The cross sections and branching ratios contained in GENESIS are summarized 
m Table EHJ For comparison with experimental data, the cocktail is divided 
by the total number of charged particles within the nominal detector accep- 
tance. It is directly related to number of produced neutral pions via the ratio 
(N n o/N ch ) = 0A4 Eng. 

The properties of the parent particles are determined by transverse momen- 
tum and rapidity distribution. The kinematic distributions of pions were mea- 
sured by WA98 HSZIEEB!, NA44 [MI!, and NA49 The transverse-momentum 
spectrum of charged pions by NA44 was used to extrapolated the transverse-mass 
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Figure A.l: Comparison of the observed particle ratios with the prediction of the 
thermal model [166 . 



spectrum by WA98 towards small transverse mass including the additional contri- 
bution of the 77 — > 37t° decay. The parameterization is documented in [§U|. The 
inverse slope parameter of the NA44 7r° transverse mass spectrum changes by 
about 5% with centrality. This dependence is not yet implemented in GENESIS 
and, therefore, contributes 5% to the inherent systematic error. WA98 quotes a 
systematic error of 10% for pairs with transverse mass above m± > 400 MeV/c 2 . 



Particle 


Decay 


a/a n o (p+A) 


a/a w o (Pb+Au) 


BR 




e + e~7 


1 


1 


1.198 x 10- 2 


V 


e + e~7 


0.053 


0.085 


5.0 x 10~ 3 


P° 


e + e~ 


0.065 


0.094 


4.44 x 10~ 5 


u 


e + e~7r 


0.065 


0.069 


5.9 x 10~ 5 


UJ 


e + e~ 


0.065 


0.069 


7.15 x lO" 5 




e + e~ 


0.0033 


0.018 


3.09 x 10" 4 


rf 


e + e _ 7 


0.009 


0.0078 


5.6 x 10" 4 



Table A.l: Relative production cross section a/a w o and branching ratio BR of 
light mesons used in new GENESIS jOHj- The relative cross sections for Pb-Au 
collisions are taken from the thermal model [166J. 
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The systematic error of the combined spectrum is estimated to be about 10%. 

The transverse-mass distributions of all other mesons are described by ex- 
ponential distributions. The inverse slope parameter T increases systematically 
with particle mass (see Fig. IA.2jl as a result of collective flow. It is parameterized 
as: 

T = 0.175 GeV + 0.115 ■ m (c = 1) . (A.l) 
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The rapidity distributions are assumed to resemble those of the negative 
hadrons measured by NA49 |170j . The width is adjusted according to the maxi- 
mum kinematic rapidity limit. 

The cross section is factorized for calculation of the final decay kinematics 
into a contribution of a QED point source and an form factor F describing the 
internal structure of the meson: 



da 
dq 2 




F{q 



2\2 



(A.2) 



point source 



The shape of the form factors are based on the measurements of the Lepton-G 
collaboration and theoratical model calculations. An extensive summary of all 
relevant form factors is presented in [OH]. According to the detailed discussion of 
the systematic errors presented in j52j, the uncertainty in the branching ratios 
and the form factors contribute about 15% below 450MeV/c 2 , 30% in the mass 
range of 450-750 Me V/c 2 , and 6% above 750MeV/c 2 . 
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Finally, the dielectrons generated are subject to detector acceptance and finite 
momentum resolution. 

A comparison of the invariant mass spectra (see Fig. 12.51 in Sec. 12.2ft shows 
the hadronic cocktail to describe the observed yield in p-Be and p-Au collisions 
well within the systematic errors of about 20%. As a result of a recent review and 
extension of the GENESIS code [HD], the cocktail has decreased by up to 20% in 
the mass region of 0.15 < m cc < 1.5 GeV/c 2 (see Fig. IA.3|) . 

The low-mass yield (m ee < 0.2 GeV/c 2 ) was previously underestimated by 35% 
as a result of a problem in the procedure used to apply the momentum resolution. 
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Figure A. 3: 

Impact of the corrected 
GENESIS acceptance filter 
on the hadronic cocktail. 
The previously used cock- 
tail |113j underestimates the 
yield in the low-mass region 
by 35%. 



invariant pair mass 



The cocktail for p-Be and p-Au collisions is much less affected (5%) due to the 
lower transverse momentum cut (p± < 50 MeV/c). With hindsight to this fact, 
the validity of the statement about the total exhaustion of the observed dilepton 
yield by the known hadronic sources remains sustained. 

The new hadronic cocktail for Pb-Au collisions at 158GeV/c per nucleon is 
plotted in Fig. IA.41 Integration of spectrum gives a total yield of (iV ee /A ir ch ) acc = 
1.21 • 10~ 5 which divides into (A^ ee /A^ ch ) acc . = 1.01 • 10~ 5 for m cc < 0.2 GeV/c 2 
and (iVee/iV c h)acc. = 2.04- 10" 6 for m ce > 0.2 GeV/c 2 . In contrast to [113], the 
yield of pairs with very high mass (m ee > 1.5 GeV/c 2 ) is slightly reduced as result 
of the upper momentum cut of p<9GeV/c imposed for the discrimination of 
high-momentum pions. 
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Figure A. 4: Genesis cocktail of the known hadronic sources for Pb-Au colli- 
sions at 158GeV/c per nucleon. Integration of spectrum gives a total yield of 
(A ee /A ch ) acc . = 1.21 • 1(T 5 which divides into (A ec /A ch )acc. = 101 • 10" 5 for 
m ee < 0.2 GeV/c 2 and (A cc /A ch ) acc . = 2.04- 1 0~ 6 for m cc > 0.2 GeV/c 2 . The ex per- 

imental momentum resolution of dp/p = \J (0.041) 2 + (0.022 • p ■ GeV _1 c) 2 was 
applied. 
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Appendix B 

Summary of rejection strategy 



Description 


Cut selection 


Comment 


acceptance 


SDD-1/2 0=141.-299. mrad 
RICH-1 0=141.-299. mrad 
RICH-2 0=141.-258. mrad 


2/3 ring accepantce 
2/3 ring accepantce 


double dE/dx 
rejection 


SDD-1 dE/dx < 1133 OR 
SDD-2 dE/dx < 1204 


5 mrad resummation 


SDD-1 

conversion 


SDD-2 dE/dx < 1204 OR 
RICH-1 sum amplitude < 1550 


8 mrad resummation 


track quality 


SDD-1-2 fi match < 1.3-1.5 mrad 
SDD-RICH-1 Q match <2 ( 7 match (p) 

RICH-1— 2 ma tch < 2 (Tmatch(p) 

RICH-2 Hough-2 amplitude > 360 
RICH-2 xVKts - 1) < 1.3 
A0 R icH-i-RicH-2 < 300 mrad 
A0sdd-pd > 5 mrad 


target disk dependence 
a = ^/l.7 2 + 1.0 2 /p 2 
a= ^/l.5 2 + 1.2 2 /p 2 
ring candidate quality 
ring fit quality 
max. B-field deflection 
min. B-field deflection 


high-px pion 
rejection 


ifelCH-l 

14.62 16.20 
A0sdd-pd > 17 mrad 


pressure dependence 
of nominal radius 


final acceptance 


77=2.1-2.65 (single track) 

p± < 0.2(0.1) GeV/c (single track) 

a ee > 35 mrad 


eq. 0=141.-243. mrad 
pair opening angle 



Table B.l: Summary of rejection cuts. The rejection cuts were applied in sequence 
as presented in the table. The order chosen was motivated by the intention to 
apply the most powerful and reliable cuts first. Most cuts are strongly correlated 
and cannot be looked at individually. All individual rejection conditions must be 
fulfilled according to a logical AND operation unless otherwise stated. 
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